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To  my  mother 


ABSTRACT 


The  time  dependent  lateral  deformations  of  the  clay 
foundation  of  an  oilsands  mine  tailings  dyke  are  analysed  in 
this  thesis.  The  lateral  deformations  analysed  were  recorded 
by  several  inclinometers  over  a  period  of  approximately  four 
years.  Based  on  the  nature  of  the  lateral  deformations  and 
piezometer  readings  in  the  clay  foundation  it  is  postulated 
that  the  movements  are  a  result  of  a  drained  creep 
mechanism. 

On  the  basis  of  laboratory  tests  the  foundation  clay  is 
normally  consolidated,  stiff,  medium  plastic  and  contains 
approximately  two  percent  organic  matter  and  four  percent 
calcium  carbonate. A  series  of  oedometer  tests  demonstrate 
that  the  clay  has  a  'low'  secondary  compressibility 
according  to  Mesr i ’ s ( 1 973 )  classification. 

Incremental  drained  creep  triaxial  tests  were  performed 
on  samples  with  the  bedding  plane  horizontal  and  at  an  angle 
of  45°. It  is  shown  that  the  creep  behaviour  of  the  clay  is 
isotropic  and  insensitive  to  normal  stress .Furthermore , it  is 
demonstrated  that  the  axial  creep  strain  rate  , after  all 
excess  pore  pressures  have  dissipated,  follows  a 
relationship  of  the  Singh-Mi tchell  (1968)  type. 

An  analysis  of  one  set  of  inclinometer  readings  is 
presented  which  illustrates  a  method  whereby  the  creep 
parameters  in  the  Singh~Mi tchell  (1968)  equation  can  be 
derived  from  field  shear  strain  versus  time  curves.  The  time 
from  the  actual  initiation  of  creep  to  the  time  at  which  the 


v 


inclinometer  was  installed  ,t,,is  identified  as  a  key 
parameter  in  the  field  creep  relationship. 
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1.  Introduction 


1 . 1  General 

Surface  mining  of  the  vast  reserves  of  the  Athabasca 
oilsands  deposit  began  in  1967  on  a  large  scale  when  the 
Suncor  (then  G.C.O.S.  -  Great  Canadian  Oil  Sands)  oilsands 
plant  began  production.  In  1978  a  second  plant,  Syncrude, 
began  production  and  others  are  planned  for  construction  in 
the  1980’s.  The  location  of  the  Athabasca  oilsands  deposit 
and  the  f orement ioned  plants  are  shown  on  Figure  1.1. 

Extraction  of  the  bitumen  from  the  sand  by  the  Clark 
hot  water  process  creates  large  volumes  of  tailings 
consisting  of  sand,  silt,  clay,  water  and  small  quantities 
of  bitumen.  For  an  oilsands  plant  producing  100,000  barrels 
of  oil  per  day  approximately  150,000  cubic  meters  of 
tailings  are  produced  daily. 

The  disposal  of  these  tailings  creates  the  need  for 
large  capacity  tailings  ponds.  Tailings  cannot  be-  disposed 
of  in  the  open  pit  mine  during  the  first  years  of  production 
so  a  separate  tailings  structure  has  traditionally  been 
built  close  to  the  plant  site.  The  tailings  structures  at 
the  Suncor  and  Syncrude  plants  have  been  constructed 
according  to  the  upstream  method  (see  Mittal  and 
Morgenstern,  1977).  Compacted  tailings  sand  comprises  the 
supporting  shell  and  the  remainder  of  the  tailings  forms  the 
beach  and  pond  of  the  tailings  disposal  structure. 
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Figure  1.1  Location  of  major  oilsands  deposits  in  Canada 
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The  tailings  dykes  at  the  Suncor  and  Syncrude  plants 
are  raised  in  stages  usually  with  a  freeboard  of  one  year  in 
tailings  output.  The  stage  construction  of  these  structures 
allows  the  performance  of  the  supporting  shells  to  be 
monitored  with  instrumentation.  Design  can  be  optimized 
through  rational  interpretation  of  the  structure's 
performance.  Such  a  policy  of  performance  monitoring  has 
been  pursued  at  Suncor fs  Tar  Island  tailings  dyke  and 
provided  the  impetus  for  this  thesis. 

In  1975  at  the  request  of  a  Design  Review  Panel, 
convened  by  Alberta  Environment  to  review  the  feasibility  of 
the  continued  raising  of  Tar  Island  Dyke,  several 
inclinometers  were  installed  to  monitor  the  lateral 
movements  of  the  clay  foundation.  The  time  dependent  lateral 
movements  displayed  by  these  inclinometers  are  the  focus  of 
this  thesis.  As  the  inclinometer  readings  were  accompanied 
by  a  plethora  of  information  on  foundation  geology, 
geotechnical  properties,  pore  pressures  and  construction 
history  this  case  history  represented  a  unique  opportunity 
for  the  investigation  of  time-dependent  lateral 
displacements  beneath  embankments. 


1.2  Lateral  Displacements  Beneath  Embankments 

It  is  of  value  to  briefly  review  the  possible  stress 
paths  in  soft  clay  foundations  in  order  to  cast  the  lateral 
deformations  beneath  Tar  Island  Dyke  in  a  conceptual 
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framework.  The  prediction  and  analysis  of  lateral 
displacements  in  soft  foundations  beneath  embankments  has 
been  afforded  much  less  attention  in  the  literature  than 
vertical  displacements.  Perhaps  this  is  because  lateral 
displacements  are  much  harder  to  predict  and  recorded  less 
often  than  vertical  displacements.  Poulos  (1972)  suggests  a 
number  of  reasons  for  their  poor  prediction.  Recently 
Tavenas  et.al  (1979)  and,  earlier,  Rutledge  and  Gould  (1973) 
documented  numerous  case  histories  concerning  lateral 
displacements  beneath  embankments. 

Tavenas  et.al  (1979)  presented  some  simple  empirical 
methods  for  predicting  lateral  displacements  of  normally 
consolidated  clay  foundations.  Figure  1.2  is  taken  from 
their  paper  and  shows  possible  stress  paths  in  clay 
foundations  during  construction  of  embankments.  It  is  of 
value  to  follow  the  stress  paths  shown  on  Figure  1.2  with  a 
view  to  understanding  the  mechanisms  causing  lateral 
deformations . 

From  O'  to  P'  the  clay  consolidates  and  the  stresses 
move  towards  the  limit  state  or  yield  surface.  Lateral 
deformations  in  this  region  of  overconsolidation  are  usually 
rapid  because  of  high  values  of  Cv.  When  the  stress  path 
intercepts  the  yield  surface  all  deformations  become 
undrained.  The  stress  path  will  follow  the  yield  surface  to 
failure  as  shown  by  Tavenas  et.al  ( op.cit .)  unless 
construction  is  halted. 
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Figure  1.2  Total  and  effective  stress  paths  under 
embankment  during  construction  (after  Tavenas  et.al ,  1979) 
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The  writer  has  added  stress  path  P’-Q’  to  Figure  1.2. 
The  foundation  clay  will  follow  this  stress  path  if  it  is 
allowed  to  consolidate  before  reaching  failure.  It  should  be 
noted  that  creep  movements  can  occur  along  any  of  the 
previously  mentioned  stress  paths.  The  importance  of  creep 
to  lateral  movements  depends  upon  the  geometry  of  the 
embankment  and  foundation,  the  creep  susceptibility  of  the 
clay  and  the  construction  sequence  of  the  embankment. 

Piezometers  in  the  foundation  clay  of  Tar  Island  Dyke 
are  recording  either  dissipation  of  excess  pore  pressures 
with  time  or  steady  state  pore  pressure  conditions.  This 
indicates  that  the  stress  path  for  the  foundation  clay  is  on 
a  line  similar  to  P’-Q'  in  Figure  1.2.  It  therefore  follows 
that  the  lateral  deformations  are  a  result  of  consolidation, 
creep  or  some  combination  of  both. 


1.3  Objective  and  Scope 

The  objective  of  this  thesis  is  the  understanding  of 
the  time  dependent  lateral  displacements  beneath  Tar  Island 
Dyke . 

Chapter  2  presents  the  construction  history,  site 
conditions  and  selected  instrumentation  readings  for  the 
clay  foundation.  Chapter  3  is  a  review  of  the  available 
literature  on  the  creep  behaviour  of  soils.  The  laboratory 
program  to  investigate  the  creep  behaviour  of  the  foundation 
clay  is  detailed  in  Chapters  4  and  5. 
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Chapter  6  examines  the  characteristics  of  the  lateral 
displacement  versus  depth  and  time  curves  from  the 
inclinometers.  Chapter  7  outlines  the  overall  conclusions  of 
the  investigation  and  presents  recommendations  for  future 


research . 


' 


2.  Tar  Island  Dyke  -  Site  Conditions  and  Instrumentation 


2 . 1  Introduction 

The  purpose  of  this  chapter  is  to  present  a  brief 
summary  of  the  construction  history  of  Tar  Island  Dyke  and 
selected  data  concerning  foundation  conditions  and 
instrumentation.  Most  of  the  material  regarding  construction 
history  is  taken  from  Mittal  and  Hardy  (1977)  and  the  Design 
Review  Panel  (1977).  The  information  regarding  foundation 
conditions  and  instrumentation  has  been  provided  by  the 
designers  of  Tar  Island  Dyke,  Hardy  Associates  Ltd  (1978). 
The  field  program  to  secure  soil  samples  for  laboratory 
testing  is  also  recounted  in  this  chapter. 


2.2  Construction  History 

Tar  Island  Dyke  was  originally  conceived  as  a  low,  12  m 
high  conventional  earthfill  dyke  of  waste  overburden  with 
sufficient  volume  to  store  tailings  until  exploited  portions 
of  the  mine  became  available  for  disposal  of  tailings. 
Construction  began  on  this  low  dyke  in  1965  and  continued 
intermittently  until  1967  when  it  was  apparent  that  tailings 
volumes  were  greater  than  anticipated  and  that  the  continued 
construction  of  a  conventional  earthfill  dyke  to  store  the 
tailings  would  be  prohibitively  expensive.  It  should  be 
stated  that  during  these  early  pioneering  days  even 
essential  tailings  dam  design  parameters  such  as  the  bulking 
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factor  could,  at  best,  only  be  estimated. 

In  November  of  1967  application  was  made  to  the  Energy 
Resources  Conservation  Board  (ERCB)  of  Alberta  to  raise  the 
dyke  with  hydraulically  placed  sand  tailings  instead  of 
continuing  with  the  conventional  earthfill  dyke.  The 
application  was  approved  and  since  1968  the  dyke  has  been 
built  with  hydraulically  placed  sand  tailings  according  to 
the  upstream  method  of  construction.  Aerial  photographs  of 
Tar  Island  Dyke  at  different  stages  of  construction  are 
shown  in  Plate  2.1.  It  is  of  benefit  to  briefly  summarize 
the  construction  procedure  for  the  dyke. 

The  dyke  consists  of  a  compacted  downstream  zone  and  an 
uncompacted  beach  zone  which  retain  a  pond  as  shown  in 
Figure  2.1.  The  compacted  zone  is  constructed  by  spigotting 
from  a  tailings  pipeline  into  a  cell  which  is  enclosed  by  a 
peripheral  berm,  1.5  m  to  2.0  m  high.  A  cell  typically  has 
dimensions  of  30  m  to  90  m  by  300  m  to  460  m.  Attached  to 
the  spigot  point  on  the  tailings  pipeline  is  a  deflector 
spoon  which  has  the  effect  of  spraying  the  exiting  tailings 
and  thereby  inhibits  the  formation  of  a  stilling  basin 
beneath  the  spigot  point.  Consequently , the  tailings  velocity 
in  the  cell  is  increased  which  prevents  premature 
sedimentation  of  fines  in  the  cell.  Sand  settles  from 
suspension  in  the  tailings  stream  in  the  cell  and  is 
subsequently  distributed  within  the  cell  by  a  caterpillar 
tractor.  The  action  of  the  tractor  is  sufficient  to  increase 
the  relative  density  of  the  sand  in  the  cell  to  above  75 
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Figure  2.1  Tar  Island  Dyke,  section  at  Station  65+00 
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percent.  The  tailings  stream  with  sand  and  fines  in 
suspension  exits  from  one  corner  of  the  cell  through  a 
spillway  box.  The  suspended  sand  and  fines  sediment  from  the 
tailings  stream  and  form  the  beach  which  has  an  average 
slope  of  approximately  8  percent  to  the  center  of  the  pond. 
A  layer  of  sludge,  composed  of  fines  and  oily  water, 
overlies  the  beach  sand  in  the  tailings  pond.  Sufficient 
storage  capacity  is  created  each  year  by  constructing  the 
compacted  section  of  the  dyke  for  only  part  of  the  year. 
During  the  rest  of  the  year  tailings  are  over-boarded 
directly  into  the  pond. 

The  f orement ioned  method  of  construction  has  been  very 
successful  with  the  exception  of  a  few  minor  downstream 
liquefaction  failures  of  the  compacted  zone  in  1972,  1973 
and  1974.  These  failures  have  since  been  avoided  with  local 
improvements  in  construction  procedure.  The  stability  of  the 
Tar  Island  Dyke  is  discussed  in  detail  in  the  Design  Review 
Panel  (1977)  report  and  will  not  be  discussed  in  this 
thesis . 

The  downstream  face  of  the  dyke  has  an  overall  slope  of 
three  horizontal  to  one  vertical.  There  are  berms  for  access 
at  approximately  13  m  vertical  intervals.  Seepage  through 
the  dyke  is  controlled  by  blanket  drains  composed  of  plant 
coke.  All  seepage  collected  in  the  drains  is  now  returned  to 
the  pond.  In  1979  the  dyke  had  reached  a  height  of  89  m  and 
ultimately  it  is  planned  to  raise  it  to  97.5m. 
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2.3  Foundation  Characteristics 

2.3.1  Geology 

Tar  Island  Dyke  rests  on  the  active  floodplain  of  the 
Athabasca  River .  The  abutments  of  the  dyke  are  the  Athabasca 
River  valley  walls  which  are  composed  of  Clearwater 
Formation  shale  and  McMurray  Formation  oilsand.  The  common 
stratigraphic  sequence  in  the  general  area  consists  of 
Recent  deposits  overlying  Pleistocene  glacial  deposits  which 
rest  unconf ormably  on  the  Cretaceous  Clearwater  and  McMurray 
Formations  which,  in  turn,  unconformbly  overlie  a  Devonian 
limestone.  The  meandering  incised  Athabasca  River  has 
completely  eroded  all  Pleistocene  and  Cretaceous  sediments 
in  the  site  specific  area.  The  foundation  stratigraphy 
consists  of  muskeg  overlying  a  Recent  deposit  of  silt  and 
clay  overlying  sand  which,  in  turn,  rests  on  the  Devonian 
limestone.  Cobbly  gravel  is  usual  just  above  the 
sand/limestone  contact.  The  clay  deposit,  which  will 
occasionally  be  referred  to  as  Tar  Island  clay, should  not  be 
confused  with  the  basal  clay  shales  of  the  McMurray 
Formation  (see  Dusseault  and  Scafe,  1979). 

The  original  foundation  investigations  for  the  dyke 
were  undertaken  in  1964  and  1965.  Since  then  there  have  been 
numerous  other  investigations  (see  summary  in  Design  Review 
Panel,  1977).  An  isopach  map  of  the  silt  and  clay  deposit 
based  on  the  1964  and  1965  borehole  logs  is  given  in  Figure 
2.2.  The  thicknesses  of  silt  and  clay  encountered  in 
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Athabasca  River 
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Figure  2.2  Isopach  map  of  silt  and  clay  deposit 
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selected  boreholes  drilled  in  1975  are  also  shown  on  the 
isopach  map. 

The  silt  and  clay  deposit  originated  by  sedimentation 
in  an  abandoned  meander  of  the  Athabasca  River.  The  basal 
sands  represent  an  Athabasca  River  bed  load  sediment  which 
was  deposited  prior  to  complete  cutoff  of  the  meander.  The 


major  source  of  sediment 

to 

the  meander  was  undoubtedly 

suspended 

silt  and  clay 

from 

overbank 

flows  when 

the 

Athabasca 

was  in  flood. 

This 

genesis 

implies  that 

the 

deposit  should  be  coarser  at  its  margins.  There  is  also 
evidence  of  essentially  still  water  deposition  in  the  form 
of  numerous  calcium  carbonate  shells  throughout  the  deposit. 

2.3.2  Geotechnical  Properties  of  Silt  and  Clay  Deposit 

The  deposit  consists  of  interbedded  silt  and  clay  with 
occasional  sand  lenses.  Silt  and  sand  predominate  towards 
the  margins  of  the  deposit.  From  34  grain  size  analyses  on 
1964  and  1965  samples  the  deposit  averages  20  percent  sand, 
63  percent  silt  and  17  percent  clay.  The  deposit  also 
contains  varying  amounts  of  coal  fragments,  calcium 
carbonate  shell  material  and  organic  material. 

Before  imposition  of  the  dyke  the  deposit  was  medium 
stiff  with  undrained  strengths  ranging  from  45  kPa  to  90 
kPa.  The  undrained  strength  over  effective  stress  ratio 
(Cu/p' )  was  0.45.  Of  course,  the  undrained  strength  has 
increased  with  time  as  the  deposit  has  consolidated  under 
the  weight  of  the  dyke. 


:  I 

* 


16 


The  liquid  limit  of  the  deposit  varies  from  24  percent 
to  65  percent  and  the  plasticity  index  varies  from  4  percent 
to  37  percent.  The  Atterberg  limit  test  results  from  the 
early  investigations  are  presented  in  Figure  2.3.  The 
plasticity  is  greatest  where  the  clay  is  thickest.  The 
compression  index ,Cc , from  the  results  of  16  oedometer  tests 
averages  0.28.  The  results  of  effective  shear  strength  tests 
performed  by  Hardy  Associates  Ltd.  (1978)  are  presented  in 
Chapter  5. 


2.4  Field  Sampling  Program 

A  sampling  program  was  undertaken  in  the  fall  of  1979 
to  obtain  samples  of  the  foundation  clay.  The  program  was 
coordinated  with  the  annual  installation  of  instrumentation 
by  the  dyke  designers  , Hardy  Assoc iates ( 1 978 )  Ltd.,  in  the 
critical  section  of  the  dyke  between  Stations  56+00  and 
65+00.  A  Mayhew  1000  wet  rotary  rig  mounted  on  a  Nodwell 
tracked  vehicle  was  used  for  all  sampling. 

The  main  requirements  for  the  samples  were  that  they  be 
representative  of  the  clay  in  the  critical  section  and, 
furthermore,  that  they  be  as  undisturbed  as  practicable.  To 
this  end,  it  was  attempted  initially  to  obtain  samples  from 
beneath  the  303  m  berm  using  10.2  cm-0  shelby  tubes. 
Unfortunately,  continual  sloughing  of  the  hole  within  the 
tailings  sand  precluded  sampling  at  the  chosen  location  on 
the  303  m  berm.  The  instrumentation  program  dictated  that 
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Figure  2.3  Atterberg  limit  test  results 
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the  rig  be  moved  to  the  291  m  berm  so  samples  were 
eventually  obtained  there.  The  inside  diameter  of  the  shelby 
tube  sampler  had  to  be  decreased  from  10.2  to  7.3  cm  as  the 
rig  did  not  have  sufficient  capacity  to  pull  or  drive  the 
larger  shelby  tubes. 

The  test  hole  log,  S79-109,  at  which  the  samples  were 
obtained  is  included  in  Appendix  A.  The  location  of  the  test 
hole  is  shown  on  Figure  2.4.  The  description  of  the  samples 
is  given  in  Chapter  5. 


2.5  Foundation  Instrumentation 


2.5.1  General 

Tar  Island  Dyke  is  heavily  instrumented  with 
piezometers  to  monitor  the  seepage  regime  in  the  tailings 
sand  and  in  the  foundation.  Foundation  piezometers  are 
concentrated  at  Stations  56+00  and  65+00  where  the  clay 
deposit  is  thickest.  Numerous  inclinometers  to  monitor 
lateral  deformation  of  the  clay  deposit  are  also  located  at 
these  two  stations.  Vertical  deformation  of  the  clay  deposit 
is  not  monitored  directly  but  is  routinely  estimated  from 
water  content  change  and  thickness  change  as  determined  by 
aperiodic  drilling  programs. 

Because  of  the  abundance  of  instrumentation  data,  only 
the  readings  from  selected  instruments  will  be  reported 
here.  Figure  2.4  shows  the  locations  of  selected 
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instrumentation  in  plan  while  Figure  2.1  presents  the  same 
instrumentation  in  section. 

2.5.2  Foundation  Piezometers 

The  readings  from  the  foundation  piezometers  at  Station 
65+00  from  the  303  m  and  the  291  m  berm  are  illustrated  in 
Figure  2.5.  The  readings  from  the  277  m  and  264  m  berms  are 
shown  in  Figure  2.6  On  every  pore  pressure  plot  a  ’steady 
state'  line  is  shown.  This  line  represents  the  estimated 
pore  pressure  regime  in  the  clay*  layer  at  the  end  of 
consolidation.  The  pore  pressure  at  the  base  of  the  clay 
layer  is  equal  to  that  in  the  basal  sand  whichcis  in  direct 
communication  with  the  Athabasca  River.  The  pore  pressure  at 
the  top  of  the  clay  is  governed  by  the  seepage  regime  in  the 
tailings  sand  and  has  been  determined  from  piezometer 
readings  or  estimated. 

It  can  be  seen  from  Figure  2.5  that  the  clay  deposit 
beneath  the  303  m  berm  has  excess  pore  pressures  and  is 
consolidating.  Recently  piezometers  have  been  installed 
upstream  of  the  303  m  berm  and  these  instruments  indicate 
that  the  clay  deposit  there  has  excess  pore  pressures  and  is 
consolidating.  From  Figures  2.6  and  2.7  it  can  be  seen  that 
the  clay  deposit  from  the  291  m  berm  downstream  is 
essentially  fully  consolidated.  It  follows  that  stresses  on 
the  foundation  caused  by  recent  dyke  construction  are  not 
significant  enough  to  induce  pore  pressures  in  the 
foundation  clay  downstream  of  the  291  m  berm. 
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Figure  2.4  Tar  Island  Dyke,  location  of  instrumentation 
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O  December,  1975 

x  December,  1976  AP  -  Pneumatic 

+  December,  1977  P  -  Standpipe 

□  December,  1978 
•  April,  1980 


Pore  Pressure,  kPa 


Figure  2.5  Summary  of,  foundation  pore  pressures,  upper-303 
berm  and  lower-291  m  berm 
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Figure  2.6  Summary  of  foundation  pore  pressures, 
berm  and  lower-264  m  berm 
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2.5.3  Inclinometers 

A  number  of  inclinometers  have  been  installed  in  the 
dyke  since  1976  to  monitor  lateral  movements  of  the 
foundation.  The  Sinco  Digitilt  inclinometer  has  been  used 
for  all  slope  indicator  installations  at  Tar  Island  Dyke 
(see  Savigny,  1980  for  description  of  inclinometer 
charac ter i st ics ) .  The  readings  from  five  inclinometers, 
numbered  S76-101,  S78-101sf  S76- 1 02 , S76- 1 03  and  S76-104, 
have  been  selected  for  review.  Their  plan  location  is  shown 
on  Figure  2.4  while  their  position  in  section  is  shown  on 
Figure  2.1. 

The  bases  of  all  inclinometers  are  founded  in  the 
limestone  to  provide  a  fixed  reference  point  for  calculation 
of  lateral  movements.  Inclinometer  readings  covering  the 
period  September,  1976  to  May,  1980  were  made  available  to 
the  writer  by  Hardy  Associates  (1978)  Ltd.  on  magnetic  tape. 
Hardy  Associates  (1978)  Ltd.  also  made  their  inclinometer 
data  reduction  computer  program  available.  Several 
modifications  were  made  to  the  program  and  it  is  of  value  to 
briefly  review  the  important  data  reduction  steps. 

As  a  first  step,  the  average  of  the  A  component  and  the 
B  component  readings  were  multiplied  by  an  instrument 
constant  to  give  the  horizontal  displacement  in  each  of  the 
A  and  B  directions.  The  vector  sum  of  the  A  and  B 
displacements  was  then  computed  to  give  the  maximum 
horizontal  displacement,  irrespective  of  orientation,  at 
each  depth.  The  tacit  assumption  was  made  in  this 
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calculation  step 

that 

the 

maximum  displacement 

is  in 

the 

normal  section 

to 

the 

longitudinal  axis  of  the  dyke. 

The 

displacements  were 

then 

cumulated  from  the 

fixed 

base 

reference  point 

in 

the 

standard  manner.  The 

program 

was 

modified  to  incorporate  computer  graphics  and  a  plot  of 
cumulative  horizontal  displacement  with  depth  for 
inclinometer  S76-101  is  illustrated  in  Figure  2.7. 

The  shear  strain  was  calculated  in  a  similar  manner  and 
is  given  by: 

=  Arcsin (A*C) *  1 00%  2.1 

where  %  =  shear  strain 

A  =  vector  sum  of  average  A  and  average  B 

component  readings 
C  =  instrument  constant 

A  plot  of  shear  strain  versus  depth  for  inclinometer  S76-101 
is  given  in  Figure  2.7.  Similar  plots  to  those  presented  for 
inclinometer  S76-101  are  included  in  Appendix  A  for  the 
other  inclinometers  enumerated  previously. 


2.6  Nature  of  Foundation  Movements 

The  inclinometer  plots  illustrate  that  the  clay  deposit 
is  experiencing  time  dependent  shear  strain.  The  basal  sand 
is  also  exhibiting  shear  strain  with  time.  However,  the 
tailings  sand  above  the  clay  deposit  does  not  show  any  shear 
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Figure  2.7  Inclinometer  S76- 1 0 1 , calibrated  January  26,1976 
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straining.  The  dyke  appears  to  be  moving  essentially  as  a 
rigid  body  on  its  foundation. 

The  time  dependence  of  the  shear  strain  in  the 
foundation  materials  may  either  be  a  consequence  of  creep  or 
of  lateral  consolidation.  It  is  considered  that  creep  is  the 
dominant  mechanism  causing  the  time  dependent  shear  strains 
since : 

1.  Shear  straining  is  occurring  in  portions  of  the 
foundation  where  any  excess  pore  pressures  have  already 
dissipated.  This  is  apparent  in  the  top  and  bottom 
portions  of  the  clay  deposit  and  in  the  full  thickness 

°  of  the  clay  deposit  downstream  of  the  291  m  berm. 

2.  A  mechanism  of  lateral  consolidation  necessitates 
compression  in  the  horizontal  direction.  However  lateral 
strains  are  in  the  order  of  0.04  percent  which  is 
insignificant  compared  to  the  shear  strains  of  1  percent 
to  3  percent  over  the  same  time  period. 

As  a  consequence  of  these  observations  the  remainder  of 
the  thesis  is  directed  towards  defining  the  creep  behaviour 
of  the  clay  deposit. 


. 


3.  A  Review  of  the  Creep  Behaviour  of  Soils 


3 . 1  Introduction 

The  creep  deformations  of  soils  are  an  important 
consideration  in  a  wide  variety  of  geotechnical  problems. 


These  problems 

range 

from 

the  time  dependent 

settlements 

of 

foundations , 

after 

all 

excess  pore 

pressures 

have 

dissipated,  to 

the 

time 

dependent  deformation  of 

soft 

embankment  foundations  which  may  eventually  fail  in  creep 
rupture  before  any  excess  pore  pressures  have  dissipated. 
The  estimation  and  importance  of  creep  deformation  in  design 
has  been  the  subject  of  considerable  research,  especially 
over  the  last  twenty  years.  Most  research  has  been  directed 
towards  understanding  the  creep  behaviour  of  soils  in  the 
laboratory.  Regretably,  there  are  few  case  histories  of 
field  behaviour,  with  the  exception  of  long  term  settlement 
of  foundations,  which  link  a  rational  design  approach  to 
actual  field  behaviour.  Of  necessity  then,  this  review  of 
creep  behaviour  is  restricted  to  a  review  of  laboratory 
investigations. 

Creep  is  defined  here  as  time  dependent  deformation 
under  a  sustained  change  in  stresses,  exclusive  of 
hydrodynamic  effects.  Creep  in  soils  depends  upon  a  host  of 
factors  which  include  time,  temperature,  soil  type,  soil 
structure,  stress  history,  stress  state  and  drainage 
conditions.  Figure  3.1  illustrates  the  wide  variation  of 
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creep  response  among  different  soil  types. 

There  are  two  approaches  to  the  investigation  of  the 
creep  behaviour  of  soils  according  to  Ladanyi  (1972).  The 
first  is  a  fundamental  approach,  which  he  termed  the 
micromechanistic  approach,  in  which  creep  behaviour  is 
related  to  events  occurring  at  the  particle  level.  The 
second  is  a  phenomenological  approach,  which  he  termed  the 
microanalyt ical  approach,  in  which  creep  behaviour  is 
related  to  macroscopic  experimental  findings.  A  brief 
outline  of  the  principal  work  in  soil  mechanics  in  relation 
to  the  fundamental  approach  is  presented  in  section  3.2. 
This  outline  has  been  included  because  it  adds  to  our 
understanding  of  soil  behaviour,  but  kept  brief  because  it 
has  not  yet  been  successful  in  providing  constitutive 
relationships  that  can  be  used  in  practice. 

The  phenomenological  approach  is  presented  in  section 
3.3  in  more  detail  than  the  fundamental  approach  as 
constitutive  relationships,  albeit  empirical,  have  evolved 
from  this  approach  which  are  used  in  practice.  The  creep 
rupture  strength  of  soils  as  determined  from  laboratory 
creep  tests  is  discussed  in  Section  3.4. 


3.2  The  Fundamental  Approach 


■ 
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3.2.1  General 

There  have  been  a  number  of  researchers  who  have 
developed  theories  to  explain  and  predict  the  creep 
behaviour  of  soils  on  a  particle  interaction  level.  Most  of 
these  theories  can  be  grouped  into  three  categories.  The 
first  category  is  that  based  upon  the  viscosity  of  the 
adsorbed  water  layers  around  clay  particles  which  was  first 
advocated  by  Terzaghi  (1941)  and  Taylor  (1942).  The  second 
is  rate  process  theory  which  was  popularized  in  soil 
mechanics  in  the  1960’s  by  a  number  of  researchers  including 
Christensen  and  Wu  (1964)  and  Mitchell  et.al  (1968).  The 
third  category  is  that  originated  by  de  Jong  (1968)  which  is 
based  upon  considerations  of  a  cavity  channel  network  model 
for  the  soil  pores.  A  brief  review  of  the  main  elements  of 
each  of  these  theories  will  now  be  given. 

3.2.2  Considerations  Based  on  Viscosity  of  Adsorbed  Water 

Clay  particles  are  surrounded  by  a  layer  of  adsorbed 
water  that  has  a  bonding  and  structure  which  is  different 
than  the  pore  water  (see  Mitchell,  1976).  The  presence  of 
these  adsorbed  water  layers  was  known  in  the  early  1940' s 
and  attempts  were  made  to  explain  secondary  consolidation  in 
terms  of  their  viscosity.  Taylor  (1942)  was  one  of  the  first 
researchers  in  soil  mechanics  to  hypothesize  that  the 
viscosity  of  the  adsorbed  water  layers  were  of  substantially 
higher  magnitude  than  that  of  the  pore  water  and, 
furthermore,  that  the  adsorbed  water  was  responsible  for  the 
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viscous  component  in  the  effective  stress-strain-time 
behaviour . 

Terzaghi  (1941)  postulated  that,  after  primary 
consolidation,  the  load  carried  by  the  soil  was  somehow 
distributed  between  grain  to  grain  contacts  and  adsorbed 
water  bonds.  Secondary  consolidation  was  thought  to  be  the 
time  dependent  transfer  of  load  from  the  adsorbed  water  to 
the  grain  to  grain  contacts.  He  also  held  that  the  viscosity 
of  the  adsorbed  water  bonds  increased  as  the  particles  moved 
together  and,  therefore,  the  load  transfer  occurred  at  an 
ever  decreasing  rate. 

Walker  (1969a)  used  this  hypothesis  to  explain  the 
independence  of  volume  creep  rate  on  stress  level.  He 
contended  that  the  increase  in  the  structural  viscosity  of 
the  adsorbed  water  bonds  was  counteracted  exactly  by  the 
increased  shear  stress. 

One  method  of  validating  the  influence  of  the  viscosity 
of  the  adsorbed  water  bonds  on  secondary  consolidation  is  to 
replace  the  adsorbed  water  and  pore  water  in  a  soil  by  an 
inert  liquid  such  as  carbon  tetrachloride  or  benzene.  Both 
Leonards  and  Girault  (1961)  and  Mesri  (1973)  have  conducted 
such  experiments  and  have  found  that  the  nature  of  the  pore 
fluid  has  some  influence  on  the  rate  of  secondary 
consolidation.  However,  both  authors  agree  that  the  presence 
of  adsorbed  water  is  not  necessary  in  order  for  secondary 
consolidation  to  occur. 
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Barden  (1965)  summarized  the  work  of  Terzaghi  and 
Taylor  and  contended  that  they  are  in  agreement  on  the 
following  points: 

1.  Primary  and  secondary  consolidation  are  part  of  a  single 
continuous  process. 

2.  The  seat  of  secondary  or  creep  effects  is  gradual 
readjustment  or  remoulding  of  the  soil  structure 
initiated  during  primary  consolidation. 

3.  The  rate  at  which  the  secondary  compression  proceeds  is 
strongly  influenced  by  the  viscous  effects  of  the 
adsorbed  water  layer. 

The  hypothesized  influence  of  the  viscosity  of*  the 
adsorbed  water  bonds  lends  some  fundamental  basis  to  the 
development  of  rheological  models  to  predict  secondary 
consolidation.  A  rheological  model  is  some  combination  of 
springs,  dashpots  and  sliders  that  has  a  mathematical 
analogy  that  fits  observed  soil  behaviour.  The  dashpots 
represent  the  adsorbed  water  bonds.  Because  the  viscosity 
hypothesis  is  tenuous  rheological  models  must  remain 
primarily  phenomenological. 

Numerous  rheological  models  have  been  proposed  for  the 
stress-strain-time  behaviour  of  soils.  Barden  (1965),  de 
Jong  (1968)  and  Mitchell  (1976)  summarize  the  majority  of 
these  models.  The  more  important  rheological  models  are 
those  by  Gibson  and  Lo  (1961),  Christensen  and  Wu  (1964)  and 
Komamura  and  Huang  (1974).  Gibson  and  Lo  ( op.cit .  )  developed 
a  rheological  model  that  accounts  for  both  primary  and 
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secondary  consolidation.  Taylor  and  Merchant  (1940)  were  the 
first  to  include  secondary  compression  in  a  consolidation 
theory  but  did  not  relate  their  theory  to  a  rheological 
model.  Christie  (1964)  demonstrated  that  the  rheological 
model  developed  by  Gibson  and  Lo  (1961)  is  equivalent  to  the 
Taylor  and  Merchant  (1940)  theory. 

Mitchell  (1976)  notes  that  most  mathematical 
expressions  from  rheological  models  are  complicated  and  none 
has  yet  been  proposed  which  is  as  simple  as  the  three 
parameter  relationship  presented  later  as  Equation  3.2. 
Simple  rheological  models  which  would  have  potential  use  are 
not  used  because  they  incorrectly  predict  infinite 
settlement.  Due  to  their  complexity  rheological  models  are 
not  used  extensively  in  creep  modelling  of  soils  and  will 
not  be  discussed  further. 

3.2.3  Rate  Process  Theory 

Rate  process  theory  is  based  upon  the  proposition  that 
atoms,  molecules  or  particles  (termed  flow  units) 
participating  in  a  time  dependent  deformation  process  are 
constrained  from  movement  relative  to  each  other  by  virtue 
of  energy  barriers  which  separate  adjacent  energy  positions. 
In  order  for  a  flow  unit,  a  soil  particle  in  this  case,  to 
move  to  a  new  equilibrium  position  it  must  acquire 
sufficient  energy  to  surmount  the  imposed  energy  barrier. 
This  energy  is  called  the  activation  energy  and  it  has  some 
frequency  of  occurrence.  A  shear  stress  distorts  the  energy 
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barriers  allowing  flow  units  to  move  preferentially  in  the 
direction  of  distortion. 

Rate  process  theory  was  first  proposed  by  Eyring  (1936) 
and  became  popular  in  the  soil  mechanics  literature  in  the 
1960  s.  Mitchell  (1976)  provides  a  compact  overview  of  its 
significance  to  different  aspects  of  soil  mechanics. 
Mitchell  Gt.sl  (1968)  proposed  the  following  equation 
derived  from  rate  process  theory  to  determine  the  creep 
strain  rate: 
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(2  X) ( KT/h ) exp ( - aF/RT )  sinh(f^/2KT)  3.1 

=  a  function  of  the  number  of  flow  units  in 
the  direction  of  deformation  and  the 
average  component  of  displacement  due  to  a 
single  unit  mounting  the  energy  barrier 
=  Boltzmann’s  constant,  1.38  x  10_1‘erg/°  K 
=  Planck's  constant,  6.624  x  10"27erg/sec 
=  °  K 

=  free  energy  of  activation  calories/mole 
=  universal  gas  constant,  1.98  cal/°K-mole 
=  average  shear  force  on  each  flow  unit, 
dynes 

=  separation  distance  between  successive 


equilibrium  positions 
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Several  authors,  including  Murayama  and  Shibata  (1964), 
Mitchell  (1964),  Christensen  and  Wu  (1964),  Wu  et.al  (1966) 
and  Krizek  et .  al (1977),  have  shown  how  shear  strength, 
consolidation  and  creep  of  soils  can  be  related  to  rate 
process  theory.  Due  to  its  complexity  rate  process  theory 
has  not  been  widely  used  in  practice  and  remains  primarily  a 
research  tool. 

3.2.4  A  Model  Based  on  a  Cavity  Channel  Network 

De  Jong  (1968)  introduced  a  concept  of  primary  and 
secondary  consolidation  that  models  the  soil  pores  as 
cavities  of  differing  compressibilities  linked  together  by 
channels  with  differing  permeabilities.  During  primary 
consolidation  pore  water  drains  from  the  accessible  cavities 
and  the  ambient  excess  pore  water  pressure  eventually  goes 
to  zero.  However,  in  a  certain  number  of  cavities  linked  to 
the  network  by  channels  of  low  permeability,  the  excess  pore 
pressure  remains  greater  than  zero.  The  slow  drainage  of 
water  from  these  cavities  gives  rise  to  secondary 
consolidation . 

De  Jong  (op.cit .)  develops  a  constitutive  relationship 
from  this  model  in  a  rigourous  stochastic  manner  and  proves 
that  the  cavity  channel  model  yields  the  same  equations  as  a 
rheological  model  also  developed  in  the  same  paper.  Due  to 
its  complexity  the  model  has  only  been  employed  by  a  few 
other  researchers  (see  Holzer  et.al,  1973). 
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3.3  The  Phenomenological  Approach 

The  phenomenological  approach  is  the  endeavor  to 
produce  empirical  relationships  that  predict  creep  behaviour 
based  solely  on  observed  creep  behaviour  in  the  laboratory. 
Singh  and  Mitchell  (1968)  contend  that  an  empirical  creep 
relationship  must  satisfy  the  following  general 
requirements: 

1.  It  must  be  applicable  to  a  reasonable  range  of  creep 
stresses . 

2.  It  must  describe  the  behaviour  of  a  range  of  soil  types. 

3.  It  must  account  for  both  linear  and  curved  relationships 
between  strain  and  the  logarithm  of  time. 

4.  It  must  contain  parameters  that  are  easily  determined. 

The  most  widely  accepted  empirical  relationship  for 
soil  is  that  by  Singh  and  Mitchell  (1968)  which  has  the 
following  form: 


where  i 
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a 


D 

t 

1 
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m 

I  =  A<FD(  1/t)  3.2 

=  strain  rate 

=  a  constant,  the  strain  rate  at  unit  time 
and  zero  deviator  stress 

=  slope  of  the  natural  logarithm  of  strain 
rate  versus  deviator  stress  curves 
=  deviator  stress 
=  time 
=  unit  time 

=  slope  of  the  natural  logarithm  of  the 
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strain  rate  versus  natural  logarithm  of 
time  curves 

The  stress  level  term  in  Equation  3.2  can  be  replaced  by: 

a  D  =  a  D  3.3 

where  a  =aDrnay. 

D  =  D/Dm<u  =  stress  level 

Dmq*  =  deviator  stress  at  failure 

In  order  to  determine  the  constants  in  Equation  3.2  it  is 
necessary  to  perform  at  least  two  creep  tests  at  the  same 
consolidation  stress  but  at  different  stress  levels.  The 
data  is  then  plotted  as  the  natural  logarithm  of  strain  rate 
versus  the  natural  logarithm  of  time  for  each  stress  level 
to  determine  m  as  shown  in  Figure  3.2.  The  same  data  is  then 
plotted  in  a  natural  logarithm  of  strain  rate  versus 
arithmetic  deviator  stress  plot  to  determine  a  and  A  as 
shown  in  Figure  3.3. 

Mitchell  (1976)  presented  a  detailed  derivation  of 
Equation  3.2.  Numerous  investigations  of  the  creep  behaviour 
of  natural  materials  since  1968  have  demonstrated  that 
Equation  3.2  applies  to  a  wide  range  of  soils  both  in 
undrained  and  drained  conditions  and,  also,  to  some  types  of 
rock  (see  da  Fontoura,  1979).  Singh  and  Mitchell  (1968)  have 
shown  that  Equation  3.2  can  be  derived,  in  part,  from  rate 
process  theory.  A  most  important  feature  of  Equation  3.2  is 
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2  Strain  rate  versus  time  during  undrained  creep  of 
Osaka  clay  (after  Murayama  and  Shibata,  1958) 


Devietor  Stress  (kg/cm2) 


3  Variation  of  strain  rate  with  deviator  stress  for 
creep  of  San  Francisco  Bay  Mud  (after  Mitchell, 
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that  it  predicts  a  continual  decrease  of  strain  rate  with 
time.  In  other  words,  steady  state  creep  does  not  exist  for 
the  materials  described  by  this  equation.  Equation  3.2  only 
applies  for  stresses  below  the  long  term  strength. 

Equation  3.2  was  derived  on  the  basis  of  laboratory 
data  taken  primarily  from  the  triaxial  test.  In  an  attempt 
to  develop  creep  relationships  that  adequately  describe  soil 
in  a  general  state  of  stress,  Kavazanjian  and  Mitchell 
(1977)  found  it  convenient  to  decompose  the  general  creep 
strain  and  stress  tensor  into  its  volumetric  and  deviatoric 
components  as  follows: 


6  =  bv/+  60 


3.4 


where  0 
Ov 


=  general  stress  tensor 
=  volumetric  stress  tensor 
=  deviatoric  stress  tensor 


and  € 

€o 


=  general  strain  tensor 
=  volumetric  strain  tensor 
=  deviatoric  strain  tensor 


3.5 


The  definition  of  the  stress  and  strain  tensors  are  as  given 
in  Fung  (1965).  They  state  that  each  strain  tensor  component 
is  related  to  each  stress  tensor  component  by  a  pseudo 
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elastic  formulation  of  the  type: 

c  =  1  /Ra  6y'  +  1/6  Adp  3.6 

=  volumetric  tensor  operator 
=  deviatoric  tensor  operator 

The  characteristics  of  the  tensor  operators  have  been 
discussed  in  detail  by  Kavazanjian  and  Mitchell  (op.cit.) 

It  is  convenient  to  review  the  creep  deformation  of 
soils  in  terms  of  the  Kavazanjian  and  Mitchell  (op.cit .  ) 
formulation  and  in  terms  of  the  drainage  conditions.  The 
volumetric  and  deviatoric  components  of  the  strain  tensor 
for  drained  creep  will  be  discussed  first.  This  is  followed 
by  a  discussion  of  undrained  creep. 

3.3.1  Drained  Creep 

Drained  creep  is  time  dependent  deformation  under 
sustained  constant  effective  stress  where  the  boundary 
conditions  during  testing  allow  volume  change.  Drained  creep 
deformation  has  a  volumetric  and  deviatoric  component  in  all 
stress  states  in  the  laboratory  with  the  exception  of 
isotropic  consolidation  in  a  triaxial  test  where  the 
deviatoric  component  is  zero. 

3.3.  1.1  Volumetric  Component 

The  volumetric  component  of  drained  creep  is 

traditionally  referred  to  in  the  literature  as  secondary 


where  R 
6 
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compression  or  secondary  consolidation.  It  is  referred  to  as 
secondary'  because  it  has  been  considered  as  deformation 
occurring  after  primary  consolidation.  This  notion  is,  of 
course,  incorrect  (see  Gibson  and  Lo,  1961)  as  creep 

deformations  must  also  occur  during  the  dissipation  of 
excess  pore  pressures.  Both  terms  will  be  used 

interchangeably  in  this  thesis.  The  writer  prefers  the  term 

drained  creep. 

A  substantial  amount  of  research  has  been  conducted  in 
the  laboratory  on  the  factors  influencing  secondary 
compression.  Most  researchers  have  chosen  to  use  the 

oedometer  for  the  investigation  of  secondary  compression 
because  of  its  simplicity  and  because  the  measurement  of 
volume  change  in  other  testing  devices  over  long  periods  of 
time  presents  difficulties.  Since  oedometer  testing  is  so 
common  it  is  of  value  to  review  briefly  the  boundary 
conditions  of  this  testing  device. 


£  2  —  £  3  =  0 


3.7 


It  follows  that: 


where  £ , _ 


e  e  1  Imarx. 

principal  strains 
volumetric  strain 
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Irnav  =  maximum  engineering  shear  strain 


The  effective  stresses  are  related  by: 
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where  6 )  _  3 
K0 


^2  =  o;  =  K o  ©!  3.9 

=principal  stresses 

=  coefficient  of  earth  pressure  at  rest 


Other  testing  equipment  such  as  the  triaxial  test,  the 
plane  strain  test  and  the  simple  shear  plane  strain  test 
have  also  been  used  to  investigate  the  volumetric  component 
of  the  drained  creep  tensor.  It  is  customary  in  all  tests  to 
determine  the  amount  of  secondary  compression  by  plotting 
the  void  ratio,  e,  against  the  logarithm  of  time.  The  slope 
of  the  curve,  shown  in  Figure  5.1,  after  the  dissipation  of 
excess  pore  pressures  is  referred  to  as  the  secondary 
compression  index  (after  Mesri  and  Godlewski,  1977)  and  is 
defined  by: 


where  a 
e 
t 


a  =  e/  log  t 

=  secondary  compression  index 
=  void  ratio 
=  time 


3.10 


The  coefficient  of  secondary  consolidation,  Ca  is  given  by: 

Ca  =  a/1 +e  o  3.11 

The  secondary  compression  index  may  be  a  function  of 
time,  stress  history,  consolidation  stress,  load  increment 
ratio,  shear  stress  level,  stress  state,  stress  path, 
duration  of  previous  stress  and  temperature.  The  influence 


“  '  :  '  ' 
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of  all  these  factors  will  now  be  discussed.  Except  where 
noted,  the  researchers  mentioned  in  the  following  discussion 
have  used  the  oedometer. 

Time 

It  is  widely  assumed,  as  first  proposed  by  Buisman 
(1936),  that  a  is  independent  of  time.  That  is,  the  slope  of 
the  void  ratio  versus  logarithm  of  time  curve  is  constant. 
Bjerrum  (1967),  in  his  Rankine  lecture,  assumed  that  Ca  was 
independent  of  time,  as  did  Walker  (1969a)  and  Kavazanjian 
and  Mitchell  (1980).  However,  Lo  (1961)  examined  the 
behaviour  of  both  remoulded  and  undisturbed  clays  and  found 
that  Ca  decreased,  remained  constant  or  increased  with  time 
depending  on  the  soil  type. 

Lo  ( op.cit .)  also  found  that  secondary  compression  went 
to  zero  with  time,  usually  within  three  weeks.  Conversely 
other  researchers,  notably  Bishop  and  Lovenbury  (1969), 
found  that  secondary  compression  from  triaxial  and  oedometer 
tests  on  undisturbed  clay  continued  for  over  three  years. 

It  is  of  value  to  further  elaborate  on  the  dependence 
of  Ca  on  time.  Equation  3.11  can  also  be  written  as: 

?  =  (Ca) ( 1/t)  3.12 

where  e  =  volumetric  strain  rate 

If  the  natural  logarithm  of  volumetric  strain  rate  is 
plotted  versus  the  natural  logarithm  of  time,  as  shown  in 


.  1 
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Figure  5.5  the  curve  can  be  expressed  by  a  relationship  of 
the  type: 

m 

*  =  A( 1/t )  3.13 
The  effect  of  consolidation  shown  in  Figure  5.5  is  neglected 
in  the  derivation  of  Equation  3.13.  If  Equations  3.12  and 
3.13  are  equated  it  is  apparent  that  Ca  is  a  constant,  i.e. 
independent  of  time,  only  if  the  parameter,  m,  is  equal  to 
one.  This  has  been  previously  recognized  by  Singh  and 
Mitchell  (1969),  Barden  and  Poskitt  (1969)  and  Walker 
(1969c)  among  others. 

Although  Ca  may  not  be  independent  of  time  most 
researchers  have  found  the  parameter,  m,  in  Equation  3.2  to 
be  independent  of  time.  As  an  example,  Tavenas  et.al  (1978) 
have  plotted  their  volumetric  creep  data  from  triaxial  tests 
on  undisturbed  St.  Alban  clay  in  the  same  form  as  Equation 
3.2  and  have  found  the  parameter,  m,  to  be  independent  of 
time.  This  suggests  that  m  is  a  more  fundamental  parameter 
than  Ca. 

Consolidation  Stress 

There  is  some  controversy  in  the  literature  concerning 
the  dependence  of  Ca  on  consolidation  stress.  However,  most 
researchers,  including  Mesri  and  Godlewski  (1977)  and 
Murayama  and  Shibata  (1964),  agree  that  Ca  increases  with 
increasing  stress  below  the  preconsolidation  pressure,  pc, 
and  reaches  a  maximum  in  the  vicinity  of  pc.  It  is  above  the 


. 
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maximum  that  the  controversy  appears.  Some  researchers, 
including  Murayama  and  Shibata  (1964),  Schiffman  et.al 
(1964)  for  one  of  the  soils  they  tested,  Ladd  and  Preston 
(1965),  Walker  (1969a),  Lo  (1961)  andShen  et.al  (1973), 
contend  that  Ca  is  a  constant  above  pc.  Others,  including 
Wahls  (1962),  Schiffman  et.al  (1964)  and  Mesri  and  Godlewski 
(1977)  maintain  that  Ca  decreases  with  increasing  stress 
above  pc.  Figures  3.4  and  3.5  illustrate  examples  of  the 
dependence  of  Ca  on  consolidation  stress.  It  is  of  interest 
to  note  the  Mesri  and  Godlewski  (  1977)  only  tested’ soils 
that  displayed  dependence  of  the  compression  index,  Cc ,  on 
consolidation  stress  above  pc.  o 

Load  I ncrement  Ratio 

The  load  increment  ratio  is  defined  here  as  the  net 
added  load  divided  by  the  total  previous  load,Ap/p.  Several 
researchers,  including  Mesri  and  Godlewski  (1977),  Leonards 
and  Girault  (1961)  and  Gibson  and  Lo  (1961)  have  documented 
that  the  shape  of  the  consolidation  curve  is  dependent  on 
the  load  increment  ratio.  This  is  illustrated  in  Figure  3.6, 
taken  from  Bjerrum  (1967).  However  most  researchers, 
including  Wu  et.al  (1978),  Bjerrum  (1967),  Barden  (1969), 

Wahls  (1962)  and  Mesri  and  Godlewski  (1977),  have  concluded 
that  Ca  is  independent  of  the  load  increment  ratio. 


Shear  Stress  Level 
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Figure  3.4  Relation  between 
secondary  compression  in 
Shibata,  1964) 


applied  stress  and  rate 
oedometer (after  Murayama 


of 

and 


Figure  3.5  Relationship  between  Ca  and  Cc  and  consolidation 
pressure  for  (a)Mexico  City  clay  (b)Leda  clay  (c )Undi stured 
New  Haven  organic  silt (after  Mesri  and  Godlewski,  1977) 
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Figure  3.6  Effect 
compression  index 


TIME  IN  YEARS 


of  load  increment  ratio  on  secondary 
(after  Bjerrum,  1967) 
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The  effect  of  shear  stress  level  on  the  magnitude  of  Ca 
has  been  investigated  only  sparingly.  This  is  a  consequence 
of  most  research  being  undertaken  with  the  oedometer  where 
the  shear  stress  level  is  a  constant.  Murayama  and  Shibata 
(1961)  found  a  shear  stress  level  dependence  of  Ca  for 
undisturbed  Osaka  clay  using  the  triaxial  test.  Walker 
(1969a)  using  the  simple  shear  plane  strain  test  and  the 
triaxial  test  found  the  volumetric  creep  rate  to  be 
independent  of  shear  stress  level.  His  data  is  illustrated 
in  Figure  3.7.  Kavazanjian  and  Mitchell  (1980)  cited  Walker 
(1969a)  to  justify  their  assumption  that  Ca  is  independent 
of  stress  level.  Yudhbir  and  Mathur  (1977)  found  that  the 
volumetric  creep  rate  was  independent  of  shear  stress  level 
up  to  a  yield  value  where  it  decreased  to  zero. 

Analysing  their  creep  data  according  to  Equation  3.2 
Tavenas  et.al  (1978)  found  a  dependence  of  volumetric  strain 
rate  in  the  triaxial  test  on  shear  stress  level  and 
expressed  it  in  the  form: 

m 

l  =  Bf  (6’  )  (t  Jt )  3.14 
where  B  =  constant 

f ( 6 ’ )  =  stress  function 

They  found,  using  different  stress  paths,  that  the  stress 
function  could  be  represented  by  the  equations  of  lines  of 
equal  volumetric  strain  rate  in  stress  space.  They  did  not 
produce  these  equations  for  the  stress  function  but  did 
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suggest  they  would  take  the  same  form  as  the  equation  of  the 
limit  state  surface  (yield  surface)  for  the  clay.  Figure  3.8 
illustrates  the  shape  of  the  stress  function. 

Stress  State 

The  effect  of  stress  state  on  the  magnitude  of  Ca  has 
been  investigated  only  sparingly.  Again,  this  is  a 
consequence  of  most  research  on  Ca  being  conducted  using 
only  the  oedometer.  Walker  (1969a)  found  a  minimal  effect  of 
stress  state  on  the  value  of  Ca  between  the  triaxial  and 
simple  shear  devices  for  tests  on  remoulded  kaolinite. 
Schiffman  et.al  (1964)  found  little  difference  between  Ca 
determined  from  isotropic  triaxial  consolidation  tests  and 
oedometer  tests.  Bishop  and  Lovenbury  (1969),  from  oedometer 
and  triaxial  tests,  and  Wu  et.al  (1978),  from  triaxial, 
plane  stain  and  simple  shear  plane  strain  tests  found  little 
dependence  of  volumetric  creep  strain  or  strain  rate  on 
stress  state.  There  appears  to  be  a  consensus  that  stress 
state  has  only  a  minimal  effect  on  volumetric  creep. Of 
course, this  conclusion  is  only  valid  for  soils  which  have 
isotropic  creep  properties. 


Stress  Path 


Tavenas  et . 
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Figure  3.7  Effect  of  stress  level  on  volumetric  creep  rate, 
left  -  triaxial  tests,  right  -  simple  shear  tests  (after 
Walker,  1969b) 


Figure  3.8  Lines  of  equal  volumetric  strain  rate  at  t  -  100 
min  in  the  stress  space  for  drained  tests  (after  Tavenas 
et.al  ,  1978) 
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Yudhbir  and  Mathur  (1977),  found  that  the  volumetric  creep 
rate  was  dependent  on  stress  path. 

Duration  of  Previous  Stress 

The  duration  of  the  previous  stress  may  have  some 
effect  on  the  subsequent  volumetric  creep  rate  according  to 
Mesri  (1973).  Bjerrum  (1967)  demonstrated  that  delayed 
compression,  i.e.  secondary  compression,  induces  a  pseudo 
preconsolidation  stress  on  the  sample.  It  has  already  been 
shown  in  this  section  that  preconsolidation  has  an  effect  on 
the  volumetric  creep.  If,  however,  the  subsequent  load  is 
higher  than  the  pseudo  preconsolidation  stress  then  the 
duration  of  the  previous  load  has  a  minimal  effect.  Wu  et.al 
(1978)  and  Barden  (1969)  found  no  effect  on  the  subsequent 
creep  of  the  sample  from  the  duration  of  the  previous 
consolidation  load. 

Temperature 

An  increase  in  temperature  will  cause  an  increase  in 
the  creep  rate  for  all  creep  conditions.  Rate  process  theory 
demonstrates  this  admirably.  However,  the  temperature  must 
be  changed  substantially  according  to  Bishop  and  Lovenbury 
(1969)and  Mesri  (1973)  before  the  effect  is  measurable. 

3.3. 1.2  Deviatoric  Component 

The  deviatoric  component  of  the  drained  creep  strain 
tensor  has  not  received  much  attention  in  the  literature. 
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Again,  this  is  a  consequence  of  the  majority  of  researchers 
using  the  oedometer  for  the  investigation  of  drained  creep 
behaviour.  It  is  revealing  to  note  that  Equation  3.8 
demonstrates  that  the  maximum  engineering  shear  strain  in 
the  oedometer  is  equal  to  the  volumetric  strain.  It  then 
follows  that  the  maximum  shear  strain  rate  must  be  dependent 
on  the  same  factors  as  the  volumetric  strain  rate  for  the 
oedometer  stress  state.  Furthermore , it  should  be  a  function 
of  shear  stress  level  but  this  cannot  be  investigated  with 
the  oedometer. 

Kavazanjian  and  Mitchell  (1980)  state  that  deviatoric 
component  of  the  drained  creep  strain  tensor  can  be 
described  by  a  relationship  of  the  type  given  by  Equation 
3.2.  The  results  of  Tavenas  et.al  (1978)  are  in  agreement 
with  this.  Their  data  fits  a  relationship  of  the  form: 


where  i 


and 


Y  =  A[  g  ( <5  ’  )](t,/t) 
=  shear  strain  rate 
=  constant 


m 


g(d')  stress  function 


where 


£  =  -  1  /3£* 

=axial  strain  rate 


=volumetric  strain  rate 


3.15 


3.16 
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Figure  3.9  shows  the  form  of  the  shear  stress  function.  It 
is  apparent  that  the  shear  strain  rate  is  a  function  of 
shear  stress  level  but, for  their  data,  not  in  the  simplified 
form  given  by  Equation  3.2.  It  is  of  interest  to  note  that 
the  value  of  m  in  Equation  3.14  is  the  same  as  that  in 
Equation  3.15. 

Walker  (1969a)  also  investigates  the  influence  of  shear 
stress  level  on  shear  creep  rate,  a  parameter  which  is 
analogous  to  Ca,  for  triaxial  and  simple  shear  tests  on 
remoulded  kaolin.  His  results  are  illustrated  in  Figure 
3.10.  A  linear  relationship  between  shear  creep  rate  and 
shear  stress  level  is  indicated  in  Figure  3.10  which  is 
independent  of  consolidation  stress.  Yudhbir  and  Mathur 
(1977)  also  have  shown  that  the  shear  creep  rate  is 
dependent  on  stress  level  and  also  on  stress  path. 

In  summary , deviator ic  creep  strains  under  drained 
conditions  are  a  function  of  stress  level, stress  path  and 
time. A  relationship  of  the  Singh-Mi tchell  type  appears  to 
fit  the  laboratory  data  of  most  researchers  reasonably  well. 

3.3.2  Undrained  Creep 

Undrained  creep  results  from  the  application  of  a 
sustained  total  stress  under  conditions  of  no  volume  change. 
It  then  follows  that  the  creep  strain  tensor  contains  a 
deviatoric  component  only.  During  undrained  creep  the  pore 
pressures  change  with  time  and,  hence,  the  effective  stress 
during  undrained  creep  is  not  constant. 
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Figure  3.9  Lines  of  equal  shear  strain  rates  at  t  =  100  min 
in  the  stress  space  for  drained  and  undrained  tests  (after 
Tavenas  et.al  ,  1978) 


Figure  3.10  Effect  of  stress  level  on  shear  creep  rate,  left 
triaxial  tests,  right  -  simple  shear  tests  (after  Walker, 
1969a) 
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The  most  common  laboratory  apparatus  for  performing 
undrained  creep  experiments  is  the  triaxial  test.  A  few 
researchers  have  also  used  the  plane  strain  device  and  the 
simple  shear  plane  strain  device.  The  prevention  of  leaks 
during  undrained  testing  in  all  three  devices  is  difficult. 
Poulos  (1964)  presented  a  comprehensive  summary  of  sources 
of  leaks  in  the  triaxial  test  which  is  of  great  value  for 
designing  undrained  creep  test  programs. 

For  undrained  creep  testing  it  is  customary  to  plot  the 
development  of  axial  strain  and  pore  pressure  with  time  for 
different  stress  levels.  Typical  examples  of  such  pore 
pressure  plots  are  shown  in  Figures  3.12  and  3.13.  Singh  and 
Mitchell  (1968),  Edgers  et.al  (1973)  and  Holzer  et.al  (1973) 
all  found  the  axial  strain  rate  to  be  adequately  represented 
by  Equation  3.2.  Tavenas  et.al  (1978)  found  that  axial 
strain  rate  can  be  described  by  a  relationship  similar  to 
Equation  3.2  but  with  a  stress  function  that  has  a  similar 
form  to  that  shown  in  Figure  3.9. 

It  is  of  interest  to  note  that  Tavenas  et.al  (1978) 
found  the  parameter,  m,  in  Equation  3.2  to  be  independent  of 
drainage  conditions  and  stress  path.  This  is  in  agreement 
with  Singh  and  Mitchell's  (1968)  contention  that  the 
parameter,  m,  is  a  material  constant.  Edgers  et.al  (1973)  in 
an  extensive  series  of  simple  shear  tests  with  no  volume 
change  and  undrained  triaxial  tests  on  undisturbed,  highly 
plastic,  very  soft,  backswamp  clay  found  that  m  was 
dependent  on  stress  state  and,  also,  a  function  of  whether 


;  d"  1 
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the  sample  was  consolidated  isotropically  or 
an iso tropically.  They  also  found  the  stress  level  term,  a, 
in  Equation  3.2  to  be  a  function  of  stress  state. 

Krizek  et  .al  (1977)  performed  a  series  of  undrained 
creep  tests  on  remoulded  kaolin  which  was  prepared  with  a 
number  of  different  microfabrics.  The  microfabrics  varied 
from  highly  random  to  highly  oriented  and  were  identified 
with  the  scanning  electron  microscope.  They  found  that  m  was 
a  function  of  void  ratio  and  anisotropy  as  shown  on  Figure 
3.11.  The  vertical  specimens  exhibited  larger  creep 
deformat  ions ,/ .e.  a  lower  value  of  m,  than  the  samples  cut 
in  the  horizontal  direction.  The  stress  level  parameter,  a, 
was  also  found  to  be  a  function  of  void  ratio  and 
anisotropy . 

Walker  (1969b)  conducted  a  series  of  undrained  triaxial 
tests  on  Leda  clay  to  determine  the  development  of 
deviatoric  strain  and  pore  pressure  with  time.  He  chose  to 
plot  these  two  variables  against  the  logarithm  of  time  and 
found  a  bilinear  relationship  for  each.  The  initial  linear 
relationship  was  attributed  to  redistribution  of  pore 
pressure.  He  proposed  a  functional  relationship  between 
axial  strain  and  pore  pressure  which  was  based  on  their 
similar  relationship  to  time .Unfortunately  an  exact 
formulation  of  this  function  was  not  given. 


Holzer  et.al  (1973)  ran  a  series  of  undrained  triaxial 
tests  on  undisturbed  San  Francisco  Bay  Mud  to  investigate 
the  nature  of  undrained  pore  pressure  response  during  creep. 
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Open  symbol* :  H  specimens 
Solid  symbols :  V  specimens 
□  Dispersed  anisotropic 
A  Flocculated  anisotropic 
O  Dispersed  and  flocculated  isotropic 


Figure  3.11  Creep  parameter, 
Krizek  et.al  ,  1977) 


m,  versus  void  ratio  (after 
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They  found  that  the  pore  pressure  response  could  be 
decomposed  into  two  parts.  The  first  was  essentially  time 
independent  and  due  to  the  application  of  an  increased 
deviator  stress.  The  second  was  time  dependent  and  due  to 
the  arresting  of  secondary  consolidation.  They  were  able  to 
decompose  the  pore  pressure  response  in  this  way  by 
determining  the  pore  pressure  response  of  an  isotropically 
consolidated  sample  at  zero  deviatoric  stress.  Their  data  is 
presented  in  Figures  3.12  and  3.13.  It  can  be  seen  from 
Figure  3.12  that  increased  time  for  secondary  consolidation 
resulted  in  a  lowering  of  the  time  dependent  pore  pressure 
response  during  shear. 


3.4  Creep  Rupture 

A  soil  may  fail  under  the  sustained  action  of  a  creep 
stress  at  undrained  strengths  substantially  lower  than  the 
strength  measured  in  ’normal’  undrained  tests.  A  failure  of 
this  type  is  referred  to  as  creep  rupture.  Mitchell  (1976) 
attributes  the  discrepancy  between  the  normal  undrained 
strength  and  the  creep  rupture  strength  to  either 
deterioration  of  cemented  bonds  by  creep  strains  or  pore 
pressure  generation  under  undrained  conditions  which 
decrease  the  effective  stresses  to  failure. 
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Time  (thousand  min) 
Consolidation  Pressure  =1.0  kg/cm2 


Figure  3.12  Excess  pore  pressure  development  following 
cessation  of  drainage  after  secondary  compression  for 
different  time  periods  (after  Holzer  et .al ,  1973) 


Figure  3.13  Excess  pore  pressure  development  during 
undrained  creep  of  San  Francisco  Bay  mud  after  consolidation 
at  1.0  kg/cm2  for  10  min  (after  Holzer  et.al  ,  1973) 
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3.4.1  Undrained  Conditions 

It  is  important  to  appreciate  that,  although  the 
undrained  creep  rupture  strength  is  substantially  less  than 
the  normal  undrained  strength,  failure  is  still  governed  by 
the  effective  strength  envelope.  Campanella  and  Vaid  (1974) 
demonstrate  this  fact  conclusively  using  a  comprehensive 
series  of  undrained  creep  test  results  on  undisturbed  Haney 
clay . 

Finn  and  Snead  (1973)  conducted  a  series  of  undrained 
creep  rupture  tests  on  undisturbed  Haney  Clay  to  determine 
the  cause  of  the  discrepancy  between  creep  rupture  strength 
and  the  ’normal’  undrained  strength.  They  found  that  the 
creep  rupture  strength  could  determined  from  a  relationship 
of  the  form: 


where  6C 

6uy 

K 

e 

n 


6c  -  6u  y  +  K  * 

=  creep  strength 
=  upper  yield  strength 
=  a  constant 

=  minimum  axial  strain  rate 
=  constant,  3  for  Haney  Clay 


3.17 


Figure  3.14  presents  Equation  3.17  in  graphical  form.  The 
significance  of  the  minimum  axial  strain  rate  is  shown  in 
Figure  3.15.  The  upper  yield  strength,  first  proposed  by 
Murayama  and  Shibata  (1958)  is  the  creep  stress  below  which 
failure  will  not  occur.  According  to  Finn  and  Snead  (1973) 
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Figure  3.14  Determination  of  upper  yield  stresses  (after 
Campanella  and  Vaid,  1974) 


Time  (min) 


Figure  3.15  Creep  rate  behaviour  of  K 0 -consol idated 
undisturbed  Haney  clay  under  axially  symmetric  loading 
(after  Campanella  and  Vaid,  1972) 
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it  is  equivalent  to  the  undrained  strength  determined  from  a 
constant  strain  rate  test  run  at  a  strain  rate  sufficiently 
slow  to  allow  complete  equilization  of  pore  pressure  within 
the  sample.  It  then  follows  that  the  discrepancy  between 
'normal’  undrained  strengths  and  creep  rupture  strength  is  a 
result  of  the  high  strain  rates,  1  %/min  to  2  %/min,  usually 
employed  for  'normal'  undrained  strength  tests. 

3.4.2  Creep  Rupture  Life 

It  is  of  interest  to  be  able  to  predict  the  time  to 
rupture  in  a  creep  test  when  the  creep  stress  is  greater 
than  the  upper  yield  strength.  The  time  period  from  the 
initiation  of  creep  to  creep  rupture  is  termed  the  creep 
rupture  life.  Saito  (1965)  presents  a  relationship  between 
creep  rupture  life  and  any  transient  strain  in  the 
'tertiary'  range  as: 


where  * 


tr 

t 

a 


?  =  a/(tr  -t)  3.18 

=  transient  strain  rate  in  the  'tertiary' 
range  at  optional  time 
=  creep  rupture  life 
=  optional  time 
=  constant 


Saito  ( op.cit .  )  demonstrates  by  graphical  means  how  Equation 
3.18  can  be  used  to  correctly  predict  creep  rupture  life. 
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Singh  and  Mitchell  (1969) 

have 

also 

presented 

an 

empirical  method  for  predicting 

the 

time 

remaining 

to 

failure  for  creep  rupture  tests.  They  plot  It,  which  they 
called  the  creep  coefficient,  versus  the  logarithm  of  time. 
This  method  applies  only  to  soils  with  the  parameter,  m, 
less  than  one.  They  found,  for  a  given  soil,  that  the  value 
of  the  creep  coefficient  reached  a  unique  value  just 
preceding  the  onset  of  failure.  They  calculated  the  time  to 
failure  from: 


where  tr 
t, 
a 

?t; 

?t, 


<tf  /t,  ) '  =  ft  /ft,  3.19 

=  time  to  failure 
=  reference  time 
=  slope  of  It  versus  log  of  t 
=  creep  coefficient  at  failure,  unique  and 
known 

=  creep  coefficient  at  reference  time 


Extending  the  work  of  Sai to ( 1 965 ) , Finn  and  Shead  (1973) 
found  that  creep  rupture  life  was  related  to  minimum  strain 
rate  during  a  creep  rupture  test  by  a  relationship  of  the 
form: 

log10t  =  0 . 75  1  -0 . 92  log  ,  0  *m'in  3.20 

where  t  =  creep  rupture  life 

?n;n  =  minimum  creep  strain  rate 
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Equation  3.20  is  similar  to  that  proposed  by  Saito  and 
Uezawa  (1961).  Finn  and  Shead  (op. ci t .) state  that  Equation 
3.20  is  valid  for  a  variety  of  soils. 

A  disadvantage  of  Equation  3.20  is  that  the  time  since 
initiation  of  creep  must  be  known  in  order  to  calculate  the 
remaining  creep  rupture  life.  Finn  and  Shead  (1973)  overcame 
this  disadvantage  by  developing  a  relationship  of  the  form: 

ttr=  C/t  3.21 

where  t±r  =  time  remaining  to  creep  rupture 
C  =  constant 

*  =  current  accelerating  creep  rate 


Campanella  and  Vaid  (1974)  have  furthered  this  approach  to 
include  the  effect  of  stress  state.  They  performed  a  series 
of  creep  tests  on  undisturbed  Haney  clay  with  the  triaxial 
device  using  both  anisotropic  and  isotropic  consolidation 
and  with  the  plane  strain  device.  They  found  that  the  total 
rupture  life  determined  from  the  isotropic  triaxial  tests 
was  four  times  longer  than  that  determined  from  K„  triaxial 
and  plane  strain  tests. This  is  illustrated  in  Figure  3.16 
which  also  shows  the  relationship  given  by  Equation  3.20. 
Campanella  and  Vaid  ( Op.cit .)  also  demonstrated  that  the 
axial  strain  at  the  minimum  strain  rate  is  independent  of 
the  creep  stress.  This  is  shown  in  Figure  3.17  which  also 
illustrates  the  much  higher  creep  strains  from  the  isotropic 
triaxial  tests. 
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Figure  3.16  Relationships  between  rupture  life 
creep  rate--normally  consolidated  undisturbed 
(after  Campanella  and  Vaid,  1974) 


?  =  <0 1  -  03)/o'ic 


Figure  3.17  Axial  strain  at  minimum  strain 
function  of  creep  stress  for  undisturbed  Haney 
Campanella  and  Vaid,  1974) 


and  minimum 
Haney  clay 
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Mitchell  (1976)  provided  a  simple  method  to  obtain  the 
constant,  C,  in  Equation  3.20  using  the  parameters  from 
Equation  3.2.  He  demonstrated  that  C  is  equal  to: 

C  =  (  1-m)ef  3.22 
where  =  strain  at  failure,  a  constant 

The  constant,  C,  can  then  be  used  in  the  following 
relationship  to  find  the  remaining  time  to  failure: 

In  t{  =  1  / (  1  -m )  |lh(C/A)-aD|  3.23 
Equations  3.20  and  3.23  can  be  used  to  predict  the  remaining 
time  to  failure  using  only  the  parameters  in  Equation  3.2. 

It  must  be  emphasized  that  the  preceding  only  applies 
to  stresses  above  the  upper  yield  strength.  Creep  rupture 
will  not  occur  below  the  upper  yield  strength  given  any 
amount  of  time. 

3.4.3  Drained  Conditions 

Strength  loss  during  creep  under  drained  conditions  has 
received  much  less  attention  in  the  literature  than  under 
undrained  conditions.  Campanella  and  Vaid  (1974)  explained 
this  by  suggesting  that  if  creep  rupture  is  a  possibility 
then  it  is  more  likely  to  happen  before  drainage  can 
strengthen  the  soil. 
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Bishop  and  Lovenbury  (1969)  performed  long  term  drained 
creep  tests,  exceeding  three  years,  at  stress  levels  near 
peak  strength  and  found  no  reduction  in  the  effective  stress 
parameters  as  a  result  of  creep  strains.  This  was 
demonstrated  for  both  London  clay,  a  strain  softening  soil, 
and  Pancone  clay,  a  normally  consolidated  soil.  Shibata  and 
Karube ( 1 969 )  and  Tavenas  et.al  (1978)  also  found  no  strength 
reduction  as  a  result  of  creep  strains  in  drained  tests. 
These  laboratory  results  cast  doubt  on  Nelson  and  Thompson's 
(1977)  hypothesis  that  the  application  of  creep  stresses  to 
strain  softening  soils  at  stress  levels  between  residual  and 
peak  strength  eventually  cause  a  strength  reduction  of  the 
soil  to  residual  as  a  result  of  deterioration  of  soil  bonds. 

The  laboratory  results  from  the  literature  demonstrate 
that  strengths  derived  from  long  term  stress  controlled 
drained  creep  tests  are  identical  to  those  obtained  from 
drained  strain  controlled  tests  on  the  same  material.  Thus, 
it  appears  that  the  concept  of  creep  rupture  does  not  apply 
to  drained  creep  conditions  for  soils. 
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4.  Laboratory  Testing  Equipment  and  Procedures 


4 . 1  Introduction 

In  order  to  assess  the  creep  behaviour  of  Tar  Island 
clay  a  series  of  drained  creep  tests  were  undertaken  using 
the  oedometer  and  the  triaxial  test.  This  chapter  is  devoted 
to  the  description  and  design  of  the  test  equipment  and  to 
an  outline  of  the  test  procedures  used  for  creep  testing. 
Recommendations  are  made  for  improvements  in  testing 
equipment . 

Available  equipment  was  used  for  the  oedometer  testing. 
The  triaxial  systems  were  designed  and  fabricated  especially 
for  this  testing  program.  The  oedometer  testing  program  was 
undertaken  first  so  that  preconsolidation  stresses  could  be 
calculated  and  used  as  a  guide  in  the  design  of  the  minimum 
pressure  requirements  for  the  triaxial  cells.  All  triaxial 
creep  tests  were  performed  at  consolidation  stresses  above 
the  apparent  preconsolidation  stress  for  the  clay.  The  creep 
strains  determined  from  the  oedometer  testing  were  also  used 
as  a  guide  to  the  required  resolution  of  the  monitoring 
systems  for  the  triaxial  tests. 

The  results  of  the  testing  program  described  in  this 
chapter  are  presented  in  Chapter  5. 
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4.2  Oedometer  Creep  Tests 

4.2.1  Testing  Equipment 

The  oedometer  tests  were  performed  in  a  temperature  and 
humidity  controlled  room  which  was  remote  from  any  source  of 
vibration.  Room  temperature  was  maintained  within  a  range  of 
20 . 5°C  ± 1 °C . 

Three  samples  were  tested  simultaneously  using  three 
floating  ring  Clockhouse  Engineering  Ltd.  (Model  J42) 
oedometers.  Stress  was  applied  to  the  top  of  the  samples  by 
placing  dead  weights  on  a  self  compensating  lever  arm  system 
which  had  an  11  to  1  load  enhancement  factor.  A  reduction  to 
the  vertical  stress  was  not  made  for  oedometer  ring  friction 
because  it  is  generally  considered  to  be  negligible.  Before 
testing,  the  compliance  of  each  oedometer  was  investigated 
to  ascertain  if  there  were  any  system  time  effects. 

For  the  compliance  testing  an  aluminum  disc  equal  in 
diameter  to  the  inside  of  the  oedometer  ring  was  used  in 
place  of  the  soil  sample.  The  disc  was  enclosed  with  a 
similar  filter  paper  and  porous  plate  arrangement  to  that 
used  during  actual  testing.  Weights  were  added  in  the  same 
increments  as  for  soil  testing.  Displacement  versus  time 
readings  were  recorded  for  each  load  increment  until  the 
difference  between  successive  one  minute  readings  was  zero. 
The  results  from  a  typical  compliance  test  are  shown  in 
Figure  4.1.  The  top  diagram  in  Figure  4.1  reveals  that  most 
of  the  compliance  occurs  within  30  seconds  of  adding  the 


. 
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Time,  min 


Figure  4.1  Typical  oedometer  compliance  test  results 
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load.  The  bottom  diagram  indicates  some  hysteresis  in  the 
loading- unloading  compliance  curves. 

The  cause  of  the  hysteresis  may  be  irrecoverable 
deformation  of  the  porous  plates.  However,  it  was  observed 
that  a  slight  change  in  the  alignment  of  the  lever  arm 
system  induced  a  change  in  the  linear  variable  displacement 
transducer  (LVDT)  output.  It  is  considered  that  most  of  the 
hysteresis  is  caused  by  this  phenomenom  which  occurs  only 
upon  removal  of  the  weights.  The  immediate  compliance 
displacement  from  the  loading  curve  was  subtracted  from  the 
apparent  immediate  deformation  of  the  soil  for  each  loading 
increment  as  part  of  the  data  reduction.  It  was  found  that 
the  compliance  accounted  for  most  of  the  apparent  immediate 
deformation  of  soil  which  is  in  agreement  with  the 
theoretical  postulate  that  immediate  deformation  of  a 
saturated  soil  under  an  oedometer  stress  state  is  zero. 

Hewlett  Packard  Model  No.  7DCT-1000  LVDT's  were  used  to 
record  displacement.  The  LVDT's  were  attached  to  the 
oedometers  with  a  bracket.  The  threaded  bottom  portion  of 
the  LVDT  extension  rod  was  screwed  into  the  top  of  the  lever 
arm  immediately  above  the  sample.  This  mounting 
configuration  reduced  the  disturbance  from  accidental 
knocking  or  vibration  on  the  displacement  readings. 

The  transducers  were  calibrated  with  a  micrometer 
before  and  after  testing.  The  transducers  had  a  linear 
travel  of  2.54  cm  with  a  concomitant  output  of  7.500  volts. 
This  corresponds  to  a  theoretical  sensitivity  for  the  LVDT 
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of  ±. 0003cm.  However ,  the  observed  sensitivity  was  in  the 
order  of  ±.001cm.  The  excitation  voltage  for  each  transducer 
was  continually  monitored  for  the  duration  of  each  test. 
This  was  necessary  since  a  change,  for  any  reason,  in  the 
excitation  voltage  alters  the  calibration  factor  for  the 
transducer . 

The  output  voltages  of  the  transducers  were  monitored 
with  time  on  a  Fluke  33003  2200A  data  acquisition  system 
which  produced  a  permanent  record  on  paper  tape  and/or  on  a 
magnetic  tape  cassette  which  was  housed  in  a  Techtran  8410. 

4.2.2  Sample  Preparation 

Samples  for  oedometer  testing  were  obtained  from  the 
shelby  tube  samples  which  were  acquired  during  the  drilling 
program  described  in  Chapter  2.  The  shelby  tube  samples, 
wrapped  in  cheese  cloth  and  waxed,  were  stored  in  a  100 
percent  humidity  room  to  prevent  moisture  loss  before 
testing . 

Three  samples  which  represented  the  range  in  soil 
properties  available  from  the  shelby  tube  samples  were 
selected  for  oedometer  testing.  The  description  of  the 
samples  is  given  in  Chapter  5.  The  samples  were  cut  to  the 
inner  dimensions  of  the  oedometer  ring  in  the  conventional 
way  of  carefully  pushing  the  cutting  edge  of  the  oedometer 
ring  into  the  sample  and  simultaneously  trimming  the  excess 
soil.  The  top  and  bottom  surfaces  of  the  soil  in  the  ring 
were  trimmed  to  produce  uniform  horizontal  planes.  Gaps  were 
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left  between  the  soil  and  the  top  and  bottom  edges  of  the 
ring  to  allow  for  proper  seating  of  the  porous  plates. 

The  oedometer  ring  was  weighed  before  and  after 
emplacement  of  the  soil  sample.  The  dimensions  of  the  soil 
sample  were  measured  and,  together  with  the  weight  of  soil, 
used  to  calculate  the  initial  total  density.  The  trimmings 
were  collected  for  index  testing  and  approximate  water 
content  determinations.  The  entire  trimming  operation  took 
approximately  one  hour.  Undoubtedly,  some  drying  of  the 
sample  occurred  during  preparation. 

4.2.3  Testing  Procedure 

The  oedometer  ring  containing  the  sample  was  taken 
immediately  to  the  testing  laboratory  after  preparation. 
Saturated  filter  papers  were  placed  in  contact  with  the  top 
and  bottom  soil  surfaces.  Porous  plates  which  had  been 
saturated  in  distilled  boiling  water  were  then  placed 
adjacent  to  the  filter  papers.  This  assemblage  was  placed  in 
the  oedometer  apparatus. 

The  enclosure  around  the  oedometer  ring  was  filled  with 
distilled  water  and  the  lever  arm  was  positioned  on  the 
sample.  The  LVDT  was  threaded  into  the  lever  arm  and  a  small 
load  was  applied  to  the  sample  to  obtain  the  initial  LVDT 
reading.  The  first  stress  increment  was  applied  and  the 
output  LVDT  voltage  was  recorded  with  time.  In  the  first 
five  minutes  the  readings  were  recorded  manually  because  the 
minimum  time  interval  between  readings  for  the  data 
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acquisition  system  was  one  minute.  Thereafter,  the  data 
acquisition  system  was  set  to  record  the  LVDT  output  voltage 
at  predetermined  time  intervals. 

Loads  were  added  to  the  samples  in  all  three  tests  in 
increment  ratios,  Ad,/ 0, ,  equal  to  one.  The  only  exceptions 
were  the  last  two  increments  for  sample  Cl  where  the 
ratio  was  0.5.  An  increment  ratio  of  one  ensured 
that  a  classical  Terzaghi  consolidation  curve  was  obtained 
which  allowed  an  unambiguous  determination  of  the  end  of 
primary  consolidation.  Mesri  and  Godlewski  (1977)  have 
discussed  the  influence  of  increment  ratio  on  the  shape  of 
the  compression-logarithm  of  time  curve  as  have  Leonards  and 
Girault  (1961). 

The  stress  at  each  increment  was  left  for  a  minimum  of 
one  day  and,  generally,  not  longer  than  one  week  before  the 
next  stress  increment  was  applied.  All  stress  increments 
were  maintained  for  a  sufficient  time  period  so  that  the 
secondary  compression  index  could  be  adequately  defined. 
Increment  durations  were  generally  longer  for  stresses  above 
the  preconsolidation  stress  where  the  creep  characteristics 
were  of  primary  interest. 

A  computer  program  was  written  to  reduce  the  data  and 
present  the  void  ratio  versus  logarithm  of  time  results  on  a 
computer  graphics  facility.  This  permitted  daily  visual 
monitoring  of  the  test  progress  and  enabled  a  quick  response 
to  any  malfunctions. 
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At  the  end  of  testing,  the  loads  were  removed  in  one 
step  without  monitoring  the  swelling  of  the  sample.  The 
sample  was  immediately  taken  from  the  oedometer  and  weighed 
for  a  water  content  determination. 


4.3  Triaxial  Creep  Tests 

4.3.1  General 

A  total  of  six  incremental  drained  creep  triaxial  tests 
were  undertaken  on  Tar  Island  clay.  Three  separate  triaxial 
systems  were  used  which  allowed  creep  tests  to  be  performed 
simultaneously  at  three  different  consolidation  pressures. 
The  first  suite  of  three  tests  used  vertical  samples  with 
horizontal  bedding  planes.  The  second  suite  used  inclined 
samples  which  had  the  bedding  plane  at  45°.  In  all  cases  the 
isotropic  consolidation  pressure  during  shear  was  in  excess 
of  the  apparent  preconsolidation  stress  determined  from  the 
oedometer  testing. 

The  triaxial  testing  equipment  and  procedures  were 
designed  to  determine  the  influence  of  anisotropy  and 
consolidation  pressure  or  water  content  on  the  creep 
behaviour  of  the  clay. 
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4.3.2  Triaxial  Testing  Equipment 


4.3.2. 1  General 

Three  identical  separate  triaxial  systems  were 
fabricated.  Each  system  consisted  of  a  triaxial  cell,  a 
hanging  yoke  with  no  mechanical  advantage  for  adding 
weights,  a  system  for  applying  cell  pressure,  a  system  for 
applying  back  pressure,  an  LVDT  to  monitor  vertical 
displacement  and  a  burette  to  provide  volume  change 
measurements.  All  triaxial  tests  were  performed  in  the  same 
laboratory  as  the  oedometer  tests.  The  room  temperature  was 
monitored  daily  and  the  results  are  given  in  Appendix  B. 

The  design  of  the  triaxial  systems  followed 
conventional  standards  and  only  the  general  aspects  of 
design  and  operation,  as  these  affect  tests  results,  will  be 
discussed  in  the  following  sections. 


4. 3. 2. 2  Triaxial  Cells 

The  cells  were  made  of  aluminum  and  were  designed  for  a 
maximum  pressure  of  3500  kPa.  The  inside  diameter  of  the 
cells  was  approximately  13  cm  and  the  inside  height  was 
approximately  23  cm.  The  pedestal  on  the  cell  base  was 
designed  for  a  sample  diameter  of  3  cm  but,  with  a  change  of 
pedestels,  the  cells  are  capable  of  housing  7.5  cm-0  x  15  cm 
samples  or  3.8  cm-0  x  7.6  cm  samples  with  lateral  strain 
gauges.  The  cell  is  shown  in  Plate  4.1 


Plate  4.1  Triaxial  cell 


Plate  4.2  Sample  with  membrane 
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The  top  plate  of  the  cell  was  connected  to  the  cylinder 
and  the  cylinder  base  plate  by  six  steel  rods.  The  pressure 
seal  at  the  cylinder/top  plate  contact  was  provided  by  an 
O-ring.  The  cylinder  arrangement  was  connected  to  the  base 
plate  by  six  steel  threaded  studs  with  wing  nuts  for  easy 
removal.  The  pressure  seal  at  the  cylinder/base  plate 
contact  was  again  provided  by  an  O-ring. 

The  base  plate  had  ports  which  allowed  access  to  the 
cell  interior  for  sample  drainage,  measurement  of  sample 
pore  pressure,  cell  pressure  and  electronics.  The 
electronics  port  can  be  used  for  a  lateral  strain  gauge  or 
an  interior  load  cell.  The  top  platechad  a  port  which  was 
used  for  filling  and  emptying  the  cell.  The  cell  fluid  was 
emptied  through  the  cell  pressure  port  by  applying  air 
pressure . 

The  sample  arrangement  in  the  cell  interior  is  shown  in 
Plate  4.2.  The  sample  sat  on  filter  paper  and  a  porous 
plate.  The  top  cap  rested  directly  on  the  sample.  Drainage 
was  permitted  from  the  sample  base  only  but  top  drainage  can 
be  added  with  only  minor  modifications  to  the  cell.  The 
sample  size  of  3  cm-0  x  6  cm  used  for  testing  was  not 
standard  so  custom  membranes  were  manufactured  at  the 
University  of  Alberta  according  to  Madeiros  (1979).  The 
membrane  was  seated  to  the  top  cap  and  bottom  pedestel  by 
two  O-rings  at  each  position.  The  O-rings  fitted  into 
grooves  and  were  tightened  against  the  membrane  with  hose 
clamps.  A  correction  to  the  vertical  stress  on  the  sample 
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due  to  membrane  strength  was  not  made  because  it  is 
generally  considered  to  be  small  (see  Bishop  and  Henkel, 
1962).  Membrane  leakage  was  never  observed.  This  is  due,  in 
part,  to  the  use  of  silicone  as  a  cell  fluid  as  recommended 
by  Poulos  (1964)  and  Lovenbury  (1969). 

Load  was  applied  to  the  sample  by  a  ram  through  the 
cell  top.  The  base  of  the  ram  was  machined  to  fit  a  ball 
bearing  which  fitted  into  the  sample  top  cap.  This 
arrangement  ensured  that  the  ram  applies  load  without 
eccentricity.  An  O-ring  and  an  oil  seal  provided  the 
pressure  seal  for  the  ram/cell  top  interface.  This  seal 
arrangement  resulted  in  a  reduction  of  external  load  due  to 
ram  friction.  The  magnitude  of  the  ram  friction  was 
determined  for  each  cell  before  testing. 

The  magnitude  of  the  ram  friction  was  calculated  by 
comparing  the  internal  load  measured  with  a  load  cell  and 
the  external  load  applied  by  the  dead  weights  on  the  yoke.  A 
typical  set  of  results  for  a  variety  of  external  loads  and 
cell  pressures  is  presented  in  Figure  4.2.  All  curves  are 
parallel  to  the  curve  passing  through  the  origin  which 
incicates  that  ram  friction  is  independent  of  the  magnitude 
of  the  external  load.  The  intercept  of  the  curves  with  the 
ordinate  is  equal  to  the  load  on  the  ram  end  from  cell 
pressure  plus  ram  friction.  By  subtracting  the  former  it  was 
found  that  ram  friction  varied  from  9  nt  to  18  nt  and  was 
independent  of  cell  pressure.  A  common  value  of  18  nt  for 
ram  friction  was  adopted  for  all  three  cells.  It  should  be 


Figure  4.2  Typical  ram  friction  test  results 
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mentioned  that  this  value  of  ram  friction  is  for  static 
conditions  while  the  magnitude  of  the  kinematic  friction  is 
probably  more  applicable  to  the  testing  conditions. 

4. 3. 2. 3  Pressure  Systems 

A  schematic  of  the  cell  fluid  pressure  system  for  each 
cell  is  shown  in  Figure  4.3.  The  cell  pressure  was  applied 
by  means  of  two  separate  systems. One  was  a  low  pressure 
system,  0  to  1380  kPa,  that  originated  from  a  compressed  air 
supply.  A  Moore  continuous  bleed  regulator  was  used  to  set 
the  desired  pressure  and  an  aluminum  closed  end  cylinder 
served  as  the  compressed  air/cell  fluid  interface.  This 
system  was  capable  of  maintaining  pressure  within  ±3  kPa 
over  sustained  periods  of  time.  A  backup  generator  was 
available  for  the  compressed  air  pump  in  the  event  of  a 
power  outage.  Plate  4.3  shows  the  aluminum  plate  which 
supports  the  regulator  and  interface.  The  back  pressure 
system  for  the  sample  pore  pressure  was  identical  to  this 
system  and  is  shown  in  Figure  4.4. 

For  cell  pressures  above  1380  kPa  a  high  pressure 
system  was  fabricated  which  was  similar  in  detail  to  other 
high  pressure  systems  used  in  the  geotechnical  laboratories 
at  the  University  of  Alberta.  It  is  illustrated 
schematically  in  Figure  4.3  and  a  portion  is  shown  in  Plate 
4.2.  The  pressure  source  was  a  17,200  kPa  nitrogen  bottle.  A 
simple  step  regulator  attached  to  the  nipple  of  the  bottle 
was  used  to  decrease  the  bottle  pressure  to  a  lower, safer 
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Mohr  Regulator 


Hoke  Regulator 


Plate  4.3  Cell  pressure  board 


Plate  4.4  Volume  change  indicators  and  back  pressure  board 
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level  of  3500  kPa  for  this  testing  program.  Downstream  of 
the  step  regulator  was  a  Hoke  920  regulator  which  was  used 
to  set  the  desired  cell  pressure.  An  EMG  Hydraulics  lnc.r 
hydropneumatic  accumulator  with  a  sliding  rubber  bladder 
served  as  the  nitrogen/cell  fluid  interface.  The  purpose  of 
the  rubber  bladder  was  to  isolate  the  explosive  gas  pressure 
from  any  system  leaks.  The  interior  of  the  accumulator  will 
corrode  if  in  contact  with  water  so  another  fluid  interface 
is  required  downstream  of  the  accumulator  if  rin  the 
future , water  is  chosen  as  the  cell  fluid. 

There  are  some  drawbacks  to  the  high  pressure  system 
just  described.  It  was  difficult  to  set  the  Hoke  regulator 
at  the  desired  pressure.  The  chosen  pressure  inevitably  fell 
by  about  10%  (in  approximately  12  hours)  before  it 
stabilized.  It  is  considered  that  this  was  a  result  of  air 
going  into  solution.  Furthermore,  the  nitrogen  bottles  lost 
pressure  due  to  inevitable  leaks  in  the  system.  The  cell 
pressure  dropped  by  approximately  7  kPa  per  week.  The 
bottles  had  to  be  changed  every  few  weeks,  taking  care  not 
to  disturb  an  ongoing  creep  test.  A  example  plot  of  pressure 
versus  time  for  the  backpressure  system  and  the  high 
pressure  system  is  included  in  Appendix  B. 

Sweglok  brass  fittings  were  used  for  the  entire  system. 
These  fittings  have  a  maximum  capacity  of  3500  kPa.  If 
pressures  above  3500  kPa  are  required  for  future  testing  all 
brass  fittings  must  be  replaced  with  stainless  steel 
fittings.  Whitey  ball  values  were  used  to  direct  flow  for 
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the  entire  system.  An  advantage  of  ball  valves  over  stem 

valves  is  that  there  is  no  volume  change  associated  with 
their  operation. 

Copper  tubing  with  a  diameter  of  0.32  cm  was  used  for 
the  cell  pressure  system  and  a  combination  of  copper  tubing 

and  polyethylene  tubing  was  used  for  the  back  pressure 
system. 

A  7000  kPa  transducer  was  used  to  monitor  cell  pressure 
and  a  2100  kPa  transducer  was  used  to  monitor  back  pressure. 

4. 3. 2. 4  Monitoring  of  Sample  Deformation 

The  axial  deformation  of  the  sample  was  measured  with  a 
Hewlett-Packard  24DCT-250  LVDT  attached  to  the  ram  as  shown 
in  Plate  4.1.  The  LVDT’s  were  calibrated  before  and  after 
testing.  No  significant  change  in  the  calibration  factors 
were  observed.  The  LVDT’s  had  a  linear  travel  of  ±0.16  cm 
with  a  concomitant  voltage  output  of  ±3.7500  volts. This 
corresponds  to  a  theoretical  sensitivity  of  ±0.000004  cm. 
The  actual  sensitivity  observed  was  closer  to  ±0.00004  cm. 

The  same  Fluke  data  acquisition  which  was  used  for  the 
oedometer  testing  was  also  used  for  this  testing.  The  data 
acquisition  unit  was  used  to  monitor  the  LVDT  voltage 
excitation  and  output  as  well  as  the  output  voltage  of  the 
two  pressure  transducers. 

The  volume  change  of  the  sample  was  measured  with  a 
burette  which  had  an  apparent  sensitivity  of  ±0.005  cc. 
Figure  4.4  presents  a  schematic  of  the  back  pressure/volume 
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system  and  Plate  4.4  shows  the  volume  change  indicators. 
Each  volume  change  indicator  consisted  of  an  outer  perspex 
cylinder  with  aluminum  end  caps  which  housed  a  10  cc  glass 
burette.  The  pressure  seal  at  the  aluminum/perspex  contact 
was  provided  by  an  O-ring.  A  dyed  kerosene/water  interface 
recorded  the  fluid  movement  in  the  burette.  The  inside  of 
the  burette  had  to  be  kerosene  wet  before  assembly  or  the 
meniscus  would  stick  to  the  glass. 

Tests  were  performed  to  determine  the  leakage  from  the 
back  pressure/volume  change  system.  An  aluminum  cylinder 
encased  in  a  membrane  in  place  of  the  sample  was  used  for 
these  tests.  Back  pressure  and  cell  pressure  were  applied 
and  it  was  found  that  the  system  leakage  varied  between  0.05 
and  0.1  cc/day.  The  possible  sources  of  leakage  are 
considered  to  be  : 

1.  Leakage  through  valves  and  Sweglok  fittings 

2.  Diffusion  through  and  creep  of  the  polyethylene  .32  cm  0 
tubing  (see  Poulos,  1964) 

3.  Incomplete  saturation  of  the  system 

Each  component  of  the  system  was  checked  and  replaced 
if  faulty.  For  the  second  suite  of  creep  tests  most  of  the 
polyethylene  tubing  was  replaced  by  copper  tubing  which 
decreased  the  leakage  rate  by  50%.  The  leakage  rate  was 
found  to  be  a  direct  function  of  the  magnitude  of  the  back 
pressure.  However,  cell  pressure  had  no  effect  on  the 
leakage  rate  which  indicated  the  sample  membranes  did  not 
leak.  Slight  changes  in  room  temperature  caused  substantial 
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changes  in  the  burette  readings. 

It  is  considered  that  the  leakage  rate  was  reduced  to 
the  minimum  possible  for  a  system  of  this  design.  Further 
minor  reductions  could  be  possible  if  the  number  of  valves 
and  fittings  and  the  length  of  tubing  were  reduced.  Tavenas 
et.al  (1978)  did  this  by  incorporating  a  burette  in  the  cell 
base.  Lovenbury  (1969)  undertook  a  series  of  long  term 
drained  creep  tests  and  found  leakage  rates  of  between  0.005 
and  0.015  cc/ day  for  similar  paraffin  burettes.  He  also 
measured  volume  change  with  a  mercury  jacket  around  the 
sample.  These  measurements  were  more  accurate  but 
reliability  problems  hampered  their  long  term  efficacy. 

The  problem  with  measuring  volume  change  with  burettes 
was  recognized  at  the  beginning  of  the  testing  program.  It 
was  considered  that  a  leakage  rate  correction  could  be 
applied  to  the  volume  change  readings.  This  view  proved  to 
be  optimistic  as  the  correction  was  usually  greater  than  the 
actual  volume  change  of  the  sample  within  one  day  after 
beginning  a  creep  increment.  Consequently,  volumetric  strain 
with  time  cannot  be  reported  in  this  thesis.  For  the  sample 
size  used  in  this  testing  program  a  leakage  rate  of  less 
than  0.001  cc/day  was  required  before  volumetric  strain  with 
time  could  be  determined  reliably  for  any  reasonable  period 
of  time.  However,  it  was  possible  to  estimate  the  total 
volumetric  strain  during  a  creep  increment  which  was  used  in 
the  calculation  of  applied  stress. 
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It  is  of  benefit  for  future  long  term  drained  creep 
research  to  review  the  volume  change  devices  for  triaxial 
testing  described  in  the  literature  and  to  assess  their 
applicability  to  creep  testing.  The  devices  are  in  two 
categories;  those  that  directly  measure  pore  fluid  exchange 
with  the  sample  and  those  that  measure  the  change  in 
diameter  of  the  sample,  usually  referred  to  as  lateral 
strain  gauges. 

The  first  category  consists  mainly  of  burettes  with 
various  adaptations  for  continuous  automatic  data  logging.  A 
few  systems  are  described  in  the  literature  where  fluid 
exchange  with  the  sample  caused  a  shift  in  a  mercury  system 
which  was  automatically  recorded.  For  descriptions  of 
various  types  of  volume  change  devices  see  Lewin  (1971), 
Mitchell  and  Burn  (1971),  Rowlands  (1972),  Klementev  (1974), 
Menzies  (1975),  Sharpe  (1978)  and  Watts  (1980).  The 
sensitivity  of  the  instruments  described  by  these  authors 
ranged  from  ±0.002  cc  to  ±0.5  cc.  None  of  these  authors 
mentioned  leakage  rates,  presumably  because  the  instruments 
were  designed  for  short  term  testing.  With  only  a  minimum 
leakage  rate  superposed  on  the  maximum  sensitivity  of  the 
instruments  it  is  apparent  that  none  would  be  suitable  for 
long  term  drained  creep  tests  of  the  type  undertaken  in  this 
testing  program. 

The  alternative  to  volume  change  indicators  are  lateral 
strain  gauges.  The  most  common  are  adaptations  of  that 
described  by  Bishop  and  Henkel  (1962).  Modern  versions  such 
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as  those  commonly  used  at  the  University  of  Alberta  employ  a 
vertically  mounted  LVDT  to  measure  displacement.  Menzies 
(1976)  describes  an  adaptation  which  employs  a  horizontally 
mounted  LVDT.  An  attempt  was  made  to  design  a  gauge  for  this 
testing  program  for  use  with  3  cm-0  samples  but  problems 
with  mounting  the  LVDT  on  a  small  gauge  and  bedding  errors 
caused  the  attempt  to  be  abandoned.  In  general,  lateral 
strain  gauges  have  sufficient  accuracy  to  be  used  for  long 
term  drained  creep  tests. 

El-Ruway ih  (1976)  described  a  lateral  strain  gauge 
developed  at  Imperial  College  which  is  superior  to  the 
conventional  design.  It  consisted  of  a  spring  stainless 
steel  strip  which  can  be  shaped  to  fit  any  sample  size. 
Lateral  displacement  was  measured  with  four  strain  gauges 
mounted  on  the  strip.  One  advantage  over  the  conventional 
design  is  that  several  gauges  can  be  mounted  on  one  sample. 
The  device  has  high  sensitivity  and  is  very  linear  but  it 
appears  difficult  to  calibrate. 

In  recent  years  a  lateral  strain  gauge  which  does  not 
contact  the  sample  has  been  described  in  the  literature  (see 
Cole,  1978  and  Kahn  and  Hoag,  1979).  The  device  consists  of 
a  source  transducer  and  an  aluminum  target  on  the  sample.  It 
operates  on  the  eddy  current  loss  principle.  The  device  has 
a  sensitivity  of  ±2.54  x  10“4mm,  is  simple  to  calibrate, 
insensitive  to  immersion  in  most  liquids  and  is  operable 
over  a  temperature  range  of  at  least  10°C  to  -6.7°C.  It  is 
recommended  that  the  utility  of  this  device  be  investigated 
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for  long  term  drained  creep  tests. 

4. 3. 2. 5  Suggested  Modifications  To  Equipment 

With  the  exception  of  the  volume  change  indicators  the 
testing  equipment  performed  as  anticipated.  However,  several 
modifications  are  suggested  here  which  would  improve  the 
test  results: 

1.  A  noncontacting  lateral  strain  gauge  should  be  adapted 
for  use  with  this  triaxial  cell  design. 

2.  The  data  acquisition  system  should  be  upgraded  to  a  unit 
which  allows  the  time  interval  between  readings  to  be 
set  for  each  channel  individually  .  This  would  preclude 

o 

having  to  filter  small  time  interval  readings  necessary 
for  one  creep  test  from  another  creep  test  which  has 
been  running  for  an  extended  period. 

3.  Thermistors  should  be  added  to  each  cell  so  that  the 
temperature  of  the  cell  fluid  can  be  continually 
monitored  with  a  data  acquisition  unit. 

4.  A  more  reliable  high  pressure  system  must  be  developed 
for  long  term  testing. 

5.  The  dead  weight  loading  system  should  be  replaced  by  a 
bellofram  so  that  axial  loads  can  be  added  more 
prec i sely . 

6.  The  transducer  manifold  system  should  be  replaced  by 
individual  transducers  for  each  cell  so  pressure  can  be 
continuously  monitored. 


92 


4.3.3  Sample  Preparation 

Test  specimens  for  triaxial  testing  were  prepared  in  a 
device  which  allowed  rotation  of  the  sample  and  is  shown  in 
Plates  4.5  and  4.6.  Both  vertical  and  inclined  samples  were 
prepared.  The  samples  were  prepared  by  rotating  and  trimming 
away  excess  material  with  a  sharp  knife.  Three  vertical 
samples  could  be  trimmed  from  one  7  cm  long  section  of  the 
shelby  tube  samples.  The  diameter  of  the  shelby  tube  samples 
was  7.3  cm.  However  it  required  23  cm  of  shelby  tube  sample 
to  prepare  one  sample  inclined  at  45°  to  the  axis  of  the 
shelby  tube. 

A  section  of  a  shelby  tube  sample  cut  at  45°  before 
final  preparation  of  an  inclined  sample  is  shown  in  Plate 
4.5.  The  maximum  diameter  of  an  inclined  sample  which  could 
be  prepared  from  a  shelby  tube  sample  of  the  f orement ioned 
diameter  an  This  constraint  determined  the  maximum  sample 
size  of  3  cm-0  x  6  cm  for  all  triaxial  testing. 

After  trimming  the  samples  were  carefully  measured  and 
weighed.  Trimmings  were  used  for  index  tests  and  water 
content  determinations. 

4.3.4  Testing  Procedure 
4 . 3 . 4 . 1  General 

Each  drained  triaxial  creep  test  involved  three 
distinct  stages.  A  B  test  was  performed  as  the  first  stage 
to  select  the  back  pressure  necessary  to  saturate  the 
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Plate  4.5  Inclined  sample  before  trimming 


Plate  4.6  Inclined  sample  after  trimming 
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sample.  In  the  second  stage  the  sample  was  subjected  to 
isotropic  consolidation  in  steps  to  the  desired 
consolidation  pressure.  Finally  the  sample  was  sheared  at 
constant  cell  pressure  by  incrementally  applying  axial  loads 
which  were  maintained  for  various  lengths  of  time.  This 
latter  step  constitutes  an  incremental  creep  test.  The  axial 
deformation  and  volume  change  were  measured  versus  time  at 
each  creep  increment.  At  the  conclusion  of  the  last  creep 
increment  the  sample  was  sheared  to  failure  by  adding  loads 
slowly  enough  to  ensure  drained  failure.  The  stress  path 
just  described ,  63  constant-o,  increasing,  was  the  only 
stress  path  studied  during  the  laboratory  investigation. 

With  the  available  equipment  three  samples  could  be 
tested  concurrently.  A  suite  of  three  vertical  samples  was 
tested  first.  The  second  and  final  suite  of  tests  used 
inclined  samples.  The  vertical  samples  were  tested  first 
because  they  could  be  replaced  if  the  test  malfunctioned  for 
any  reason.  Due  to  the  shortage  of  shelby  tube  samples  the 
inclined  samples  could  not  be  replaced  if  the  test 
malfunctioned.  Modifications  to  the  testing  technique  as  a 
result  of  experience  gained  from  the  first  suite  of  tests 
led  to  improved  test  procedures  for  the  second  suite  of 
tests . 

A  package  of  computer  programs,  developed  especially 
for  this  testing  program,  were  used  to  monitor  the  progress 
of  each  test  daily. 
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4 . 3 . 4 . 2  B  Tests 

The  samples  were  saturated  with  back  pressure  at  the 
start  of  each  test.  The  back  pressure  was  applied  under  a 
slightly  higher  cell  pressure  for  24  hours  to  allow  the 
sample  to  saturate  before  performing  a  B  test.  The  B  tests 
were  performed  according  to  conventional  standards.  The  cell 
pressure  was  increased  by  70  kPa  and  the  pore  pressure  was 
recorded  with  time  until  it  stabilized.  This  usually  took 
about  five  minutes.  The  calculated  value  of  B  was  usually  in 
the  range  of  0.92  to  0.98.  It  is  considered  that  all  the 
samples  were  saturated  and  the  B  parameter  was  marginally 
less  than  one  because  of  sample  stiffness. 

4. 3. 4. 3  Isotropic  Consolidation 

Each  sample  was  consolidated  in  steps  to  the  desired 
cell  pressure  with  cell  pressure  increments,  a<5  3/03  ,  equal 
to  one.  In  the  first  suite  of  samples  the  volume  change  with 
time  was  monitored  using  the  burettes  only.  This  proved  to 
be  unsatisfactory  for  volumetric  strain  rate  calculations 
because  of  burette  leakage,  although  reliable  estimates  of 
the  total  volume  change  could  be  made. 

For  the  second  suite  of  samples,  both  volume  change  and 
axial  displacement  with  time  were  measured.  The  axial 
displacement  was  monitored  with  an  LVDT  connected  to  the 
ram.  The  ram  was  seated  to  the  sample  by  an  axial  load  which 
was  just  in  excess  the  ram  friction  and  the  load  on  the  ram 
end  due  to  cell  pressure.  The  ram  had  to  be  reset  with  the 


. 
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sample  drainage  line  closed  each  time  the  consolidation 
pressure  was  increased. 

4. 3. 4. 4  Incremental  Creep  Tests 

Axial  loads  were  added  to  the  sample  in  steps.  At  each 
step  or  increment  vertical  displacement  and  volume  change 
were  measured  with  time.  Three  increments,  referred  to  here 
as  extended  increments,  were  chosen  for  each  test  where  the 
stress  was  maintained  constant  for  periods  in  excess  of 
10,000  mins.  Since  the  cross  sectional  area  of  the  sample 
continuously  increased  during  shear  loading  small  loads  were 
added  to  the  sample  to  maintain  constant  stress  for  the 

extended  increments  only.  The  changing  area  was  calculated 
according  to: 

A  =  A0  ( 1 ~  c  3 ) 2  4,1 

where  A  =  area 

Ac  =  sample  cross  sectional  area  at  the 
beginning  of  present  increment 
e3  =  lateral  strain  (always  negative  for  this 
stress  path) 

The  lateral  strain  is  given  by: 


£  3  = 


ev-e  ,/2 


4.2 


MM*  t  ■  >  ^ •  .1 
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where  -  volumetric  strain  (estimated  from  volume 

change  readings,  based  on  sample  volume  at 
the  beginning  of  the  increment) 
c,  =  axial  strain  (based  on  sample  height  at 
the  beginning  of  the  increment) 

Equation  4.1  is  valid  only  if  it  is  assumed  that  the 
soil  is  cross-anisotropic.  Equation  4.1  can  be  used  directly 
with  a  lateral  strain  gauge.  Despite  the  approximations  used 
to  estimate  ,  the  test  results  indicate  that  the  stress 
was  maintained  constant  during  the  extended  increments. 

Axial  loads  were  added  to  the  sample  in  small  enough 
increments  to  preclude  undrained  failure.  The  increase  in 
axial  stress  at  any  stress  level  just  necessary  to  cause 
undrained  failure  is  given  by: 


A6,  =  6\ ( 1-sin0f )-6'3 ( 1+sin0' )/( 1-2A)sin0’-1 )  4.3 


where  a, 

o; 


<5; 


0’ 

A 


=  change  in  axial  stress  necessary  to  cause 
undrained  failure 

=  effective  axial  stress  at  the  beginning 
of  the  increment 

=  effective  radial  stress  at  the  beginning 
of  the  increment 

=  effective  angle  of  shearing  resistance 
=  pore  pressure  parameter 


The  derivation  of  Equation  4.3  is  given  in  Appendix  B. 


. 
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Equation  4.3  can  be  extended  to  include  cohesion.  In  order 
to  use  Equation  4.3  it  is  necessary  to  know  0'  and  the 
variation  of  A  with  stress  level.  The  value  of  0'  for  this 
testing  program  was  estimated  initially  from  strain 
controlled  shear  tests  by  Hardy  Associates  (1978)  Ltd  (see 
Design  Review  Panel,  1977).  The  value  of  A  was  estimated  for 
each  stress  level .Equation  4.3  was  used  as  a  check  at  the 
beginning  of  each  stress  increment  to  ensure  that  the 
applied  load  would  not  cause  undrained  failure.  It  can 
easily  be  shown  from  Equation  4.3  and  a  simple  consideration 
of  the  stress  path  that  Aot ,  the  change  in  axial  load 
necessary  to  cause  undrained  fa i lure , decreases  with 
increasing  axial  load. 

Loads  were  maintained  for  periods  of  1500  to  5000  mins 
at  the  intervening  increments  between  the  extended 
increments.  The  stress  was  not  maintained  constant  at  the 
intervening  increments.  At  the  conclusion  of  the  third 
extended  increment  loads  were  added  until  the  sample  failed. 
These  loads  were  added  slowly  enough  so  that  the  failure  was 
completely  drained. 

The  creep  test  just  described  is  referred  to  as  an 
incremental  creep  test.  A  possible  alternative  to  this 
procedure  is  the  single  increment  creep  test  where  the 
desired  load  is  added  in  one  step.  This  type  of  test  is 
commonly  employed  in  undrained  creep  testing  of  soils  and 
creep  testing  of  rock.  A  single  increment  creep  test  cannot 
be  used  in  a  drained  creep  examination  of  normally 
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consolidated  soils  at  higher  stress  levels  because  the 
sample  would  fail  in  undrained  shear  due  to  the  generation 
of  pore  pressures  before  exhibiting  drained  creep.  The 
influence  of  the  incremental  testing  procedure  on  the  test 
results  will  be  discussed  in  Chapter  5. 


5.  Experimental  Results 


5.  1  Introduction 

The  results  from  the  oedometer  test  program  ,the 
triaxial  test  program  and  an  investigation  of  the  soil 
composition  are  presented  in  this  chapter.  The  testing 
program  was  conducted  primarily  to  investigate  the  creep 
behaviour  of  the  clay.  Other  information  regarding  shear 
strength  and  compressibility  was  also  acquired  during  the 
course  of  the  laboratory  program  and  is  detailed  herein.  It 
will  be  shown  in  Chapter  6  that  the  relevant  state  of  stress 
in  the  portion  of  the  clay  deposit  to  be  analysed  is  one  of 
simple  shear.  Simple  shear  laboratory  tests  were  not 
conducted  because  Wu  et.al  and  Edgers  et.a/  have  that  there 
are  only  minor  differences  between  the  creep  parameters  from 
triaxial  and  simple  shear  tests. 

It  is  demonstrated  that  the  creep  test  results  from  the 
triaxial  program  fit  a  Singh-Mi tchell  relationship.  The 
influence  of  stress  state,  octahedral  stress,  test  procedure 
and  anisotropy  on  the  creep  parameters  in  the  Singh-Mi tchell 
equation  are  discussed.  The  applicability  of  the  laboratory 
results  to  the  field  problem  under  consideration  is  reviewed 
in  Chapter  6. 
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5.2  Soil  Composition 

5.2.1  Index  Tests 

The  specimens  of  Tar  Island  clay  were  very  stiff, 
medium  to  highly  plastic  and  contained  varying  amounts  of 
silt,  sand,  organic  material,  coal  fragments  and  calcium 
carbonate  shell  material.  A  series  of  standard  index  tests 
were  performed  on  each  sample  and  the  results  are  presented 

in  Table  5.1.  A  grain  size  analysis  for  each  sample  is 
included  in  Appendix  B. 

Horizontal  bedding  was  evident  throughout  the  entire 
soil  horizon  sampled.  Bedding  planes  were  marked  by  organic 
material  concentrated  in  thin  dark  seams,  less  than  .10  cm 
thick  at  irregular  intervals,  and  by  variations  in  material 
properties.  The  organic  material  also  occurred  disseminated 
throughout  the  soil.  The  bedding  planes  are  evident  in  Plate 
C.3  in  Appendix  C. 

The  plasticity  index  ranged  from  20  to  30  percent  and 
the  liquid  limit  varied  from  44  to  54  percent.  The  average 
natural  water  content  of  the  soil  tested  was  approximately 
27  percent  which  is  significantly  less  than  the  natural 
water  content  of  approximately  45  percent  before  dyke 
construction.  The  difference  in  water  content,  of  course, 
reflects  consolidation  of  the  clay  due  to  the  imposed  dyke 
load. 
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Table  5 . 1 


Summary  of  index  test  results 


Sample 

No. 

Depth 

(m) 

w-% 

P 

(  i  ) 

Density 
(mg/m3 ) 

%-Silt 

%-Clay 

PL-% 

(  2  ) 

LL-% 

(  3  ) 

Cl 

61.3 

28.2 

2.68 

2.00 

62 

34 

C2 

60.0 

30.3 

2.71 

1.98 

60 

38 

- 

- 

C3 

60.5 

32.8 

2.66 

1.89 

63 

32 

— 

— 

T 1 

61.5 

26.0 

2.74 

2.02 

60 

28 

23.  1 

44.6 

T2 

61.5 

26.4 

2.74 

1.96 

60 

28 

23.  1 

44.6 

T4 

59.6 

25.9 

2.64 

1.99 

62 

34 

23.5 

50.0 

I  1 

59.  1 

24.2 

2.70 

1.98 

66 

28 

23.9 

49.8 

12 

59.0 

26.5 

2.69 

1.94 

58 

36 

24.3 

53.7 

13 

60.  1 

25.0 

2.71 

1.99 

58 

34 

23.4 

51.0 

C-oedometer  test  samples 
T,I-triaxial  test  samples 

(1  ^specific  gravity 
( 2  5 -plastic  1 imi t 
( 3  *  -liquid  limit 
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The  average  specific  gravity  of  the  soil  tested  was 
2.70  and  the  average  density  was  1.97  Mg/m3.  Approximately 
30%  of  the  clay  was. clay  size  material  with  the  majority  of 
the  remainder  in  the  silt  size  range. 

Although  sampling  was  undertaken  only  in  one  location 
over  a  limited  depth,  a  comparison  with  the  index  tests  in 
Chapter  2  indicate  that  the  samples  tested  adequately 
represent  the  average  properties  of  the  entire  clay  deposit. 

5.2.2  Mineralogy 

An  X-ray  diffraction  analysis  was  performed  by  the 
Alberta  Research  Council  to  determine  the  mineralogy  of  the 
clay  size  fraction.  The  majority  of  the  clay  size  fraction 
is  quartz  and  the  remainder  is  composed  of  clay  minerals. 
The  clay  mineral  assemblage  consists  of  40  percent  illite,35 
percent  montmor illoni te , 20  percent  kaolinite  and  5  percent 
chlorite.  The  relative  percentage  of  the  clay  minerals  is 
based  on  the  'peak  height  method'.  From  the  scanning 
electron  microscope  analysis  described  in  the  next  section 
the  particles  larger  than  2  microns  are  predominantly 
quartz . 

Two  bulk  samples  were  also  examined  for  organic 
content ,  pH  and  CaC03 .  The  laboratory  procedures  used  to 
determine  these  quantities  are  described  by  McKeague  (1978). 
The  organic  content  is  approximately  1.9%  by  weight,  from 
one  test  only,  although  it  can  be  locally  much  higher 
because  of  concentration  in  thin  seams.  The  pH  of  the  soil 
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is  7.9,  an  average  of  two  tests. 

The  total  percentage  of  CaC03  by  weight  in  the  soil  is 
4  percent,  an  average  of  two  tests.  The  quantitative 
gravimetric  method  which  consists  of  measuring  the  amount  of 
C02  liberated  after  pouring  HC1  on  the  soil  was  used  to 
determine  the  CaC03  content.  It  is  not  possible  to 
distinguish  among  CaC03  in  solution  in  the  pore  water,  CaC03 
precipitated  on  soil  particles,  perhaps  contributing  to 
cementation  of  the  soil,  and  CaC03  in  the  form  of  fossil 
shells  by  this  method.  Because  of  the  consequences  of  soil 
cementation  on  the  creep  and  shear  strength  behaviour,  a 
scanning  electron  microscope  study  was  undertaken  in  cin 
attempt  to  distinguish  the  form  CaC03  takes  in  the  soil. 

5.2.3  Scanning  Electron  Microscope  Analysis 

A  Cambridge  Stereoscan  150  scanning  electron  microscope 
(SEM)  operated  by  the  Department  of  Entomology  at  the 
University  of  Alberta  was  used  to  examine  the  microfabric  of 
the  soil.  A  Kevex  elemental  analytical  mass  spectrometer  was 
used  in  conjunction  with  the  SEM  to  determine  the 
predominate  elements  in  the  soil.  The  spectrometer  has  the 
capability  of  analysing  the  entire  screen  of  the  SEM  or  any 
selected  spot  on  the  screen. 

The  samples  for  the  SEM  were  obtained  by  cutting 
approximately  1  cm  cubes  from  a  fresh  fracture  surface  that 
was  produced  by  manually  breaking  a  shelby  tube  sample. 
Samples  were  obtained  from  both  horizontal  and  vertical 
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fracture  surfaces.  Sample  moisture  was  removed  by  immersing 
the  samples  in  liquid  nitrogen  at  195. 8°C,  then  liquid  freon 
and  finally  sublimating  the  sample  in  a  dessicator. 

The  samples  were  examined  at  magnifications  ranging 
from  25X  to  10,000X.  Plate  5.1  illustrates  the  fabric  of  the 
soil  at  different  magnifications.  The  elemental  analysis 
shown  on  the  spectrometer  charts  are  for  the  SEM  micrographs 
which  show  the  4  micron  scale  designation. 

The  low  magnification  micrographs  reveal  that  the  soil 
is  dense  with  little  pore  space.  At  higher  magnifications  it 
is  evident  that  silica  grains  predominate.  Clay  platelets 
occur  in  groups  or  masses  surrounding  the  silica  grains.  In 
a  number  of  micrographs  silica  grains  are  joined  by 
connectors  which  may  be  bonds  (see  Collins  and  McGown, 
1974) . 

The  mass  spectrometer  results  reveal  that  the 
predominate  element  in  the  soil  is  silica  with  lesser 
amounts  of  aluminum,  iron,  calcium  and  potassium.  The 
primary  sources  of  the  silica  are  the  particle  grains  as 
confirmed  by  several  spot  analysis.  The  mineralogy  of  the 
clay  mineral  assemblage  accounts  for  most  of  the  elements 
present.  Aluminum  and  silica  are  present  in  all  the  clay 
minerals.  Iron  and  potassium  are  constituents  of  illite  and 
chlorite.  Calcium  may  occur  in  the  montmor i lion i te .  Spot 
analyses  with  the  mass  spectrometer  revealed  that  iron  and 
calcium  are  constituents  of  most  of  the  connectors  between 
soil  grains. These  connectors  may  provide  some  cohesion  but 
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Mass  Spectrometer  Results 


Plate  5.1  Scanning  electron  micrographs 


107 


this  can  only  be  confirmed  with  strength  tests. 


5.3  Oedometer  Test  Results 


5.3.1  General 

Oedometer  tests  were  undertaken  as  a  first  step  in  the 
laboratory  creep  study  to  investigate  the  creep  behaviour  of 
the  soil  in  a  simple  way  and,  secondly,  to  ascertain  some 
conventional  engineering  parameters  of  the  soil  such  as  the 
coefficient  of  consolidation  and  the  preconsolidation 
stress.  The  conventional  analysis  of  the  oedometer  test  data 
is  given  in  this  section  and  the  creep  behaviour  of  the  soil 
is  discussed  in  the  subsequent  two  sections. 

Table  5.2  summarizes  the  oedometer  test  results  for 
sample  Cl .  Figure  5. 1  presents  the  void  ratio  versus 
logarithm  of  time  curves  for  each  stress  increment  of  sample 
Cl.  The  test  results  for  samples  C2  and  C3  are  presented  in 
similar  fashion  in  Appendix  B. 

In  order  to  compile  the  void  ratio  versus  time  curves 
as  shown  in  Figure  5.1,  the  initial  void  ratio  was 
calculated  from  the  initial  water  content  and  the  specific 
gravity.  The  initial  water  content  and  the  specific  gravity 
was  determined  on  the  basis  of  the  entire  sample,  not  on  the 
preparation  trimmings.  The  void  ratio  curves  were  then 
calculated  from  the  displacement  readings,  allowing  for  the 
compliance  of  the  equipment. 
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Table  5.2  Oedometer  test  results,  sample  Cl 


ncr . 
No. 

Stress 

(kPa) 

Elapsed 

Time 

(min ) 

Void 

Ratio 

Change 

Ae 

Cv<”  Ca 

( cm!/sec ) 

m( 2 ’  Unit 

Strai 
Rate ( 3 
( %/min 

1 

47.9 

3093 

.0113 

— 

.0004 

1  .  14 

.049 

2 

98.2 

2230 

.0096 

- 

.0009 

1.08 

.064 

3 

194.0 

2763 

.0144 

- 

.0010 

1.04 

.065 

4 

386.4 

2890 

.0240 

- 

.0015 

1.03 

.095 

5 

770.6 

2563 

.0382 

.0024 

.0027 

1.05 

.171 

6 

1512 

9018 

.0628 

.  0022 

.0038 

1.08 

.307 

7 

2997 

5387 

.0754 

.0013 

.0047 

1.06 

.333 

8 

4482 

4301 

.0472 

- 

.0051 

1.08 

.322 

9 

6000 

6707 

.0316 

- 

.  0054 

1.01 

.273 

( 1 ’Calculated  only  at  stresses  above  the  preconsolidation 
stress  where  possible 

( 2 ’Slope  of  logarithm  of  volumetric  strain  versus  logarithm 
of  time  curves 

( 3 ’Intercept  of  logarithm  of  volumetric  strain  versus 
logarithm  of  time  at  t=1  min,  ignoring  consolidation  effects 
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Figure  5.1  Void  ratio  versus  logarithm  of  time 
, sample  Cl 
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Void  ratio  versus  logarithm  of  vertical  effective 
stress  curves,  e-log  p  curves,  have  been  compiled  for  all 
three  samples  in  order  to  determine  the  compression  index, 
Cc,  and  the  preconsolidation  stress,  pc.  The  results  for 
sample  Cl  are  given  in  Figure  5.2  and  the  results  for 
samples  C2  and  C3  are  included  in  Appendix  B. 

Out  of  interest  the  e-log  p  curves  are  plotted  for 
three  different  time  intervals;  that  corresponding  to  the 
end  of  primary  consolidation,  the  actual  duration  of  the 
test  which  is  a  variable  and  , finally, an  extrapolated  time 
of  10,000  minutes.  The  value  of  Cc  is  shown  on  the  plots  and 
for  the  longer  time  curves  ranges  from  0.25  to  0.31  with  an 
average  of  0.28.  The  value  of  pc  is  quite  sensitive  to  the 
time  interval  and  ranges  from  480  kPa  to  900  kPa  with  a  mean 
of  650  kPa.  Since  the  clay  is  normally  consolidated  the  mean 
value  should  be  close  to  the  effective  stress  imposed  on  the 
soil  by  the  dyke  load  at  the  time  and  depth  of  sampling  . 

The  coefficient  of  consolidation,  Cv,  is  tabulated  for 
stresses  above  the  preconsolidation  stress  of  650  kPa  on 
Table  5.2  for  sample  Cl  and  on  similar  tables  for  samples  C2 
and  C3  in  Appendix  B.  The  magnitude  of  Cv  ranges  from  0.0024 
cm2/sec  to  0.0001  cm2/sec  and  generally  decreases  with 
increasing  stress  above  the  preconsolidation  stress. 
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Figure  5.2  Void  ratio  versus  vertical  stress , sample  Cl 
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5.3.2  Coefficient  of  Secondary  Consolidation 

It  is  possible  to  analyse  the  creep  behaviour  of  a  soil 
from  laboratory  test  results  in  a  number  of  ways.  The  most 
common  for  oedometer  test  results  is  to  determine  the 
coefficient  of  secondary  consolidation,  Co,  according  to 
Equation  3.11  and  plot  it  against  the  vertical  stress.  This 
section  is  devoted  to  an  analysis  of  this  type.  Another 
possible  analytical  route  is  to  plot  the  logarithm  of 
vertical  strain  rate  versus  the  logarithm  of  time  and 
proceed  with  an  analysis  according  to  Singh  and  Mitchell 

(1968).  The  next  section  is  devoted  to  this  alternate 
approach. 

For  each  void  ratio  versus  logarithm  of  time  curve  Ca 
was  determined  graphically  assuming  that  the  slope  of  the 
curve  is  independent  of  time  in  the  secondary  range.  Figure 
5.3  shows  the  variation  of  Ca  with  vertical  stress  for  all 
three  samples.  There  is  some  tendency  for  Ca  to  level  off 
with  increasing  stress,  especially  where  the  load  increment 
ratio  drops  below  one.  This  is  consistent  with  other  work, 
such  as  that  by  Mesri  and  Godlewski  (1977)  and  Muryama  and 
Shibata  (1964)  who  found  that  Ca  was  a  constant  at  stresses 
above  one  to  two  times  the  preconsolidation  stress. 

Mesri  (1973)  presented  a  global  correlation  between 
water  content  and  the  coefficient  of  secondary  compression, 
an  analogous  form  of  Ca,  which  is  given  here  as  Figure  5.4. 
The  oedometer  test  results  are  plotted  on  Figure  5.4  and  fit 
his  correlation  quite  well.  According  to  Mesri's  ( op.cit .) 
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Figure  5.4  Coefficient  of  secondary  compression  for  natural 
soil  deposits , after  Mesri,1973 
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creep  classification  system  for  the  oedometer  stress  state, 
the  clay  test  here  has  a  'low'  secondary  compressibility. 

5.3.3  Vertical  Strain  Rate 

As  mentioned  previously  the  oedometer  data  can  also  be 
analysed  in  terms  of  the  logarithm  of  vertical  strain  rate 
versus  the  logarithm  of  time.  The  results  from  the  test  on 
sample  C2  are  plotted  in  this  manner  in  Figure  5.5.  Similar 
plots  are  included  in  Appendix  B  for  samples  Cl  and  C3.  The 
vertical  strain  rate  was  calculated,  after  da  Fontoura 
(1980),  using  a  moving  five  point  linear  regression  analysis 
on  the  strain  versus  time  data.  The  vertical  strain  was 
calculated  on  the  basis  of  the  sample  height  at  the 
beginning  of  each  increment.  Analysing  the  data  in  this 
fashion  is  not  new.  it  was  first  done  by  Wilson  and  Lo 

(1965)  and,  later,  by  Lovenbury  (1969)  and  Tavenas  et.al 
(  1978)  . 

Each  strain  rate  versus  time  curve  in  Figure  5.5  can  be 
broken  into  two  sections  based  on  the  shape  of  the  curve  and 
time.  In  the  first  section  which  is  at  times  less  than  t100, 
as  determined  from  the  void  ratio  versus  time  curves,  the 
strain  rate  continually  decreases  with  time  in  a  nonlinear 
manner.  The  time  dependent  deformation  is  a  result  of  both 
transient  consolidation  and  creep  effects. 

In  the  second  section  of  the  curve  at  times  greater 
than  t, oo  the  strain  rate  decreases  continually  with  time 
but  in  a  linear  fashion.  The  deformation  is  considered  to  be 
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Figure  5.5  Vertical  strain  rate  versus  time  , sample  C2 


entirely  a  result  of  creep  effects.  A  linear  regression 

analysis  was  used  to  best  fit  a  line  through  the  data  points 
of  this  part  of  the  curve. 

The  slope  of  these  best  fit  lines  is  equivalent  to  the 
creep  parameter,  m,  in  Equation  3.2.  The  value  of  m  for  all 
three  samples  ranges  from  1.01  to  1.29  with  a  mean  of  1.12 
(standard  deviat ion=0 . 07 ) .  The  values  of  m  for  each  test  for 
each  sample  are  included  in  the  summary  tables  for  the 
oedometer  tests  and  are  plotted  in  Figure  5.11.  The 
variation  in  m  is  small  and  is  attributable  to  sample 
variations.  From  this  data  it  is  apparent  that  m  is 
independent  of  vertical  stress  and  load  increment  ratio.  It 
is  assumed  that  m  is  independent  of  the  duration  of  the 
transient  consolidation  effects. 

It  is  obvious  from  Figure  5.5  that  an  increase  in  the 
vertical  stress  causes  an  increase  in  the  strain  rate.  To 
investigate  the  influence  of  vertical  stress,  the  strain 
rate  at  time  equals  one  minute  (unit  strain  rate)  for  the 
best  fit  lines  has  been  plotted  against  vertical  stress  in 
Figure  5.6.  This  figure  shows  that  the  unit  strain  rate 
increases  with  stress  at  lower  stresses  and  then  appears  to 
level  off  at  higher  stresses  especially  when  the  load 
increment  ratio  drops  below  one.  It  is  assumed  that 
transient  consolidation  does  not  effect  the  position  of  the 
strain  rate  curves,  that  is,  the  magnitude  of  the  unit 
strain  rate.  The  best  fit  lines  in  Figure  5.6  can  be 
described  by  a  relationship  of  the  Singh-Mi tchell  type  as 
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follows : 


where  ?, 
A 


t 

m 


m 

*1  =  A(1/t)  5.1 

=  vertical  or  volumetric  strain  rate 
=  unit  strain  rate,  a  function  of  vertical 
stress  and  load  increment  ratio 
=  t  ime 

=  1.12,  a  constant 


The  difference  between  Equations  3.2  and  5.1  is  the 
lack  of  a  shear  stress  level  term  in  Equation  5.1.  This  is  a 
result  of  the  stress  state  in  the  oedometer  test.  Regardless 
of  the  vertical  stress  the  shear  stress  level  is  a  constant 
in  the  oedometer  and  it's  effect  is  therefore  included  in  A. 
It  follows  from  this  that  the  dependence  of  m,  if  any,  on 
shear  stress  level  cannot  be  evaluated  with  oedometer  tests. 

It  must  also  be  added  that  since  m  is  not  precisely 
equal  to  one  then  the  assumption  of  Ca  being  independent  of 
time  in  Section  5.3.2  is  incorrect.  Ca  decreases  with  time. 


5.4  Triaxial  Creep  Test  Results 

5.4.1  General 

The  results  of  two  suites  of  drained  triaxial  creep 
tests  are  presented  in  this  section.  Both  suites  consist  of 
three  samples,  each  consolidated  under  different  isotropic 
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consolidation  pressures.  The  samples  in  one  suite  were 
prepared  with  the  bedding  plane  horizontal  and  are  numbered 

T 1 '  T2  and  T4*  The  samples  in  the  other  suite  were  prepard 
with  the  bedding  plane  inclined  at  an  angle  of  45°  and  are 
numbered  II,  12  and  13.  The  suite  of  inclined  samples  was 
tested  at  45°  to  investigate  the  anisotropy  of  the  creep 
behaviour  but  also  to  have  the  maximum  shearing  stress  along 
the  bedding  plane.  The  relevance  of  this  stress  orientation 
to  the  field  problem  is  discussed  in  Chapter  6. 

All  samples  were  consolidated  isotropically  to 
pressures  above  the  average  preconsolidation  stress 
determined  from  the  oedometer  testing  program.  The  stress 
range  below  the  preconsolidation  stress  was  not  investigated 
in  shear  since  only  the  normally  consolidated  creep 
behaviour  of  clays  is  of  interest  here  and  differs  from  the 
overconsolidated  behaviour.  The  final  isotropic 
consolidation  pressures  during  shear  were  varied  from  740 
kPa  to  1500  kPa  between  samples  which  allowed  an 
investigation  of  the  influence  of  normal  effective  stress  or 
water  content  on  the  creep  behaviour.  This  is  of  importance 
in  a  field  problem  of  this  type  where  the  effective  stresses 
vary  substantially  in  the  lateral  and  vertical  directions. 

The  test  results  from  the  isotropic  consolidation  stage 
are  presented  in  the  next  section.  The  results  from  the 
incremental  creep  tests  are  then  presented.  Following  is  a 
discussion  of  the  various  factors  which  influence  the  creep 
parameters  in  the  Singh-Mi tchell  equation  derived  from  this 
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testing  program. 

5.4.2  Isotropic  Consolidation 

All  samples  were  consolidated  isotropically  in  pressure 
increments,  AO;/o3,of  one  to  the  selected  pressure  for 
testing.  It  is  generally  considered  that  a  sample 
consolidated  anisotropically  along  the  same  stress  path  it 
has  experienced  in  the  field  will  yield  results  that  are 
more  indicative  of  the  field  behavior  than  an  isotropically 
consolidated  sample.  For  this  field  situation  consolidating 
anisotropically  would  have  brought  the  sample  to  the  60  % 
stress  level  at  the  beginning  of  the  incremental  creep  test 
which  would  have  precluded  the  study  of  the  creep  behaviour 
below  this  stress  level.  Since  the  proper  definition  of  the 
stress  level  term  in  the  Singh-Mi tchell  equation  requires 
testing  at  stress  levels  below  60  %,  the  samples  were 
consolidated  isotropically. 

The  isotropic  consolidation  results  for  the  vertical 
and  inclined  samples  are  presented  in  Tables  5.3  and  5.4 
respectively.  The  vertical  strain  and  vertical  strain  rate 
versus  time  curves  are  presented  in  Figure  5.7  for  the  last 
consolidation  step  of  sample  12.  Similar  plots  are  contained 
in  Appendix  C  for  the  remainder  of  the  consolidation  steps 
for  sample  12  and  the  corresponding  curves  for  samples  II 
and  13.  The  volume  changes  for  the  inclined  samples  given  in 
Table  5.4  are  calculated  from  the  axial  strain  readings 
assuming  that  the  sample  is  isotropic.  Axial  displacement 
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Table  5.3 


Summary  of  isotropic  consolidation 
vertical  samples 


results , 


Sample  o; 
No. 

( kPa ) 


Volume  Height  Volume  Time 
Change 

(cm3)(,)  (cm)  (cm3)  (min) 


B  Test 
Result 

(  4  ) 


T 1 


T2 


T4 


0 

0 

6.39 

165 

-0.9 

6.35 

559 

-1.29 

6.29 

1  103 

-1.09 

6.24 

0 

0 

5.83 

(  5  ) 

+  5 . 07 (  2  > 

6.07 

190 

-3.25 

5.93 

380 

-0.54 

5.91 

758 

-0.79 

5.87 

0 

0 

6.15 

187 

+  0.27 

6.16 

377 

-0.60 

6.  13 

722 

-0.74 

6.10 

1500 

-1.38 

6.04 

0(  3  ) 

+  1.15 

6.09 

827 

-1.51 

6.02 

1507 

-0.85 

5.98 

45.17  -  0.97 

44.27  1176 

42.97  5592 

41.88  4328 

40.93  -  0.98 

46.00  1116 

42.75  5523 

42.21  4433 

41.42  4450 

44.34  -  0.97 

44.61  4312 

44.01  2740 

43.28  1441 

41.89  4074 

43.04  613 

41.53  1485 

40.70  5734 


(  1  5  volume  change  readings  are  calculated  including  a  leakage 
correction 

( 2  5  burette  seriously  malfunctioned,  replaced,  reading  only 
approximate 

(3)cell  pressure  reduced  to  zero  due  to  leakage  of  cell 
fluid  through  a  crack  in  cell  top,  replaced,  sample 
reconsolidated 

(  4 ’ backpressure  is  207  kPa  for  all  samples 
(  5  Sample  swelled  at  0  kPa 


' 
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Table  5.4 


Summary  of  isotropic  consolidation 
inclined  samples 


results , 


Sample 

No. 

6; 

( kPa ) 

Volume 
Change 
(cm3)'  1  5 

Height 

(cm) 

Volume 

(cm3 ) 

Time 

(min ) 

B  Test 
Result 

1 1 

0 

0 

6.70 

48.63 

— 

0.92 

365 

1  .  14 

6.65 

47.49 

5790 

743 

0.85 

6.61 

46.64 

4017 

12 

0 

0 

6.81 

49.86 

— 

0.98 

334 

1.05 

6.76 

48.81 

2600 

672 

0.77 

6.72 

48.04 

4145 

- 

1386 

1 .80 

6.64 

46.24 

4549 

13 

0 

0 

6.39 

47.35 

- 

0.96 

271 

2.28 

6.29 

45.08 

1405 

548 

0.85 

6.25 

44.23 

1370 

109  0 

1.25 

6.19 

42.98 

4539 

I  1  Calculated  from  axial  strain  assuming  deformation  is 
isotropic 

( 2  Absolute  value  of  slope  of  logrithm  of  axial  strain  rate 
versus  logarithm  of  time 

( 3  backpressure  is  310.3  kPa  for  13  and  413.7  kPa  for  II  and 
12 


Vertical  Strain  Rate,  %/min 
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Figure  5.7  Vertical  strain  and  strain  rate  versus  logarithm 
of  time  , sample  I2f  consolidation  pressure= 1 386  kPa 
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readings  were  not  taken  with  time  for  the  vertical  samples. 

The  volume  changes  given  in  Table  5.3  for  the  vertical 

samples  are  from  the  burette  readings  and  have  been 

corrected  for  leakage. 

The  axial  strain  in  Figure  5.7  is  calculated  on  the 
basis  of  the  sample  height  at  the  beginning  of  the 
consolidation  increment.  The  strain  rate  curve  in  Figure  5.7 
shows  the  same  characteristic  shape  as  the  oedometer  strain 
rate  curves.  The  creep  portion  of  the  strain  rate  curves  for 
all  the  inclined  samples  have  much  steeper  slope  than  the 
oedometer  samples.  Deformation  due  to  creep  under  isotropic 
stresses  appears  to  cease  with  time  especially  at  the  lower 
consolidation  pressures. 

5.4.3  Incremental  Creep 

The  results  of  the  incremental  creep  tests  are 
presented  in  this  section.  A  summary  of  the  test  results  for 
sample  T4  is  given  in  Table  5.5.  Vertical  strain  versus  time 
for  each  creep  increment  of  sample  T4  is  given  in  Figure  5.8 
.  The  axial  strain  is  calculated  on  the  basis  of  the  sample 
height  at  the  end  of  isotropic  consolidation.  The  stress 
levels  included  in  Table  5.5  and  on  Figure  5.8  are  the  ratio 
of  the  shear  stress  during  the  increment  divided  by  the 
shear  stress  at  failure.  Similar  summary  plots  and  tables  to 
those  presented  for  sample  T4  are  contained  in  Appendix  C 
for  all  the  other  samples. 
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Table  5.5  Summary  of  incremental  creep  test  results, 

sample  T4 


Incr . 
No. 

Deviator 
Initial 
( kPa ) 

Stress 
Final 
( kPa ) 

Stress 

Level 

(%) 

Time 

(min ) 

Axial 

Strain 

<%) 

Cumulative 
Axial  Strain 
(%) ( 2 ’ 

1 

350 

350 

13.2 

1303 

0.34 

0.34 

2 

781 

777 

29.4 

14541 

1.32 

1.66 

3 

1070 

1068 

40.3 

1331 

0.43 

2.09 

4 

1365 

1358 

51.4 

12747 

1.82 

3.91 

5 

1727 

1717 

65.  1 

1200 

1  .  14 

5.05 

6 

1982 

1947 

74.7 

2821 

2.94 

7.99 

7 

2147 

2143 

80.9 

14538 

2.23 

10.22 

8 

2328 

2312 

87.7 

4194 

1.02 

11.24 

9 

2434 

2395 

91.7 

7340 

2.26 

13.50 

10 

2550 

2485 

96.  1 

3317 

3.48 

16.98 

1  1 

2582 

2539 

97.3 

1121 

2.24 

19.22 

1  2  (  1  ’ 

2654 

2619 

100.0 

130 

1.73 

20.95 

Total  duration  of  test  =  64583  min,  c5 ;  =  1 507  kPa 


( 1 ’Failure 

( 2 ’Strain  at  the  end  of  the  increment 


Vertical  Strain,  % 


Figure  5.8  Vertical  strain  versus  time  , sample  T4 
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Figure  5.9  Vertical  strain  rate  versus  time  , sample 
I  3-left , sample  t4-right) 
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The  strain  rate  has  been  calculated  for  the  incremental 
creep  results  in  the  same  manner  as  for  the  oedometer  test 
results.  A  plot  of  the  logarithm  of  strain  rate  versus  the 
logarithm  of  time  is  shown  in  Figure  5.9  for  samples  T4  and 
13.  The  strain  rate  results  for  the  extended  increments  only 
are  shown  on  Figure  5.9.  Stress  levels  for  the  extended 
increments  are  also  indicated  on  Figure  5.9.  The  strain  rate 
plots  for  the  other  samples  are  included  in  Appendix  C. 

The  strain  rate  plot  in  Figure  5.9  has  the  same 
characteristics  as  the  strain  rate  plot  for  the  oedometer 
test  shown  in  Figure  5.5.  The  strain  rate  decreases 
continuously  with  time  for  the  length  of  the  test.  The 
strain  rate  for  the  initial  portion  of  the  curves  decreases 
continuously  but  in  a  nonlinear  fashion.  This  strain  is  a 
result  of  both  consolidation  and  creep  effects.  In  the 
latter  portion  of  the  curve  the  strain  rate  decreases 
continuously  in  a  linear  fashion  and  is  entirely  due  to 
drained  creep. 

Other  authors ,  notably  Lovenbury  (1969)  and  Tavenas 
et.al  (1978)  have  also  recorded  transient  consolidaton 
effects  on  strain  rate  versus  time  plots  for  drained 
triaxial  creep  tests.  As  an  example,  the  strain  rate  versus 
time  for  Pancone  Clay  from  Lovenbury  ( op.cit .  )  is 
illustrated  in  Figure  5.10.  A  strain  instability  at 
approximately  40  to  100  days  is  clearly  evident  on  Figure 
5.10.  It  is  considered  that  the  duration  of  the  transient 
consolidation  effect  depends  upon  the  coefficient 
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Figure  5.10  Strain  rate  versus  time  for  Pancone  clay  (after 
Lovenbury,  1969) 
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consolidation,  Cv,  the  sample  drainage  conditions,  the 
sample  size  and  the  load  increment  ratio,  i<5,/6,.  A  drained 
creep  incremental  test  should  be  designed  so  that  each  of 

these  factors  is  manipulated  to  minimize  the  duration  of  the 
transient  effect. 

For  the  purposes  of  the  interpretation  of  the  test 
results  reported  herein  it  is  assumed  that  the  duration  of 
the  transient  effect  does  not  influence  the  slope  of  the 
linear  portion  of  the  strain  rate  versus  time  curves  or  the 
position  of  the  curves,  that  is,  the  magnitude  of  the  strain 
rate  for  any  given  time  after  the  transient  effect.  To  this 
end  a  best  fit  straight  line  based  on  a  linear  regression 
analysis  has  been  drawn  through  the  linear  portion  of  the 

curves  and  extrapolated  to  lesser  times  ignoring 
the  transient  effect. 

A  summary  of  the  creep  parameters  for  the  strain  rate 
plots  for  all  samples  is  given  on  Table  5.6.  The  slope  of 
the  lines,  m,  is  plotted  against  octahedral  stress,  the 
first  stress  invariant,  for  all  triaxial  samples  and  the 

i 

oedometer  samples  in  Figure  5.11.  The  average  value  of  m  for 
all  tests  is  1.13.  The  unit  strain  rate  for  all  triaxial 
samples  is  plotted  versus  stress  level  in  Figure  5.12.  It  is 
observed  that  the  slope  of  the  curves  are  similar  with  the 
exception  of  sample  12  and  average  approximately  3.0.  The 
average  value  of  the  unit  strain  rate  at  zero  stress  level 
is  approximately  0.18  %/min  with  a  range  of  0.54  %/min  to 
0.04  %/min. 


Table  5.6  Summary  of  incremental 
samples 
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creep  test  results, all 


Sample 

No. 

Stress 

Level 

(%) 

<5oct 

( kPa ) 

m 

Unit  Strain 
Rate 
( %/min ) 

Duration  of 
I ncrement 
(min ) 

T 1 

23.6 

1446 

1.15 

0.28 

16742 

49.9 

1827 

1.19 

0.84 

15754 

67.4 

2081 

1.19 

1.8 

27255 

T2 

26.7 

880 

1.09 

0.11 

17171 

53.8 

1004 

1.10 

0.33 

17418 

76.9 

1110 

1.11 

0.70 

204  18 

T4 

29.4 

1767 

1.15 

0.44 

14541 

51.4 

1962 

1.11 

0.59 

12747 

80.9 

2223 

1.10 

1.0 

14538 

I  1 

36.3 

881 

1.15 

0.25 

6800 

68.1 

1001 

1.21 

0.60 

1431 

83.4 

1059 

1.21 

0.88 

13189 

12 

34.2 

1618 

1  .26 

1 .14 

7103 

63.3 

1816 

1.20 

1.80 

10180 

88.0 

1984 

1.20 

5.20 

8845 

13 

32.4 

1277 

1.00 

0.14 

7123 

61.4 

1438 

1  .06 

0.45 

10110 

84.2 

1564 

1.06 

0.70 

9915 

' 
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Figure  5.11  Creep  parameter,  m,  versus  octahedral  stress 
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Figure  5.12  Unit  strain  rate  versus  stress 
tests 
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From  Figures  5.11  and  5.12  it  is  apparent  that  the 
creep  behaviour  of  this  clay  fits  a  Singh-Mi tchell 
relationship  because: 

1.  The  strain  rate  versus  time  decreases  linearly  on  a 
In-In  plot. 

2.  The  creep  parameter,  m,  is  independent  of  stress  level. 

3.  The  strain  rate  increases  exponentially  with  stress 
level  for  a  given  time. 

If  the  average  parameters  from  all  tests  are  used  the 
vertical  strain  rate  versus  time  can  be  described  by  a 
relationship  of  the  form: 

?,  =  .  18eiB(  1/t )  3  5.2 
Although  Equation  5.2  holds  for  the  average,  it  is  of  value 
to  discuss  the  factors  which  influence  the  creep  parameters, 
A,  a  and  m  in  Equation  5.2  in  more  detail. 

5.4.4  Factors  Influencing  Creep  Parameters 
5.4.4. 1  Stress  State 

The  influence  of  stress  state  on  the  creep  parameters 
can  be  evaluated  from  a  comparison  of  the  oedometer  and 
triaxial  test  results.  The  value  of  the  creep  parameter,  m, 
as  shown  on  Figure  5.11  is  insensitive  to  stress  state.  The 
unit  strain  rate  for  the  oedometer  test  results  are  plotted 
on  Figure  5.12  for  an  assumed  stress  level  for  the  oedometer 
stress  state  and  this  soil  of  61  percent.  Again  the 
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oedometer  results  are  approximately  similar  to  the  triaxial 
test  results. 

It  appears  from  this  limited  comparison  that  the 
influence  of  stress  state  on  the  creep  parameters  is  small 
compared  to  sample  variations.  This  conclusion  cannot  be 
extrapolated  to  other  stress  states  such  as  plane  strain 
simple  shear  without  additional  testing  however.  The 
influence  of  stress  state  will  be  discussed  further  in 
Chapter  6. 


5. 4. 4. 2  Anisotropy 

The  average  values  of  m  for  the  vertical  and  inclined 
samples  are  1.13  and  1.15  respectively.  The  average  values 
of  a  for  the  vertical  and  inclined  samples  are  3.18  and  2.74 
respectively.  It  is  apparent  that  any  effects  due  to 
anisotropy  are  masked  by  sample  variation.  Any  effect  of  the 
possible  cohesion  discussed  in  section  5.2.3  does  not 
manifest  itself  in  the  creep  behaviour  of  the  clay. 


5. 4. 4. 3  Octahedral  Stress 

The  value  of  the  parameter,  m,  as  shown 
is  insensitive  to  octahedral  stress  for  both 
and  oedometer  tests.  The  value  of  the  stress 
for  the  triaxial  tests  does  not  show 
relationship  with  isotropic  consolidation  str 


in  Figure  5.11 
the  triaxial 
level  term,  a, 
a  consistent 
ess  for  either 


the  inclined  or  vertical  samples.  The  same  is  true  for  the 
parameter  A.  The  equivalent  parameter  A  in  the  oedometer 
tests  shows  some  tendency  to  a  constant  value  above  the 
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preconsolidation  stress. 

It  may  therefore  be  stated  that  the  creep  parameters  in 
Equation  5.2  are  not  influenced  by  the  magnitude  of  the 
octahedral  stress. 

5. 4. 4. 4  Testing  Procedure 

The  effect  of  the  incremental  testing  procedure  on  the 
creep  parameters  in  Equation  5.2  cannot  be  definitely  stated 
since  it  is  not  possible  to  compare  the  results  with  single 
increment  tests.  As  discussed  in  Chapter  4  single  increment 
tests  cannot  be  used  to  examine  drained  creep  at  high  stress 
levels  since  the  soil  will  fail  in  undrained  shear  before 
drainage  occurs.  An  advantage  of  an  incremental  creep  test 
is  that  the  creep  behaviour  of  a  single  soil  sample  can  be 
examined  at  a  variety  of  stress  levels.  The  influence  of  the 
variation  in  properties  between  soil  samples  is  thus 
eliminated. 

The  key  variable  which  can  be  used  to  manipulate  the 
results  from  an  incremental  creep  test  is  the  stress 
increment  ratio,  AO,/o,  .  The  magnitude  of  ao,/o,  is  governed 
by  Equation  4.3  and  must  decrease  with  increasing  c, .  The 
ratio,  A<5,/0,  ,  governs  the  duration  of  transient 
consolidation,  the  relative  importance  of  previous  creep 
strains  on  the  strains  during  the  present  increment  and  the 
influence  of  pseudo-preconsolidation  on  the  strains  in  the 
succeeding  increments. 
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Transient  Consolidation 

Tavenas  et.al  (1978)  in  a  comprehensive  series  of 
drained  creep  tests  postulate  that  the  stress  increment 
ratio  must  be  kept  less  than  0.3  in  order  to  diminish  the 
relative  effect  of  primary  consolidation  strains  to  creep 
strains.  The  effect  of  consolidation  was  probably  not 
pronounced  for  the  results  of  Tavenas  et.al  ( op.Cit .) 
because  the  soil  tested,  St.  Alban  clay,  is  overconsolidated 
at  the  low  stresses,  less  than  50  kPa,  used  during  their 
testing  program.  It  is  considered  that  decreasing  the 
drainage  path  is  more  effective  than  manipulation  of  stress 
increment  ratio  in  reducing  the  effect  of  consolidation. 

Previous  Creep  Strains 

It  is  not  clear  how  the  previous  creep  strains  affect 
the  creep  strains  at  succeeding  increments.  Mitchell  et.al 
(1969)  have  examined  this  problem  for  undrained  creep  tests 
on  soil  and  have  related  the  creep  strain  rates  between 
increments  by  a  relationship  of  the  form: 

?(D,,t  )/?(D,,t  )  =  (t  /t  )m+(t  /t  -t  ’  )  5.3 


a  =  In (Y+ 1 )/ (D2”D, ) 


5.4 


m 

=  (  (tj  -t1  )/U  )  [  *  (D2  ,tj  )  /  *  ( D , 


m 

1 1  -  <  1 1  /tj  >  ] 


Y 
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5.5 

The  variables  in 

Equation  5.3  are  defined  in 

Figure  5.13. 

Equation 

5.3 

was 

derived  from  Equation 

3.2  with  the 

assumption 

that 

the 

net  strain  in  the  second 

increment  is 

independent 

of 

the 

accumulated  strain  in 

the  previous 

increment . 

Using 

this  equation  they  find  that  the  creep 

parameter , 

a , 

is 

the  same  for  multiple  increment  tests  as 

for  single  increment  tests. 


Pseudo- Preconsolidation 

As  a  sample  experiences  drained  creep  the  effective 
stress  state  remains  constant  but  the  yield  surface  must 
move  away  from  the  effective  stress  point  of  the  sample 
because  of  the  decrease  in  volume  of  the  sample.  The  stress 
increment  ratio  must  be  high  enough  to  overcome  this  pseudo 
preconsolidation  or  the  soil  may  behave  in  an 
overconsolidated  manner.  This  effect  can  be  reduced  by 
limiting  the  duration  of  the  intervening  stress  increments 
between  the  extended  increments  to  the  minimum  necessary  to 
dissipate  the  pore  pressure. 


5.5  Shear  Strength  Results 

5.5.1  General 

Although  the  primary  objective  of  the  incremental  tests 
was  to  evaluate  the  creep  behaviour  of  the  clay,  the  results 
also  yield  the  stress-strain  and  strength  properties  of  the 


Strain 
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Figure  5.13  Strain  versus  time  relat ionship-f or  step-creep 
test  (after  Mitchell  et.ai,  1969) 
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clay.  The  stress  strain  behaviour  is  discussed  in  the  next 
section.  This  is  followed  by  a  discussion  of  the  strength 
properties  which  includes  a  comparison  with  the  strain 
controlled  test  results  from  Hardy  Associates  (1978)  Ltd. 
(see  Design  Review  Panel,  1977). 


5.5.2  Stress-Strain  Behaviour 

The  results  from  the  incremental  tests,  as  shown  in 
Figure  5.7,  can  also  be  plotted  in  the  conventional  manner 
of  deviatoric  stress  versus  strain.  Such  a  plot  is  shown  in 
Figure  5.13  for  sample  13.  Similar  plots  for  all  other  tests 
are  included  in  Appendix  C. 

The  stress  and  strain  at  the  beginning  and  the  end  of 
each  increment  on  Figure  5.13  are  shown  and  these  points  are 
joined  by  straight  lines.  This  straight  line  path  represents 
the  stress  strain  path  followed  by  the  sample  during  the 
test.  A  smooth  line  is  also  drawn  through  the  final 
stress-strain  point  of  each  extended  increment.  The  smooth 
line  represents  the  long  term  stress  strain  behaviour  of  the 
clay . 

The  stress  strain  curves  demonstrate  that  the  material 
is  ductile  with  strains  to  failure  in  excess  of  10  percent. 
The  clay  continually  compresses  during  shear.  The  volume 
change  curve  levels  off  at  higher  stresses. 

All  samples  failed  along  distinct  shear  planes  and  as 
with  all  stress  controlled  tests,  failure  was  catastrophic. 
It  is  known  that  the  clay  is  strain  softening  (see  Design 


aV/V,  %  Deviator  Stress  ( cr1  -  03 ),  kPa 


Figure  5.14  Stress-strain  curve  , sample  13 
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Review  Panel,  1977)  but  it  was  not  possible  to  evalute  this 
characteristic  in  a  stress  controlled  test.  Photographs  of 
all  failed  specimens  are  included  in  Appendix  C.  The 
inclined  specimens  failed  along  bedding  and,  usually,  along 

obvious  organic  layers.  The  vertical  samples  failed  across 
bedding . 

It  must  be  stated  the  sample  T2  was  not  taken  to 
failure  although  the  final  vertical  strain  approached  14  %. 
It  is  estimated  that  the  final  imposed  stress  is  within  5  % 
of  the  failure  stress. 

5.5.3  Strength  Envelope 

A  p-q  plot  of  the  shear  strength  results  are  presented 
in  Figure  5.15  for  both  the  creep  tests  and  the  strain 
controlled  tests  from  Hardy  Associate  Ltd.  The  effective 
angle  of  shearing  resistance,  0’,  for  all  tests  is  26.5°. 
The  effective  cohesion  intercept,  c?,  is  zero. These  results 
are  based  on  a  best  fit  line  through  all  of  the  data  points 
shown  on  Figure  5. 15.  It  is  apparent  that  the  shear  strength 
results  determined  from  this  testing  program  are  consistent 
with  previous  work. 

An  angle  of  24°  will  be  used  for  the  computations  in 
Chapter  6  as  this  value  is  indicative  of  the  shear  strength 
for  the  portion  of  the  soil  to  be  analysed. 


■  • 
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Figure  5.15  P-q  plot,  from  creep  tests  and  previous  work 
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5.6  Concluding  Remarks 

The  creep  behaviour  of  Tar  Island  clay  can  be  described 
by  a  Singh-Mitchell  relationship.  The  strain  rate 
continually  decreases  with  time  and  is  proportional  to 
stress  level  for  a  given  time.  The  effects  of  transient 
consolidation  and  the  incremental  test  procedure  were 
assumed  to  be  negligible  in  fitting  the  creep  data  to  the 
Singh-Mitchell  relationship.  The  effects  of  stress  state, 
octahedral  stress  and  anisotropy  on  the  creep  parameters  are 
so  minor  that  they  cannot  be  detected  due  to  sample 

variations.  It  has  also  been  shown  that  Ca  is  not  a  constant 
and  decreases  with  time. 

The  effective  angle  of  shearing  resistance  obtained 
from  the  stress  controlled  tests  correlates  well  with  the 
previous  strain  controlled  tests  by  Hardy  Associates.  This 
demonstrates  that  0'  does  not  decrease  due  to  creep  strains 
over  the  duration  of  the  stress  controlled  tests. 


1  3 


6.  Analysis  of  Lateral  Creep  Movements  in  the  Tar  Island 

Dyke  Clay  Foundation 


6 . 1  Introduction 

The  lateral  movements  of  the  clay  foundation  of  Tar 
Island  Dyke  are  examined  in  this  chapter.  As  discussed  in 
Chapter  2  it  is  considered  that  the  time  dependent  movements 
displayed  by  the  inclinometers  are  a  consequence  of  creep 
rather  than  of  lateral  consolidation.  The  lateral  creep 
movements  of  the  clay  foundation  at  Station  65+00  are  shown 
in  Figure  6.1. 

A  simple  analysis  of  the  lateral  movements  of  one 
inclinometer  is  presented  in  this  chapter.  The  analysis 
treats  the  boundary  conditions  in  the  immediate  vicinity  of 
the  inclinometer  as  uniform  and  quasi-static.  The  form  of 
the  creep  movements  is  similiar  for  all  inclinometers.  A 
more  rigorous  creep  analysis  which  would  treat  the 
inclinometer  movements  collectively  requires  the  use  of  the 
finite  element  method  because  of  the  complex  geometry  and 
boundary  conditions. 

The  creep  problem  must  be  considered  as  drained  since 
pore  pressures,  where  they  are  in  excess  of  hydrostatic, 
have  dissipated  significantly  over  the  duration  of  the  creep 
movements  analysed.  It  will  be  assumed  as  per  Kavazanjian 
and  Mitchell  (1980)  that  the  creep  strain  tensor  can  be 
decomposed  into  a  volumetric  and  deviator ic  component.  There 
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Figure  6.1  Lateral  creep  movements  of  the  Tar  Island  Dvke 
clay  foundation,  Station  65+00 
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is  some  justification  for  this  assumption  in  that  the  form 
of  the  shear  strain  versus  depth  profiles  are  similar  for 
regions  of  the  foundation  that  have  different  rates  of 
excess  pore  pressure  dissipation. 

Lateral  or  horizontal  strains  with  time  have  been 
negligible ,  in  the  order  of  0.04  %,  so  volume  change  of  the 
clay  foundation  must  be  accommodated  primarily  in  the 
vertical  direction.  Vertical  displacement  readings  for  the 
clay  foundation  have  not  been  recorded  for  Tar  Island  Dyke 
so  the  volumetric  component  of  the  creep  strain  tensor 
cannot  be  analysed.  However,  total  volumetric  strains  from 
the  the  beginning  of  construction  can  be  calculated  from 
water  content  changes  as  determined  from  drilling. Using  this 
method  volumetric  strains  in  excess  of  20  %  have  occurred. 
Only  the  deviatoric  component  of  the  creep  strain  tensor 
will  be  analysed  here. 

It  will  be  assumed  for  the  purpose  of  this  analysis 
that  the  creep  deformations  can  be  described  by  an  equation 
of  the  Singh-Mi tchell  type.  One  justification  for  this 
assumption  is  that  the  laboratory  creep  of  the  clay 
described  in  Chapters  4  and  5  followed  such  a  relationship. 

The  analysis  then  is  concerned  with  defining  the 
various  parameters  in  the  Singh-Mi tchell  equation.  In  order 
to  do  this  it  is  first  necessary  to  determine  the  stress 
state  in  the  foundation  and  cast  the  Singh-Mi tchell  equation 
in  a  form  applicable  to  this  stress  state. 
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6.2  State  of  Stress  in  Foundation  Clay 
6.2.1  Introduction 

Tar  Island  Dyke  has  been  constructed  in  stages  from 
1965  and  by  1980  had  reached  an  elevation  of  327  m.  It  is 
intended  to  ultimately  raise  the  dyke  to  335  m.  The  dyke 
crest  elevation  and  pond  level  elevation  rise  with  time  are 
shown  in  Figure  6.2.  The  effective  stresses  in  the 
foundation  clay  with  time  are  a  function  of  the  increase  in 
dyke  load  with  time  and  the  dissipation  of  excess  pore 
pressure  with  time.  Since  it  is  considered  that  effective 
stress  governs  the  creep  deformation  of  the  soil  it  is 
necessary  to  determine  these  effective  stresses  for  the 
period  of  the  inclinometer  readings.  Two  approaches,  both  of 
which  use  elastic  constitutive  relationships,  to  the 
solution  of  this  problem  will  be  described  here.  The  use  of 
the  finite  element  method  is  necessary  in  both  approaches 
because  of  the  complex  boundary  conditions  and  incremental 
nature  of  construction. 

One  approach  to  this  problem  is  to  determine  the 
effective  stresses  directly.  Variations  in  this  approach 
have  been  outlined  by  Eisenstein  (1974),  Law  (1975),  Nobari 
and  Duncan  (1972)  and  Lee  and  Idriss  (1975).  An  advantage  to 
this  'effective  stress'  approach  is  that  material  properties 
from  drained  or  undrained  laboratory  shear  tests  with  pore 
pressure  measurements  can  be  used  directly  in  the  finite 
element  analysis.  This  procedure  is  necessary  if 
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Figure  6.2  Dyke  crest  and  pond  level  rise  with  time 
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displacements  other  than  those  arising  from  creep  are  to  be 
predicted  accurately.  Since  only  stresses  are  required  for 
this  part  of  the  analysis  a  simpler  approach  was  adopted. 

This  approach  employs  the  finite  element  method  to 
determine  the  total  stresses.  The  effective  normal  stresses 
are  then  calculated  by  subtracting  known  pore  pressures 
measured  in  the  field  from  the  total  stresses.  This  method 
is  accurate  for  the  determination  of  the  effective  vertical 
stress  in  the  clay  but  is  less  accurate  for  the 
determination  of  the  shear  stress  and  the  horizontal 
effective  stress  since  both  are  highly  dependent  on  the 
ratio  of  the  material  moduli  and  the  value  of  Poisson’s 
ratio.  The  next  section  describes  the  finite  element  stress 
analysis  and  the  approximations  used  to  overcome  the 
difficulties  of  selecting  the  material  parameters. 

6.2.2  Total  Stress  Analysis 
6 . 2 . 2 . 1  General 

The  finite  element  program,  NLCP ( nonl inear  construction 
program),  which  was  developed  by  Simmons  (1981), was  used  to 
obtain  the  total  stress  configuration  in  the  foundation 
clay.  The  program  is  reviewed  in  detail  by  Simmons  (op.cit .) 
and  only  a  brief  summary  of  the  features  of  the  program  used 
for  this  analysis  will  be  discussed  here. 
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The  program  employs  quadrilateral  elements.  Both  linear 
and  nonlinear  elastic  constitutive  relationships  can  be 
used.  Materials  can  be  modelled  as  either  isotropic  or 
cross-anisotropic.  Computer  graphics  support  programs  were 
available  for  plotting  displacements  and  mesh  outlines. 

6. 2. 2. 2  Boundary  Conditions 

The  finite  element  mesh  used  for  the  program  is  shown 
on  Figure  6.3.  Plane  strain  was  assumed  for  the  analysis. 
The  foundation  stratigraphy  is  an  idealized  version  of  that 
at  Station  65+00.  The  sand/limestone  interface  formed  the 
lower  boundary  of  the  grid  and  was  assigned  a  zero 
displacement  boundary  condition. 

A  total  unit  weight  of  1.97  Mg/m3  was  used  for  both  the 
compacted  zone  and  beach  zone  of  the  tailings  dyke.  Total 
unit  weights  for  all  materials  are  listed  in  Table  6.2.  The 
construction  boundary  between  these  tailings  zones  is  shown 
in  Figure  6.3.  The  densities  were  taken  to  be  equivalent  on 
the  basis  of  an  extensive  series  of  in  situ  density  tests  in 
Tar  Island  Dyke  reported  by  Mittal  and  Hardy  (1977),  which 
revealed  that  for  depths  greater  than  30  m  the  density  of 
the  beach  zone  was  equal  to  that  of  the  compacted  zone. 

The  upstream  boundary  of  the  tailings  dyke  was  modelled 
as  a  vertical  surface  with  two  degrees  of  freedom  and  a 
supporting  thrust.  The  magnitude  of  the  thrust  was 
calculated  on  the  basis  of  a  heavy  liquid  with  a  unit  weight 
equal  to  that  of  the  beach  zone.  It  must  be  emphasized  that 
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Figure  6.3  Tar  Island  Dyke  stress  analysis,  finite  element 
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the  thrust  does  not  add  to  the  shear  stresses  in  the  clay 
foundation  but  simply  supports  the  assumed  vertical  upstream 
boundary.  The  actual  configuration  of  the  upstream  slope  was 
not  used  because  stresses  and  displacements  in  this  portion 
of  the  dyke  were  not  relevant  to  the  analysis.  As  a  check  on 
this  assumption  a  zero  x-displacement  upstream  boundary 
condition  was  used  in  one  analysis  instead  of  the  stress 
boundary  condition.  It  was  found  that  the  stresses  in  the 
region  of  interest  were  not  affected  by  the  boundary 
condition  change. 

Dyke  loads  were  incrementally  applied  to  the  foundation 
in  the  finite  element  analysis  according  to  an  idealized 
version  of  the  actual  dyke  construction  history  as  portrayed 
in  G . C . 0 . S . ( Gr ea t  Canadian  Oil  Sands,  now  Suncor)  Drawing 
No.  TL-PL-009 ' E ' ,  entitled  Construction  History.  Six 
increments  were  used  and  their  timing  and  configuration  are 
shown  on  Figures  6.2  and  6.3.  The  finite  element  program, 
NLCP,  uses  the  nodal  displacements  and  stresses  calculated 
in  the  lower  increments  as  input  to  the  next  increment.  The 
total  stresses  in  the  foundation  before  dyke  construction 
were  calculated  according  to: 

6  v  =  »  D  6.1 

where  =  total  vertical  stress 


s$7  =  total  unit  weight 

D  =  depth 


. 
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and: 

6h  =  6v  ( I-  (^y^si n<fi‘)  6.2 

where  G)d  =  total  horizontal  stress 
^6  =  submerged  unit  weight 

0  -  effective  angle  of  shearing  resistance 

The  incremental  analysis  provided  an  estimate  of  the 
stresses  in  the  foundation  at  all  stages  of  construction. 

6. 2. 2. 3  Parametric  Study 

Two  types  of  elastic  analyses  were  performed  to 

determine  the  foundation  stresses.  The  first  analysis  was 
undertaken  using  a  linear  elastic  constitutive  relationship 
and  the  second  analysis  used  a  nonlinear  elastic 

constitutive  relationship.  A  constant  Poisson's  ratio  of 
0.49  was  used  in  all  analyses  except  for  one  linear  elastic 

analysis  where  a  value  of  0.35  was  used.  Further  parameteric 

analysis  should  treat  the  value  of  Poisson's  ratio  in  the 
sand  as  less  than  0.49  . 

A  summary  of  the  material  parameters  used  is  presented 
in  Tables  6.1  and  6.2.  One  way  of  comparing  the  results  of 
the  different  analysis  is  to  plot  the  horizontal  shear 
stress  in  the  clay  layer  just  beneath  the  tailings  sand/clay 
interface  along  the  downstream  slope.  Such  a  plot  is  shown 
in  Figure  6.4.  The  shear  stresses  are  those  at  the  end  of 
the  sixth  increment.  It  was  found  from  the  linear  elastic 
analyses  that  the  shear  stresses  depend  on  the  ratio  of  the 
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Table  6.1  Summary  of 
elastic  analyses 


material  constants  used  for  linear 


Analysis 

Number 


Poisson's  Ratio 

P 


Modulus  Ratio 
Es/Ec 


0.49 


2 


0.49 


10 


3 


0.35 


10 


Table.  6.2  Summary  of  material  constants  used  for  nonlinear 
elastic  analyses 


Material 

Type 

Total  Unit 
Weight  (Kn/m3) 

K (  1  > 

n  (  1  ) 

0' 

Rf  (  1  > 

P 

Compacted 
Tailings  sand 

19.25 

2000 

.54 

36° 

0.91 

0.49 

Clay 

19.32 

350 

0 

24° 

0.8 

0.49 

(  i  ) 


Defined  in  Appendix  D 


) 
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Berm  Elevations 


Figure  6.4  Horizontal  shear  stress  in  the  upper  half  of  the 
foundation  clay  along  the  downstream  slope 
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Figure  6.5  Total  vertical  stress  with  time  in  the  upper  half 
of  the  foundation  clay  for  the  291  m  and  277  m  berms 
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Young's  modulus  of  the  sand  to  the  clay,  Es/Ec,  not  on  the 
magnitudes  of  the  respective  moduli.  This  is  not  true,  of 
course,  for  the  elastic  displacements.  Increasing  the  moduli 
ratio  from  5  to  10  decreased  the  shear  stresses  in  the  clay 
by  approximately  25  percent  as  illustrated  in  Figure  6.4. 

For  the  linear  elastic  analysis  the  moduli  ratio  was 
not  altered  between  increments.  Allowing  volume  change  in 
the  clay  by  decreasing  Poisson’s  ratio  did  not  affect  the 
horizontal  shear  stress  substantially  beneath  the  291  m  berm 
where  the  creep  analysis  will  be  performed.  It  did  however 
decrease  the  shear  stresses  towards  the  toe  of  the 
embankment.  It  is  considered  that  stresses  downstream  of  the 
264  m  berm  are  not  reliable  for  all  of  the  analysis  shown  in 
Figure  6.4. 

The  nonlinear  analysis  was  performed  to  investigate  the 
effect  of  changing  the  moduli  ratio  between  increments  while 
keeping  Poisson's  ratio  constant.  The  hyperbolic  parameters 
listed  in  Table  6.2  for  the  clay  are  from  the  drained  tests 
discussed  in  Chapter  5.  The  hyperbolic  parameters  for  the 
tailings  sand  are  based  on  the  properties  of  a  similar  sand 
reported  in  the  literature  by  Duncan  and  Chang  (1970).  It 
must  be  emphasized  that  the  writer  is  not  advocating  the  use 
of  drained  parameters  in  a  total  stress  analysis.  The 
drained  parameters  were  used  because  it  is  considered  that 
the  tangent  moduli  ratio  is  applicable,  not  the  absolute 
value  of  the  moduli.  The  horizontal  shear  stress  results  are 
presented  in  Figure  6.4  and  are  within  the  range  of  the 
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linear  elastic  analyses. 

Figure  6.5  illustrates  the  increase  in  vertical  total 
stress  with  time  for  the  foundation  clay  beneath  the  291  m 
and  277  m  berms.  It  is  evident  that  the  vertical  stress  is 
not  sensitive  to  the  selection  of  the  material  parameters. 

6.2.3  Rationale  for  Simple  Shear 

Computed  shear  stresses  with  time  in  the  foundation 
clay  below  the  291  m  and  277  m  berms  are  presented  in 
Figures  6.6  and  6.7  respectively.  The  maximum  shear  stress 
as  well  as  the  horizontal  shear  stress  are  shown  in  these 
figures.  It  is  apparent  that  a  Poisson’s  ratio  of  0.49 
forces  the  principal  stresses  to  rotate  approximately  45°, 
that  is,  the  horizontal  and  maximum  shear  stresses  are 
approximately  equal.  The  principal  stresses  for  a  Poisson's 
ratio  of  0.35  do  not  rotate  45°  but  do  rotate  more  than  15°. 

A  criterion  for  simple  shear  is  that  the  principal 
stresses  rotate  during  shear  (see  Edgers  et.al.,  1973  and 
Gale,  1981).  Because  the  stress  analyses  indicate  rotation 
of  principal  stresses  and  because  the  form  of  the  creep 
strains  from  the  inclinometers  are  similar  to  that 
determined  from  laboratory  simple  shear  tests  it  is 
considered  that  the  stress  state  applicable  to  the  field 
situation  is  simple  shear.  Shear  stresses  are  not  dependent 
upon  the  magnitude  of  the  pore  water  pressure  so  the  shear 
stresses  from  the  total  stress  analysis  can  be  used  in  the 
creep  analysis.  This  of  course  neglects  the  effect  of 
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Figure  6.6  Shear  stress  with  time  in  upper  half  of  clay 
foundation  below  291  m  berm 
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Figure  6.7  Shear  stress  with  time 
foundation  below  277  m  berm 
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changing  material  parameters  during 
stresses.  The  effective  vertical 
material  parameters. 


consolidation  on  shear 
stress  is  insensitive  to 


6.3  Formulation  of  In  Situ  Creep  Equation 
6.3.1  Introduction 

A  creep  equation  based  on  the  Singh-Mi tchell 
relationship  is  formulated  in  this  section  for  the  shear 
strain  versus  time  curves  from  inclinometer  S76-103.  The 
methodology  for  obtaining  shear  strains  from  inclinometer 
readings  is  discussed  in  Chapter  2.  The  resultant  shear 
strain  versus  depth  plots  for  a  number  of  inclinometers  is 
included  in  Appendix  A.  These  shear  strain  versus  depth 
plots  can  also  be  displayed  as  shear  strain  versus  time  for 
each  depth.  Such  a  plot  is  shown  for  inclinometer  S76-103  in 
Figures  6.8  and  6.9.  The  formulation  derived  here  will  be 
developed  in  terms  of  the  curves  shown  in  Figures  6.8  and 
6.9. 

Inclinometer  S76-103  was  chosen  for  the  analysis 
because  samples  for  the  laboratory  testing  program  were 
obtained  nearby  and,  furthermore,  because  of  the  proximity 
of  a  reliable  piezometer.  The  same  analytical  method  to  be 
presented  here  for  S76-103  can  also  be  applied  individually 
to  the  other  inclinometers. 


. 
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o 


o 


Figure  6.8  Shear  strain  versus  time , inclinometer  S76-103, 
above-upper  half  of  clay,  below-lower  half  of  clay 


Shear  Strain, 


163 


o 


Depth,  m 

B-B 

71.93M 

®-e 

71.32M 

An* 

70.71M 

-t-f 

70. 1 0M 

X-X 

69.49M 

68. 88M 

4-f 

68. 28M 

JHC 

67.67M 

2-Z 

67.06M 

CM 


CD 


CM 


Calibration  Date  =  November  2,  1976 


Time,  days 


moo 


Figure  6.9  Shear  strain  versus  time, 
basal  sand 
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In  Figure  6.8  it  is  apparent  that  the  shear  strains  at 
a  depth  of  56.08  m  do  not  have  the  same  form  as  the  other 
curves.  The  shear  strains  are  much  larger  than  at  other 
depths  and  appear  to  have  accelerated  with  time.  The  larger 
strains  occur  immediately  below  the  tailings  sand/clay 
contact  and  are  evident  in  the  other  inclinometer  plots  at 
the  same  stratigraphic  position.  The  Singh-Mi tchell  equation 
cannot  be  applied  to  these  strains. 

6.3.2  Mult iaxial  Form  of  Singh-Mi tchell  Equation 

The  Singh-Mitchell  equation  provides  a  relationship 
between  axial  strain  rate,  stress  level  and  time  for 
triaxial  test  results.  In  order  to  use  this  relationship  for 
the  field  problem  it  is  necessary  to  cast  it  in  a  multiaxial 
form  and  then  extract  the  equation  applicable  to  a  simple 
shear  stress  state. 

The  generalization  of  a  uniaxial  creep  relationship  to 
a  general  multiaxial  relationship  has  been  outlined  by 
Odqvist  (1966).  His  approach  is  based  on  some  similarities 
between  creep  and  plastic  deformation.  Emery  (1971)  applied 
the  work  of  Odqvist  ( op.cit .)  to  develop  a  multiaxial 
relationship  for  thawed  soils  as  has  Sego  (1980)  for  frozen 
soils.  The  Singh-Mitchell  equation  is  cast  in  a  general 
multiaxial  form  according  to  the  method  of  Odqvist  (1966)  in 
Appendix  D.  If  simple  shear  boundary  conditions  are  applied 
to  the  general  equation  the  relationship  between  horizontal 
shear  stress  and  shear  strain  rate  becomes: 


‘ 
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i  '  v  '  ^ 

'Sxy  =  shear  strain  rate 


6.3 


triaxial  stress  state  creep  constants 
=  horizontal  shear  stress 


It  is  apparent  that  Equation  6.3  is  of  the  same  form  as 
the  triaxial  creep  equation  presented  earlier  as  Equation 
3.2.  This  theoretical  observation  has  practical  verification 
in  the  literature  where  creep  equations  describing  the  axial 
strain  rate  from  triaxial  creep  tests  have  been  shown  to  be 
of  the  same  form  as  equations  describing  the  shear  strain 
rate  from  simple  shear  tests.  Examples  are  given  by  Wu  et.al 
(  1978)*,  Edgers  et.al  (  1973)  and  Sego  (  1980). 

To  develop  a  field  creep  equation  the  constants  in 
Equation  6.3  must  now  be  determined  from  the  shear  strain 
versus  time  curves  shown  in  Figures  6.8  and  6.9. 

6.3.3  Derivation  of  Stress  Level  Term 

Equation  6.3  predicts  that  the  shear  strain  rate  is  an 
exponential  function  of  the  horizontal  shear  stress.  The 
horizontal  shear  stress  can  be  normalized  to  stress  level  in 
Equation  6.3  by  making  some  assumptions  about  the  state  of 
stress  in  simple  shear.  This  is  done  in  Appendix  D  and  the 
resulting  normalized  creep  equation  is: 


6.4 


Where  A5  =  3/2  A 
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6\/  =  effective  vertical  stress 
0  -  effective  angle  of  shearing  resistance 

It  is  apparent  in  Figures  6.8  and  6.9  that  the  creep 
strain  rate,  i.e.  the  slope  of  the  shear  strain  versus  time 
curves,  decreases  with  increasing  depth  in  the  clay 
foundation.  It  is  considered  that  the  decrease  in  strain 
rate  is  a  consequence  of  the  decrease  in  stress  level  from 
the  top  to  the  bottom  of  the  clay  layer.  The  stress  level 
decreases  with  depth  because  of  the  pore  pressure  decrease 
with  depth  in  the  clay  horizon  as  shown  in  Figures  2.5  and 
2.6.  One  method  of  determining  the  parameter,  a,  is  to 
compare  the  strain  rates  at  two  depths  where  the  stress 
levels  are  known. 

Proceeding  in  this  fashion  the  shear  strain  rates  at 
two  different  depths__are  given  by: 

=  Ai 

0*0 

&.  =  Ai 

6.6 

where  A,  =  strain  rate  at  zero  deviator  stress  at 

t  =  t, 


The  time,  t,  ,  in  Equations  6.5  and  6.6  is  the  time  period 
from  the  effective  initiation  of  the  creep  of  the  clay  layer 
to  the  time  at  which  the  inclinometer  was  installed.  A 
procedure  for  estimating  tt  will  be  outlined  in  the  next 
section . 
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If  t  is  equal  to  t,  then  the  term  (t./t)™  reduces  to 
one  in  Equations  6.5  and  6.6.  The  shear  strain  rates,  f,  and 
y,,  become  the  strain  rates  at  t,  or,  in  other  words,  the 
strain  rates  at  the  moment  of  installation  of  the 

inclinometer.  If  Equation  6.5  is  divided  by  Equation  6.6  the 
result  is:  __  _ 

/e  6^7 

The  parameter  A,  cancels  if  it  is  assumed  that  the  clay 
has  identical  creep  properties  at  both  depths.  Solving  for  a 
yields : 

S  =  (Ini)  6>8 

estimate  of  a.  from  field  data  Equation  6.7  requires 
the  strain  rates  at  two  depths  for  which  the  stress  level  is 
known . 

A  sample  calculation  of  a  for  inclinometer  S76-103  is 
given  in  Appendix  D.  The  value  of  a  calculated  is  5.1.  This 
compares  favourably  with  the  average  value  of  ci  of  3  from 
the  laboratory  testing  program. 

6.3.4  Derivation  of  t,  and  m 

It  is  customary  when  presenting  the  Singh-Mi tchell 
equation  as  is  done  in  Equation  6.4  to  include  unit  time  as 
the  numerator  in  the  time  variable.  The  coefficient,  A,  is 
the  corresponding  strain  rate  at  zero  stress  level  at  unit 
time.  For  this  field  case  it  is  not  known  when  the  clay 
began  to  creep  under  its  present  stress  level  so  A  cannot  be 
determined  directly  at  unit  time.  It  is  therefore  convenient 


. 
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to  replace  the  unit  time  numerator  by  the  time  corresponding 
to  the  period  between  initiation  of  creep  at  this  stress 
level  and  the  installation  of  the  inclinometer.  Equation  6.4 
is  then  rewritten  as: 


=  At  6  6.9 

t,  =  duration  of  creep  period  to  the  time  at 
which  the  inclinometer  was  installed 
A,  =  shear  strain  rate  at  zero  deviator  stress 
at  time  t, 


The  parameter  A,  is  easily  determined  from  the  following 
since  the  shear  strain  rate  and  stress  level  are  known  at 
time  t ,  . 


A, 


6.10 


where  y, 


=  shear  strain  rate  at  t, 


A  procedure  for  determining  y,  and  stress  level  is  presented 
in  Appendix  4. 

Implicit  in  Equation  6.8  is  the  assumption  that  the 
creep  of  the  clay  under  the  current  stress  level  is 
independent  of  accumulated  creep  strains  under  previous 
stress  levels.  This  concept  has  been  discussed  in  Chapter  5 
for  incremental  laboratory  drained  creep  tests  and  is 
considered  to  be  valid  for  this  analogous  field  situation. 


By  making  this  assumption  it  is  possible  to  derive  a 
closed  form  solution  for  t,  that  is  independent  of  m.  Two 
expressions  are  derived  for  t,  in  Appendix  D.  One  is  valid 


. 
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for  m  equal  to  one,  the  other  is  more  general  and  valid  for 
m  of  any  magnitude.  Values  of  t,  were  determined  by  using 
strain  rates  calculated  from  the  polynomial  regression 
curves  shown  in  Figure  D.l.  From  the  expression  for  m  equal 
to  one  values  of  t,  were  obtained  ranging  from  500,000  to 
3,000,000  minutes.  The  general  expression  for  t,  did  not 
yield  reasonable  values  for  t, .  This  occurred  because  the 
general  expression,  although  mathematically  correct,  cannot 
account  for  the  random  component  of  the  field  strain  rates. 

To  further  delineate  the  value  of  t ,  a  purely  empirical 
approach  was  adopted.  Strains  were  calculated  for  values  of 
t,  ranging  from  100,000  to  3,000,000  minutes  for  the  initial 
strain  rates  and  stress  levels  for  depths  of  67.06  m  and 
62.18  m  in  inclinometer  S76-103.  The  magnitude  of  m  was 
varied  between  0.7  and  1.3  for  each  value  of  t, .  This  range 
in  m  was  used  because  it  represents  all  reasonable  values 
for  soils  according  to  Mitchell  (1976).  Shear  strains  were 
calculated  from  one  of  the  following  equations.  For  m  not 
equal  to  one  strains  are  given  by: 

Y  =  Y.+  6.1, 

i  — m  ^  ' 

where  y,  =  accumulated  strain  at  t, 

At  =  net  time  past  t, 

For  m  equal  to  one  strains  are  given  by: 

Y  =  Yi  +  Ai  (In  6.12 

Equations  6.10  and  6.11  are  obtained  by  integrating  Equation 
6.8  with  respect  to  time.  It  was  assumed  in  the  integration 
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of  Equation  6.8  that  the  stress  level  is  independent  of 
time.  This  is  certainly  true  at  the  base  of  the  clay  layer 
where  the  pore  pressures  are  constant  but  not  strictly  true 
at  depth  62.18  m  where  piezometer  AP-75-18A  has  recorded 
substantial  excess  pore  pressure  dissipation  with  time. 
However  the  decrease  in  pore  pressure  only  reduces  the 
stress  level  by  6  percent  so  the  stress  level  was  assumed 
constant  for  the  purpose  of  integration. 

Time ,  t,,  was  then  selected  on  the  basis  that  it  must 
be  the  same  for  both  depths  and  that  it  must  result  in 
strains  that  correctly  match  those  recorded  for  both  depths. 
The  tfalue  of  t,  determined  in  this  way  was  approximately 
1,800,000  minutes  or  about  3.5  years.  The  shear  strain 
versus  time  curves  produced  using  this  value  of  t,  in 
Equations  6.11  and  6.12  are  presented  in  Figure  6.10. 

Inclinometer  S76-103  was  installed  in  late  1976.  A 
value  of  t,  of  3.5  years  predicts  that  the  creep  of  the  clay 
was  initiated  at  this  stress  level  in  mid-1973.  This  appears 
reasonable  since  this  time  corresponds  to  the  large 
increment  in  construction  in  1973  as  shown  in  Figure  6.2  and 
the  resulting  large  increase  in  shear  stress  shown  in  Figure 
6.6. 

As  shown  in  Figure  6.10  the  value  of  m  cannot  be 
determined  precisely  from  this  approach.  A  value  of  1.0  for 
m  gives  a  reasonable  fit  to  both  curves.  However  it  is 
possible  to  back  calculate  the  magnitude  of  y,  in  Equation 
6.11  given  a  t,  of  3.5  years.  For  values  of  m  greater  than 
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o 


Figure  6.10  Actual  and  predicted  field  shear  strain  versus 
time  curves 


o 


Figure  6.11  Field  shear  strain  versus  time  curves  and  curves 
predicted  from  laboratory  equation 
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1.0  high  values  of  y ,  are  obtained,  more  reasonable  values 

of  y,  are  obtained  if  m  is  taken  to  be  in  the  range  0.7  to 

1.0. 

It  is  of  interest  to  calculate  the  shear  strain  versus 
time  curves  predicted  by  the  laboratory  test  results  given  a 
value  of  t,  of  3.5  years.  This  is  shown  in  Figure  6.11.  It 
is  apparent  that  Equation  5.2  underestimates  the  field  shear 
strains.  The  discrepancy  is  considered  to  be  primarily  a 
result  of  the  difference  in  the  stress  level  parameter,^ 
' not  in  the  value  of  t, .  Size  effects  between  the  laboratory 
and  field  are  not  considered  to  be  the  seat  of  the 
discrepancy . 


6.4  Concluding  Remarks 

It  has  been  demonstrated  that  the  shear  strain  versus 
time  curves  for  two  depths  in  inclinometer  S76-103  follow 
similar  patterns  to  those  observed  in  laboratory  creep 
tests.  The  key  similarities  are: 

1.  Shear  strain  rate  decreases  with  time 

2.  Shear  strain  and  shear  strain  rate  are  exponential 
functions  of  stress  level. 

3.  A  significant  load  increment  ratio  during  stage 
construction  of  the  dyke  in  1973  caused  creep  strain 
rates  to  increase  and  then  decay  in  a  manner  independent 
of  the  previously  accumulated  creep  strains. 
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The  shear  strain  versus  time  curves  have  been  modelled 
with  a  S 1 ngh-Mi tchell  relationship  by  making  various 
assumptions  about  the  boundary  conditions  and  the  state  of 
stress  in  the  foundation  clay.  It  has  been  shown  that  the 
time, t,, can  be  approximated  by  assuming  that  the  creep  of 
the  clay  fits  a  Singh-Mi tchell  relationship. The  best 
estimate  of  t,  is  approximately  3.5  years.  Simple  procedures 
have  been  outlined  for  deriving  the  other  creep  parameters 
in  the  Singh-Mi tchell  relationship  from  the  field  data.  The 
stress  level  term,  a, is  approximately  5. The  damping  term  ,m, 
is  in  the  range  0.7  to  1.0. This  approximation  for  m  is 
tenous  since  it  is  based  on  a  back  calculation  which  assumes 
a  value  of  t,  . 
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7.  Conclusions  and  Recommendations 


7 . 1  Introduction 

This  thesis  has  examined  the  lateral  creep  deformation 
of  the  clay  foundation  of  Tar  Island  Dyke.  An  analytical 
technique  has  been  developed  which  enables  creep  parameters 
to  be  obtained  from  field  inclinometer  data.  Extensive  field 
data  concerning  foundation  geology  and  soil  properties 
, foundation  pore  pressures  and  dyke  construction  have  been 
synthesized.  A  series  of  oedometer  and  drained  triaxial 
creep  tests  have  been  performed  on  undisturbed  clay 
foundation  samples  at  different  sample  orientations,  cell 
pressures  and  stress  levels.  A  review  of  the  creep 
literature  on  soils  has  been  presented. 

Conclusions  and  recommendations  for  further  study 
arising  from  this  research  program  are  presented  in  the 
following  sections. 


7.2  Conclusions 

7.2.1  General 

Tar  Island  clay  is  a  normally  consolidated  Recent 
deposit  which  originated  by  infilling  of  an  abandoned 
channel  of  the  Athabasca  River.  The  clay  has  deformed 
vertically  and  horizontally  due  to  the  continued 
construction  of  Tar  Island  Dyke  which  began  in  1965.  At 
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present,  in  situ  piezometer  records  indicate  that  the  clay 
has  fully  consolidated  downstream  of  the  291  m  berm  and  that 

significant  excess  pore  pressures  are  dissipating  upstream 
of  this  berm. 

Inclinometers  installed  from  the  303  m  berm  downstream 
show  substantial  time  dependent  lateral  deformation. 
Horizontal  shear  strains  (in  some  instances)  are  in  excess 
of  5  %  over  a  four  year  period  while  horizontal  strains 
average  0.04  %  over  the  same  period.  These  small  horizontal 
strains  indicate  only  minor  compression  in  the  horizontal 
direction  and  are  indicative  of  a  creep  mechanism  rather 
than  lateral  consolidation  for  the  shear  strains.  A  further 
justification  for  postulated  creep  mechanism  is  that  lateral 
movements  resulting  from  horizontal  shear  straining  occur 
where  the  clay  is  fully  consolidated. 

Vertical  displacements  of  the  clay  layer  have  not  been 
monitored  directly  but  indirect  estimates  of  vertical 
displacements  can  be  obtained  from  water  content  changes  and 
changes  in  the  thickness  of  the  deposit  as  recorded  by 
drilling.  Indirect  estimates  such  as  these  reveal  that  the 
clay  at  some  locations  has  compressed  in  excess  of  3  m  over 
the  life  of  the  project. 

7.2.2  Laboratory  Creep  Tests 

Both  the  oedometer  and  drained  triaxial  test  results 
demonstrate  that  the  axial  creep  strain  rate  decreases 
according  to  a  power  law  with  time.  The  average  value  of  the 
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power  law  exponent,  m,  is  1.13  and  is  independent  of  the  two 
stress  states  used.  The  drained  triaxial  test  results 
demonstrate  that  m  is  independent  of  sample  orientation, 
stress  level  and  normal  stress. 

The  triaxial  test  results  show  that  the  axial  creep 
strain  rate  is  an  exponential  function  of  stress  level.  This 
result  together  with  the  power  law  dependence  of  strain  rate 
on  time  fulfill  the  requirements  of  a  Singh-Mi tchell  model 
for  creep  strains.  Transient  consolidation  effects  at  the 
beginning  of  the  stress  increments  disrupt  the’  decay  of 
creep  strain  rates  with  time  and  have  been  ignored  in 
fitting  the  Singh-Mi tchell  model  to  the  creep  data.  c> 

Shear  strength  as  determined  from  the  drained  triaxial 
creep  test  results  are  consistent  with  results  from  strain 
controlled  tests.  The  effective  angle  of  shearing 

resistance,  0,  is  26.5°  for  Tar  Island  Clay  and  is 

independent  of  sample  orientation. 

7.2.3  In  Situ  Lateral  Strains 

The  shear  strain  versus  time  curves  from  inclinometer 
S76-103  indicate  that  the  shear  strain  rate  decreases  with 
time  at  all  depths  in  the  clay  horizon  except  at  the 

clay-tailings  sand  contact.  At  this  location  the  shear 

strain  rate  has  accelerated  in  the  past  and  now  appears  to 
be  steady.  It  is  considered  that  the  clay  at  this  depth  is 
in  a  state  of  local  shear  failure.  Large  strains  are  evident 
in  the  clay  downstream  of  S76-103  at  the  same  stratigraphic 
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position  but  the  shear  strain  rates  here  are  decreasing  with 
time . 

The  shear  strain  rates  for  S76-103  are  a  function  of 
stress  level.  Preliminary  indications  are  that  the  shear 
strain  rate  is  an  exponential  function  of  stress  level.  More 
investigation  is  needed  to  confirm  this  hypothesis. 

Central  to  the  understanding  of  the  in  situ  creep 
behaviour  of  the  clay  is  the  estimation  of  the  time  period, 
t,»  that  the  clay  has  been  creeping  prior  to  the 
installation  of  the  inclinometer.  In  order  to  make  an 
estimate  from  this  field  data  of  t,  it  was  first  necessary 
to  assume  a  creep  model  for  the  soil.  It  was  assumed  that 
the  in  situ  creep  behaviour  could  be  described  by  a 
Singh-Mi tchell  relationship.  Furthermore,  it  was  assumed 
that  the  creep  strain  rates  at  the  current  stress  level  were 
independent  of  the  accumulated  creep  strains  at  the  previous 
stress  levels.  On  this  basis  t,  is  estimated  to  be 
approximately  3.5  years.  This  estimate  of  t,  appears  to  be 
reasonable  since  it  correlates  with  a  major  episode  of  dyke 
construction  during  1973. 


7.3  Recommendations  For  Further  Study 
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7.3.1  Drained  Laboratory  Creep  Experiments 

Further  research  is  required  to  delineate  the  influence 
of  stress  state  on  the  creep  behaviour  of  undisturbed  soils. 
A  complete  evaluation  of  the  influence  of  stress  state  would 
involve  oedometer,  triaxial,  plane  strain  and  simple  shear 
tests.  It  would  be  of  benefit  for  practioners  to  be  able  to 
estimate  parameters  for  the  Singh-Mi tchell  equation  for 
field  situations  that  involve  complex  stress  states  on  the 
basis  of  oedometer  tests. 

Of  particular  importance  to  the  estimation  of  the  in 
situ  creep  behaviour  of  clay  foundations  beneath  stage 
constructed  embankments  is  the  load  increment  ratio  required 
to  make  the  creep  strain  rates  independent  of  the  previously 
accumulated  creep  strains.  It  is  considered  that  this  load 
increment  ratio  is  a  function  of  stress  level  and  the 
duration  of  the  previous  stress  level.  It's  investigation  at 
high  stress  levels  is  complicated  by  the  small  load 
increments  required  to  cause  undrained  failure. 

Creep  tests  on  undisturbed  soils  of  long  duration 
similar  to  those  of  Lovenbury  (1969)  would  be  of  benefit  to 
the  long  term  prediction  of  creep  strains.  Of  paricular 
interest  are  the  strain  instabilities  at  extended  creep 
periods  observed  by  Lovenbury  (1969)  and  later  by  Tavenas 
et.al  (1978).  Strain  instabilities  have  not  been  displayed 
by  the  Tar  Island  Dyke  inclinometers. 


;  •-  ■ 


179 


7.3.2  Lateral  Creep  Strains  Beneath  Tar  Island  Dyke 

The  analytical  approach  described  in  Chapter  6  can  be 
applied  individually  to  the  other  Tar  Island  Dyke 
inclinometers.  The  magnitude  of  t ,  should  increase  in  the 
downstream  direction.  It  would  be  of  interest  to  install 
inclinometers  at  higher  berms  where  the  clay  has  been  loaded 
recently  by  dyke  construction  and,  of  most  interest,  because 
an  a  priori  assumption  would  not  have  to  be  made  about  the 
creep  model  which  best  fits  the  clay,  to  install  an 
inclinometer  in  a  position  in  the  clay  deposit  where  t, 
would  be  zero  at  some  time  in  the  future. 

The  next  logical  step  in  this  analysis  would  be  to 
treat  the  creep  displacements  of  the  clay  collectively  by 
solving  the  entire  boundary  value  problem.  This  would 
involve  the  use  of  a  finite  element  program  that  had  coupled 
facilities  for  properly  determining  the  effective  stresses 
in  the  clay  foundation. 
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Figure  A. 5  Inclinometer  S76- 1 03 , cal ibrated  November 
2, 1976, shear  strain 


I 

200 


E 

o 


c 

0) 

E 

8 

o 

a 

ca 

O 

ffl 

C 

o 

N 

•|_ 

o 

X 


Figure  A. 6  Inclinometer  S76- 1 03 , cal ibrated  November 
2, 1976, horizontal  displacement 
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Figure  A. 8  Inclinometer  S75- 1 04 , cal ibrated  September 
1 6 , 1 976 , hor i zontal  displacement 


Appendix  B 

Gram  Size  Curves, Scanning  Electron  Micrographs  and 

Oedometer  Test  Results 


B.1  Derivation  of  Equation  4.3 

For  a  given  change  in  vertical  total  stress  the 
effective  stresses  immediately  after  loading  arei 

£>2>£  -  63l-Aa6  1  td  1 


where 


6ii 

Gif 


=  6a  +  A  6i  (l-A')  b  .  2 

=effective  radial  stress  at  the  beginning 
of  the  increment 

=effective  radial  stress  at  the  end  of  the 
increment 

=effective  axial  stress  at  the  beginning  of 
the  increment 

=effective  axial  stress  at  the  end  of  the 
increment 


It  follows  that: 

sirup1  =  6/c  f  A6t  ~  A6i  A,-63c  +A6i  A 

6'\i  +  A6 i  -A6(  A  -t  63c-  A6i  A 

Solving  for  a6  ,  : 


AO,  =  0! ( 1-sin0 ’ )-o; ( 1+sin0 ’ )/( 1-2A) sin0 ’ -  1 )  B.4 
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Table  B.1  Oedometer  test  results,  sample  C2 


ncr . 

No. 

Stress 

( kPa ) 

Elapsed 

Time 

(min ) 

Void 

Ratio 

Change 

Ae 

Cv(  1  > 

(cm2/sec ) 

Ca 

m(  2  5 

Unit 
Strai 
Rate ( 3 
( %/min 

1 

31.3 

1700 

0(  4  > 

— 

— 

— 

2 

62.5 

2620 

.0055 

- 

.0030 

1.19 

.059 

3 

125.5 

5331 

.0157 

- 

.0010 

1.23 

.21 

4 

250.0 

4483 

.0236 

- 

.0016 

1.16 

.19 

5 

500.0 

2362 

.0315 

- 

.0021 

1.09 

.16 

6 

1000.0 

2937 

.0463 

.00065 

.0030 

1.11 

.30 

7 

1996.0 

401  1 

.0719 

.00045 

.0052 

1.12 

.69 

8 

4000 

1  1502 

.0928 

.00025 

.0060 

1.15 

.77 

( 1  Calculated  only  at  stresses  above  the  preconsolidation 
stress  where  possible 

( 2 Clope  of  the  logarithm  of  volumetric  strain  versus 
logarithm  of  time  curves 

(3)Intercept  of  logarithm  of  volumetric  strain  versus 
logarithm  of  time  at  t=1  min,  ignoring  consolidation  effects 
( 4  Cample  swelled  slightly 


. 

' 
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Table  B.2  Oedometer  test  results,  sample  C3 


ncr. 

No. 

Stress 

( kPa ) 

Elapsed 

Time 

(min) 

Void 

Ratio 

Change 

Ae 

Cv(  1  > 

(cm2/sec ) 

Ca 

m(  2  } 

Unit 
Strai 
Rate ( 3 
( %/min 

1 

39.30 

2691 

.0157 

— 

.0003 

— 

2 

78.0 

4555 

.  0202 

— 

.0015 

- 

- 

3 

158.9 

2482 

.  0296 

- 

.0020 

1.29 

0.71 

4 

314.0 

2676 

.  0375 

- 

.0025 

1.19 

0.41 

5 

624.0 

4093 

.  0493 

- 

.0033 

1.18 

0.56 

6 

1251 

4687 

.  0743 

.00015 

.0052 

1.15 

0.72 

7 

2502 

4180 

.0963 

.00010 

.0071 

1.17 

1.15 

8 

5027 

2865 

.0686 

.00011 

— 

1.16 

_ 

( 1 1 Calculated  where  possible 
preconsolidation  stress 
(2)Slope  of  the  logarithm  of 
logarithm  of  time  curves 
(  3  intercept  of  logarithm  of 
logarithm  of  time  at  t=1 
consolidation  effects 


at  stresses  above  the 
volumetric  strain  versus 


volumetric  strain  versus 
min,  ignoring  transient 


t 


• 

' 

Temperature, 
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Figure  B.1  Room  temperature  versus  time 
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Figure  B.3  Grain  size  analysi s , samples  Cl ,  C2  and  C3 
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Figure  B.4  Grain  size  analysis , samples  T1,  T2  and  T4 
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Figure  B.5  Grain  size  analys i s , samples  II,  12  and  13 
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Figure  B.6  Void  ratio  versus  time,  sample  C2 
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Figure  B.7  Void  ratio  versus  time,  sample  C3 
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Appendix  C 

Triaxial  Creep  Test  Results 
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Table  C. 1  Summary  of  incremental  creep  test  results 

sample  T1  ' 


Incr . 
No. 

Deviator 
Initial 
( kPa ) 

Stress 
Final 
( kPa ) 

Stress 

Level 

(%) 

Time 

(min ) 

Axial 

Strain 

(%) 

Cumulative 
Axial  Strain 
(%) ( 2  > 

1 

181 

181 

8.3 

20 

0.10 

0.10 

2 

326 

326 

15.0 

1875 

0.27 

0.37 

3 

515 

514 

23.6 

16742 

0.58 

0.95 

4 

681 

681 

31.2 

5636 

0.33 

1.28 

5 

892 

891 

40.9 

1339 

0.35 

1.63 

6 

1087 

1077 

49.9 

15754 

1.51 

3.  14 

7 

1263 

1259 

58.0 

5885 

0.56 

3.70 

8 

1467 

1467 

67.4 

27255 

1.57 

5.27 

9 

1530 

1530 

70.3 

1334 

0.06 

5.33 

10 

1595 

1592 

73.2 

5692 

0.33 

5.66 

1  1 

1649 

1645 

75.7 

3019 

0.27 

5.93 

12 

1695 

1688 

77.8 

4573 

0.60 

6.53 

13 

1720 

1718 

78.9 

1  104 

0.15 

6.68 

14 

1781 

1764 

81.8 

5918 

1.29 

7.97 

15 

1804 

1804 

82.8 

4224 

0.39 

8.36 

16 

1881 

1881 

86.4 

1517 

0.39 

8.75 

17 

1959 

1959 

90.0 

5572 

3.41 

12.16 

18 

1998 

1998 

91.8 

4149 

0.37 

12.53 

19 

2030 

2030 

93.2 

280 

0.07 

12.60 

(cont ’ d ) 
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Table  C . 1  (cont ’ d) 


Incr . 
No. 

Deviator 
Initial 
( kPa ) 

Stress 
Final 
( kPa ) 

Stress 

Level 

(%) 

Time 

(min ) 

Axial 

Strain 

(%) 

Cumulative 
Axial  Strain 
(%) ( 2  ’ 

20 

2049 

2047 

94.1 

350 

0.08 

12.68 

21 

2089 

2085 

95.9 

808 

0.26 

12.94 

22 

2126 

2104 

97.6 

1475 

1.70 

14.64 

23 

2143 

2085 

98.4 

3130 

3.42 

18.06 

24 

2116 

21145 

97.2 

270 

0.09 

18.15 

25 

2152 

2141 

98.8 

1020 

0.61 

18.76 

26 (  1  > 

2178 

2155 

100.0 

1430 

2.87 

21.63 

Total  duration  of  test=120371  mins,  O3 

=  1103'  kPa 

(  1  failure 

( 2 ’Strain  at  the  end  of  the  increment 


. 

' 
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Table  C.2  Summary  of  incremental  creep  test  results, 

sample  T2 


Incr 

No. 

.  Deviator 
Initial 
( kPa ) 

Stress 
Final 
( kPa ) 

Stress 

Level 

(%) 

Time 

(min ) 

Axial 

Strain 

(%) 

Cumulative 
Axial  Strain 
(%) ( 2  > 

1 

186 

186 

13.5 

1645 

0.14 

0.14 

2 

366 

365 

26.7 

17171 

0.48 

0.62 

3 

490 

490 

35.7 

2736 

0.14 

0.76 

4 

615 

614 

44.8 

2686 

0.44 

1.20 

5 

738 

740 

53.8 

17418 

0.86 

2.06 

6 

864 

863 

63.0 

2574 

0.37 

2.42 

7 

986 

970 

71.9 

1181 

1.18 

3.60 

8 

1055 

1082 

76.9 

20418 

1.89 

5.49 

9 

1  122 

1121 

81.7 

1187 

0.06 

5.55 

10 

1161 

1  157 

84.6 

4603 

0.45 

6.00 

1  1 

1  187 

1  186 

86.6 

1395 

0.21 

6.21 

12 

1216 

1206 

88.6 

5911 

1.05 

7.26 

13 

1243 

1238 

90.6 

4159 

0.55 

7.81 

14 

1312 

1299 

95.6 

1457 

1.25 

9.06 

1  5 ( 

1  >  1372 

1307 

100.0 

380 

5.30 

14.36 

Total 

duration  of 

test  = 

84921  mins, 

758  kPa 

( 1 5 Failure  did  not  occur,  sample  strained  vertically  5.3% 
( 2 )Strain  at  the  end  of  the  increment 


' 
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Table  C.3  Summary  of  incremental  creep  test  results, 

sample  II 


Incr 

No. 

.  Deviator 
Stress 
( kPa ) 

Stress 
Final 
( kPa ) 

Stress 

Level 

(%) 

Time 

(min ) 

Axial 

Strain 

(%) 

Cumulative 
Axial  Strain 
(%) ( 2  5 

1 

233 

233 

20.5 

505 

0.48 

0.48 

2 

413 

412 

36.3 

6800 

0.85 

1.33 

3 

594 

593 

52.2 

1280 

0.55 

1.88 

4 

775 

767 

68.1 

1431 

1.68 

3.56 

5 

949 

948 

83.4 

13189 

3.44 

6.99 

6 

1064 

1064 

93.9 

2850 

0.47 

7.47 

7 ( 

1  >  1138 

1131 

100.0 

212 

0.63 

8.10 

Total 

Duration  of 

Test  = 

26267  mins,  Oa  = 

743  kPa 

(  1  failure 

( 2  Strain  at  the  end  of  the  increment 


• 
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Table  C.4 


Summary  of  incremental  creep  test 
sample  12 


results , 


Incr 

No. 

.  Deviator 
Initial 
( kPa ) 

Stress 
Final 
( kPa ) 

Stress 

Level 

<%) 

Time 

(min ) 

Axial 

Strain 

(%) 

Cumulative 
Axial  Strain 
(%) ( 2 ’ 

1 

356 

356 

17.5 

1331 

0.47 

0.47 

2 

696 

691 

34.2 

7103 

1.74 

2.21 

3 

894 

892 

43.9 

1648 

0.42 

2.63 

4 

1079 

1073 

53.0 

1190 

0.96 

3.59 

5 

1290 

1288 

63.3 

10180 

2.53 

6.12 

6 

1531 

1525 

75.2 

1350 

0.67 

6.79 

-  7 

1648 

1640 

80.9 

310 

0.62 

7.41 

8 

1729 

1693 

84.9 

1380 

2.72 

10.13 

9 

1744 

1740 

85.6 

1070 

0.28 

10.41 

10 

1793 

1785 

88.0 

8845 

0.82 

11.23 

1  1 

'1821 

1820 

89.4 

1461 

0.06 

11.29 

12 

1895 

1892 

93.0 

1264 

0.13 

11.42 

13 

1966 

1963 

96.5 

340 

0.13 

11.55 

1  4  ( 

2)  2037 

2024 

100.0 

330 

0.72 

12.27 

Total 

Duration  of 

Test  = 

37802  mins,  Ol= 

1386  kPa 

( 1 ’Failure 

( 2 ’Strain  at  the  end  of  the  increment 


'6iv so;  -  iyi 
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Table  C.5  Summary  of  incremental  creep  test  results, 

sample  13 


Incr 

No. 

.  Deviator 
Initial 
( kPa ) 

Stress 

Final 

(kPa) 

Stress 

Level 

(%) 

Time 

(min ) 

Axial 

Strain 

(%) 

Cumulative 
Axial  Strain 
(%) ( 2  > 

1 

291 

291 

17.5 

1281 

0.45 

0.45 

2 

539 

539 

.  32.4 

7123 

0.07 

0.52 

3 

695 

693 

41.7 

1667 

0.44 

0.96 

4 

834 

832 

50.1 

1240 

0.42 

1  .38 

5 

1022 

1021 

61.4 

10110 

1.69 

3.07 

6 

1205 

1201 

72.4 

1318 

0.47 

3.54 

7 

1290 

1287 

77.5 

310 

0.34 

3.88 

8 

1350 

1339 

81.0 

1400 

1.10 

4.98 

9 

1402 

1397 

84.2 

9915 

1.41 

6.39 

10 

1466 

1464 

88.0 

1461 

0.14 

6.53 

1  1 

1525 

1522 

91.6 

1264 

0.25 

6.78 

12 

1588 

1583 

95.3 

330 

0.34 

7.12 

13 

1621 

1599 

97.4 

1095 

1.58 

8.70 

14 

1636 

1626 

98.3 

410 

0.78 

9.48 

1  5 ( 

15  1665 

1648 

100.0 

230 

1.18 

10.66 

Total 

Duration  of 

Test  = 

39154  mins,  6\= 

1090  kPa 

( 1 } Failure 

( 2 )Strain  at  the  end  of  the  increment 


9  a  •  .  i 
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Figure  C.1  Vertical  strain  and  strain  rate  versus  time, 
sample  II,  consolidation  pressure=365  kPa 
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Figure  C.2  Vertical  strain  and  strain  rate  versus  time, 
sample  II,  consolidation  pressure=743  kPa 
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Figure  C.3  Vertical  strain  and  strain  rate  versus  time, 
sample  12^  consolidation  pressure=334  kPa 
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versus  time, 


Figure  C.4  Vertical  strain  and  strain  rate 
sample  12,  consolidation  pressure=672  kPa 
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Figure  C.5  Vertical  strain  and  strain  rate  versus  time, 
sample  13,  consolidation  pressure=271  kPa 
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Figure  C.6  Vertical  strain  and  strain  rate  versus  time, 
sample  13,  consolidation  pressure=548  kPa 
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Figure  C.7  Vertical  strain  and  strain  rate  versus  time, 
sample  13,  consolidation  pressure= 1 090  kPa 
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Figure  C.8  Vertical  strain  versus  time, 


sample  T1 


Vertical  Strain,  % 
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Figure  C.9  Vertical  strain  versus  time, 


sample  T2 
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Figure  C.10  Vertical  strain  versus  time,  sample  II 
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Figure  C . 1  1  Vertical  strain  versus  time,  sample  12 
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Figure  C.12  Vertical  strain  versus  time, 


sample  13 
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Figure  C.14  Vertica 
12,  right-sample  II 
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Figure  C . 1 5 
T2-lower 


Stress 


strain  curves , sample 


T 1 -upper , sample 
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Figure  C.16  Stress  strain  curve,  sample  t4 


AV/V,  %  Deviator  Stress  (  a-j  —  (t3 ),  kPa 


Figure  C.17  Stress  strain  curve,  sample  II 


Figure  C.18  Stress  strain  curve,  sample  12 
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Plate  C . 1 


Sample  T1 -upper, 


sample  T2-lower, 


after  failure 
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Plate  C . 2  Sample  T4-upper, 


sample  11-lower,  after  failure 


I 
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Plate  C . 3  Sample  I2-upperf  sample  13-lower,  after  failure 


Appendix  D 

Derivation  of  Formulas  for  Analysis 


D.1  Definition  of  Nonlinear  Elastic  Material  Parameters 

Nonlinear  elastic  constitutive  relationships  in  soil 
mechanics  can  be  derived  from  the  observation  that  stress 
strain  curves  can  be  described  by  the  equation  for  a 
hyperbola.  This  was  first  reported  by  Kondner  (1963)  and 
later  popularized  by  Duncan  and  Chang  (1970).  Reference 
should  be  made  to  this  latter  publication  for  a  more 
detailed  explanation  than  is  given  here. 

A  stress  strain  curve  can  be  represented  by  an  equation 
for  a  hyperbola  of  the  form: 


D.  1 


where 


£  =  strain 

E  =  initial  tangent  modulus 

/-6b)  =  asymptotic  value  of  deviator  s 


stress  as 


represented  by  hyperbolic  equation 


The  compressive  strength  of  the  soil  is 
( 6 , -d3 ) ul t  by: 


related  to 


D.  2 
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where  Rf  =  failure  ratio 


The  initial  tangent  modulus  varies  with  confining  stress 
according  to: 


El 


=  Kpa(6^fpa) 


h 


D .  3 


where  K 
n 

pa 


=  modulus  number 
=  modulus  exponent 

=  atmospheric  pressure  in  chosen  units 


The  Mohr  Coulomb  failure  criterion  is  given  by: 


(61-65)1--  Zc  COS^(±Z63  Slhjfi1 

r-s/h#' 


By  differentiating  Equation  D.1  with  respect  to  strain  and 
substituting  Equations  D.2,  D.3  and  D.4  the  tangent  modulus 
is  given  by:  Z  n 

£t  =//-  k, Pd[ / '<5a)  D  5 
/  2c  costf+zees/nfSj  (■ 

Poisson’s  ratio  can  be  described  by  an  equation  which 
accounts  for  its  dependence  on  normal  stress  and  shear 
stress  (for  contractant  materials  only)  as  per  Duncan  et .31 
(1978).  Poisson's  ratio  was  taken  as  a  constant  since  this 
analysis  was  assumed  undrained. 


D.2  Transformation  of  Singh-Mi tchell  Equation  To  a 
Multiaxial  State  of  Stress 

The  following  formulation  follows  that  of  Odqvist 
(1966).  Odqvist  begins  his  formulation  by  drawing  from  the 
theory  of  plasticity.  He  defines  the  rate  of  dissipation  of 


.3 

* 


’  ^ 


247 


energy  in  a  creep  test  by: 

5  U/y  -  &Sij  Vcj  d.6 

where  SW y  =  change,  in  rate  of  dissipation  of  energy 
5Sy  =  increment  of  stress  deviator  tensor 
=  strain  rate  tensor 

Equation  D.6  provides  a  relationship  between  the  stress 
deviator  tensor  and  strain  rate  tensor.  He  continues  by 
stating  that  W  is  a  flow  potential  and  dependent  on  some 
scalar  function.  He  calls  the  scalar  function  the 
'effective'  stress  and  defines  it  in  terms  of  the  second 


stress  invariant  as: 

6e  5=1  3/^  Scj  Sij 

-  6n  +  €£2.  +  633  "  6//  62 z  ~  622-  633  -633  6// 


D .  7 


y  3^6/3.  y-  623  +-  63/) 


Equation  D.6  can  now  be  manipulated  as  follows: 

w..  -  dw  *  _  _W  J_  *  d  6e  _  5£/'D.8 

d  <i6e  ciStj  d6e  26e  dScj  d6e  6e 

The  preceding  assumes  coaxiality  of  the  tensors  of  strain 

rate  and  stress.  It  also  assumes  no  volume  change  and 
i sotropy . 

For  the  case  of  uniaxial  stress: 

5//  -  2  611/3 

6e- 

Sy$—  “'6^/3 

Now  Equation  D.8  will  be  manipulated  so  that  it  reduces  to 
the  Singh-Mi tchell  equation  for  a  uniaxial  stress  state. 

V/i  =  3A  dU/  Z_  &1L  =  4<S  d . 9 

3  6(1 


This  reduces  to: 


. 
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o(  6<S 


D.  10 


aw  ^  Ae^n/J. Y 

Substituting  D.10  into  D.8  we  have: 

*r  H 

Equation  D.11  is  a  general  multiaxial  equation  relating  each 


Vej  -  %  4  e*6e/< 


D.  1  1 


term  in  the  strain  rate  tensor  to  each  term  in  the  stress 
deviator  tensor.  The  relationship  desired  for  the  analysis 
is  one  between  horizontal  shear  strain  and  horizontal  shear 
stress  for  a  simple  shear  stress  state.  It  will  be  assumed 
that  a  pure  shear  stress  state  is  representative  of  a  simple 
shear  stress  state  at  zero  normal  stress.  The  assumptions  of 


pure  shear  result  in  the  following: 

6i-th  ,6z-o ,  6^  =  'b, 

5//  '-533  “  &  )  5x-  O 

The  effective  stress  as  given  by  Equation  D.2  reduces  to  the 


D.  12 


following : 


6>e  =2  'th 


D.  13 


If  this  is  now  substituted  into  Equation  D.7  the  following 


results : 


Reducing : 


"  u 


C/e)  ^ 


D.  14 


D.  15 


D . 3  Definition  of  Stress  Level  For  Simple  Shear  Stress 
Conditions 

In  order  to  define  the  stress  level  in  simple  shear  it 
is  necessary  to  make  some  assumptions  about  the  state  of 
stress  (see  Ladd  and  Edgers,  1972  and  Gale,  1981).  One 
assumption  is  that  the  horizontal  plane  is  the  failure 
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plane,  another  is  that  the  horizontal  plane  is  the  plane  of 
maximum  shear  stress.  These  two  assumptions  converge  as  the 
effective  angle  of  shearing  resistance  decreases.  For  an 
angle  of  shearing  resistance  of  24°  of  the  foundation  clay, 
the  maximum  difference  between  these  assumptons  in  defining 
the  shear  stress  level  is  less  than  10  percent.  It  will  be 
assumed  that  the  horizontal  plane  is  the  failure  plane. 

The  tangent  of  the  angle  of  shearing  resistance  for  a 
simple  shear  test  is  then  defined  as: 

iqrnft1-  'tmqi / 6V  d.16 


where  ^  x  =  shear  stress  on  horizontal  plane  at 

failure 

6V  =  effective  normal  stress  on  horizontal 
plane  at  failure 


The  shear  stress  level  is  defined  by: 

D  —  'hlo /  'fcmCDC.  D.17 


where  D 


=  shear  stress  level 

=  shear  stress  on  the  horizontal  plane 


The  stress  level  term  in  Equation  6.3  then  becomes: 

D  =  'kh/ 

It  is  of  interest  to  compare  the  stress  level 
from  a  triaxial  and  a  simple  shear  test.  In  a  triaxial 

Scr-  oC  i/2  (6,-63)  f 


D.  18 
terms 
test : 
D.  19 


In  a  simple  shear  test:  ^ 

^5  -  oC  6v/  hr><f> 
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D .  20 

If  the  approximation  is  made  that  <5vtan0  ’<=  1  /2  (6 ,  -63 )  f  then: 

=  d-T  0,21 

D.4  Sample  Calculation  of  a 

The  shear  strain  versus  time  curves  for  inclinometer 
S76-103  will  be  used  to  determine  a  value  of  a  from  Equation 
6.7.  Two  depths  in  the  clay  layer  are  required  for  the 
calculation  of  a  at  which  the  pore  pressures , hence  the 
stress  levels,  are  known.  The  obvious  choices  for 
inclinometer  S76-103  are  the  base  of  the  clay  layer  which  is 
at  a  depth  of  .67.07  m  and  the  depth  corresponding  to  the 
piezometer  AP-75-18A  of  62.18  m. 

The  stress  level  at  the  base  of  the  clay  is  given  by: 

D/  -  t/7  / 6v  fantfi' 

-  124  kP a  / (1323  kPq  -  IIOK~Pq)  fan  2iv 

-  0-40 

The  stress  level  at  62.18  m  is: 

2)^  =  J-75 H2o tPa.  - 35o  k?®)  0-4° ^ 

~  0‘5/ 

In  order  to  use  Equation  6.7  the  strain  rates  at  t,  are 
required.  This  was  done  by  passing  the  best  fit  second  order 
polynomial  through  the  shear  strain  versus  time  curves  and 
then  calculating  the  shear  strains  at  t, .  Inserting  these 
values  into  Equation  6.7  results  in: 
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—7. 


^  =  /n_  (5$ * %/^ )  -  In  (/Q4-xi0~^/ol*i)n)  D_ 

- 4o-.5! 

-  5.1 


22 


D.5  Procedure  For  Determining  t, 

The  procedure  detailed  here  for  the  determination  of  t, 
is  based  on  the  assumption  that  m  is  constant,  i.e.  the 
slope  of  the  natural  logarithm  of  strain  rate  versus  the 
natural  logarithm  of  time  is  a  straight  line.  The  equation 
of  such  a  line  is  given  by: 

In  Y%  ■=  JnX -  /v  In  ^  j  d.  23 

where  y,  =  shear  strain  rate  at  time  t, 

■y  2  =  shear  strain  rate  at  time  At2+t, 


If  m  is  equal  to  1  then  t,  can  be  found  fr 


om: 


ti  =  Atz 


-7 


D .  24 


D .  25 


If  m  is  not  equal  to  one  then  two  equations  are  needed  to 
solve  for  TI.  Another  expression  similar  to  D.23  can  be 
written  for  the  same  line  as: 

In  Yi=  In  Y, -m/n  ( 

Solving  for  m  in  Equation  D.25  and  substituting  this 

expression  into  Equation  D.23  yields: 

In  Yj,  -  In  Y,  __  Jn_  ((At$ +ti)/ti) 

lni?.-lnYi  In  ((Atzrt/)/t/) 

The  left  hand  side  of  Equation  D.26  is  a  constant  and  will 


D .  26 


D .  27 


be  called  P.  Rearranging  D.26  yields: 

/  Atl+tl  \P  _  At?,+ti 

(  t,  )  tt 

A  solution  for  t,  is  a  root  of  Equation  D.27.  The  roots  of 

Equation  D.27  can  be  obtained  by  using  a  numerical  technique 

such  as  Newton-Raphson .  To  proceed  in  this  way  Equation  D.27 


. 

£  .  lo  no  t 


■,  u 
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Figure  D.1  Regression  analysis  of  shear  strain  versus 
at  depths  67.06  m  (upper)  and  62.18  m  (lower) 
inclinometer  S76-103 


t  ime 
for 
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must  be  rearranged  as  follows: 

-p 

(X+')-/Cx-/=  o  D.28 

where  X  =  At2/t, 

K  =  At  3 /At  2 
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